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A MITOCHONDRIAL FACTOR PRODUCING 
“LOOSE-COUPLING” OF RESPIRATION* 


By LEMar F. REMMErRTT AND ALBERT L. LEHNINGER 
DEPARTMENT OF PHYSIOLOGICAL CHEMISTRY, JOHNS HOPKINS SCHOOL OF MEDICINE, BALTIMORE 


Communicated November 10, 1958 


Carefully prepared liver or kidney mitochondria yield high rates of respiration 
and maximum P:QO ratios in the presence of adenosine diphosphate (ADP) and 
phosphate but respire only very slowly in the absence of ADP, as was first shown 
by Lardy and Wellman.' Such ‘‘tight-coupling” of respiration has since been 
observed in many laboratories and has led to the concept? that the ADP level is 
an important determinant of the rate of respiration. On the other hand, some 
types of mitochondrial preparations, although capable of yielding high P:O ratios 
in the presence of ADP, do not require the presence of ADP for high rates of respira- 
tion. An example is provided by the fly sarcosomes described by Sacktor and 
Cochran.’ Phosphorylating respiration of the latter type may be designated as 
“loosely coupled.” There is very little information available concerning the 
factors responsible for ‘tightness’? or “looseness” of respiratory coupling, the 
mechanism of these effects, or their physiological significance. 

In this communication we wish to report that a soluble protein-containing fraction 
obtained from rat liver mitochrondria is capable of fully releasing the respiration 
of mitochondrial subfragments from its dependence on the presence of ADP, 
without lowering the P:O ratio greatly when ADP is available. This mitochondrial 
factor, designated here as “R factor’ for its respiration-releasing activity, is thus 
capable of converting a “tightly coupled” phosphorylating respiration to a ‘loosely 
coupled” system. 

The mitochondrial fractions containing the R-factor activity were prepared as 
follows: Freshly isolated rat liver mitochondria, following washing with 0.25 M 
sucrose, were suspended in 0.05 .W sucrose—0.02 1 KCl and subjected to sonic 
disruption at 9 kilocycles at 0° for 60 minutes. The sonicate was adjusted to pH 
7.4, clarified by sedimentation at 115,000 X g for 30 minutes, dialyzed against 
several changes of 0.03 1 KCI for 40 hours, again adjusted to pH 7.4, and clarified 
with a high-speed centrifugation. This dialyzed extract, or a protein fraction 
precipitating from it at pH 5.5, was used in the experiments described. The 
action of the R factor was assayed on membrane fragments of rat liver mitochondria 
capable of phosphorylating electron transport;‘~* under the conditions used in 
these experiments the oxidation of 6-hydroxybutyrate (BOH) by oxygen is ordi- 
narily stimulated between 60 and 400 per cent by the presence of an excess of ADP.? 


l 
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Figure 1 shows that addition of increments of the R-factor preparation to a 
suspension of mitochondrial fragments oxidizing BOH in the absence of ADP 
am produces increasing rates of respira- 
ADP_=<=_ONP tion, until a plateau is reached, at 
which further addition of R= faetor 


mM 
oO 
Oo 


produces no further increase in rate 
of oxidation of BOH. This plateau 
coincides closely with the maximum 
rate of respiration which is given by 
the presence ot an excess of ADP as 
phosphate acceptor; it coincides also 


A. My MOLES BOH 


with the maximum rate of respiration 

given by the presence of concentra- 

L i tions of  2,4-dinitrophenol (DNP) 

2 04 which uncouple phosphorylation com- 

R FACTOR (ML.) pletely. t-factor preparations thus 
Fic. 1.—The test system contained 0.015 are capable of full release of respiration 


phosphate, pH 6.5, 0.01 M DL-BOH, and mito- from dependence on ADP when added 
chondrial fragments (286 ug of N) in a total ; . E 
volume of 2.0 ml., incubated 17 minutes at 22°. to the submitochondrial preparations. 
The R-factor preparation contained 2.56 mg. of The R-factor preparations, when 
N per ml. The data are corrected for the i G 
amount of BOH oxidized by the mitochondrial tested in the system alone in the 
fragments alone; respiratory release by R factor gbsenee of the mitochondrial sub- 
is compared with that given by ADP (0.0024 17) fi a aes 
or DNP (1 X 10-5 M fragments, show no ability to oxidize 
BOH. 


The stimulation of respiration caused by R factor is qualitatively different from 








RESPIRATORY STIMULATION 


that caused by DNP. The data summarized in Table 1 compare the action of 


TABLE 1 
Errects or R Facror aNd DNP ON RESPIRATION AND PHOSPHORYLATION * 


Respiration 
mymoles BOH 
Addition ADP Oxidized) 


R factor 

R factor 

DNP, 3 X 10°°M 

DNP, 3 X 10-5M | 2: 0.02 

* The test system contained 0.015 M phosphate (7.2  105¢.p.m. of P82), pH. 6.5, 0.01 M DL-BOH, and mito 

chondrial fragments (220 ug. of N) in a total volume of 2.0 ml. Where added, the ADP concentration was 0.0024 
M; the R factor contained 384 we. of N. The reaction period was 15 minutes at 22°. 
such concentrations of the R factor and of DNP as just produce full release of 
respiration, tested in the presence and absence of ADP. It is seen that both R 
factor and DNP stimulate respiration of the ADP-less system to the level given by 
ADP alone. However, in the presence of excess ADP, the R factor causes very 
little uncoupling of phosphorylation, whereas DNP uncouples nearly completely. 
It. is therefore evident that the addition of the R factor has converted the relatively 
“tightly coupled” system of Experiment 1, which requires ADP for a maximum 


rate of respiration, into the “loosely coupled” phosphorylating system of Experi- 
ment 2, which has no dependence on ADP for respiration. The P:O ratios in the 
“tightly coupled” and “loosely coupled” systems are approximately the same. 
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On the other hand, DNP releases respiration but produces an “uncoupled” system. 

The respiration released by the presence of the R factor is completely inhibited by 
antimycin A in low concentrations, as is the respiration in the presence of ADP.‘ 
Tests of the site of action of the factor in the respiratory chain indicate that it 
stimulates electron transfer not only in the span BOH to cytochrome c¢* but also 
in the span cytochrome ¢ to oxygen, tested with ascorbate as non-enzymatic 
reductant of cytochrome c*. It therefore appears likely that the factor can act at 
all three coupling sites in the respiratory cham. The R factor also stimulates the 
oxidation of BOH by freshly prepared ntact rat liver mitochondria in the absence 
of ADP. On the other hand, the respiration of non-phosphorylating mitochondrial 
fragments similar to Keilin-Hartree preparations is not stimulated by the R factor. 
The ability to release respiration in the submitochondrial fragments in the absence 
of ADP, without severe uncoupling when ADP is present, is relatively unique; it is 
not shown by DNP, Dicumarol, gramicidin, ethylenediaminetetraacetate, a- 
tocopherol, vitamin Ky, menadione, ribopolynucleotides, serum albumin, coenzyme 
A, or thyroxine; nor can antimycin A-sensitive respiration be released by the 
addition of known coenzymes such as diphosphopyridine nucleotide. The action 
of the R factor is not prevented by addition to the test system of Mg**, the lack 
of which can cause loss of respiratory control and decline in P:O ratios in intact 
mitochondria. '° 

The characteristic activity of the R factor in increasing respiration is not occa- 
sioned by the simultaneous presence of ADP and an intense ATP-ase activity in the 
preparations. The exhaustively dialyzed material contains no ADP or ATP 
detectable enzymatically or chromatographically. After precipitation at pH 5.5 
the R-factor preparations have no ATP-ase activity and stimulate the ATP-ase 
of the mitochondrial fragments only moderately. The data in Table 2 demonstrate 


TABLE 2 
NoN-INVOLVEMENT OF ADP IN RESPIRATORY STIMULATION BY R. Factor* 
Hexokinase Respiration Phosphorylation 
and myumoles BOH mumoles 
Glucose Oxidized) ATP 
202 0 
199 27 
R factor, 0.4 ml. 103 0 
R factor, 0.4 ml. 4 422 3 
ADP, 8 X 10°*M 289 36 
ADP, 8 X 10°* M + 634 
ADP, 2.4 X 10°? M 5 704 
ADP, 2.4 X 10°? M { 5: 838 
* The test system contained 0.015 M phosphate (1.03 * 106¢.p.m. of P®?), pH 6.5, 0.01 M DL-BOH, and mito 


of N) in a total volume of 2.0 ml Where added, the R factor contained 1056 » of ug 


Vi: 


chondrial fragments (2624 
the hexokinase concentration was 0.5 mg 270 K.M, units'3) per ml.; the glucose concentration was 0.01 


Incubated 15 minutes at 22 


that very large excesses of highly purified yeast hexokinase and glucose added to 
membrane fragments whose respiration has been, stimulated by an excess of the R 
factor in the absence of added ADP causes no trapping of esterified phosphate. 
However, the amount of hexokinase added is capable of trapping phosphate efficiently 
at levels of ADP as low as8 X 107° /. It is evident that the release of respiration 
by the R factor is not caused by catalytic action of traces of ADP coupled with 
intense ATP-ase activity; presumably the characteristic action of the R factor 
occurs at some stage in the coupling mechanism prior to formation of ATP. 
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The R-factor preparations contain a mixture of factors which affect the coupling 
of respiration with phosphorylation, since uncoupling of phosphorylation can also 
be produced when great excesses of the preparation are added to the standard 
test system, far beyond the amoung required to release respiration maximally. 
Furthermore such great excesses of the R factor also stimulate ATP-ase activity! 
and inhibit the ATP-P,** exchange reaction.'? This fact suggested that the R- 
factor preparations may contain the isoédctane-soluble uncoupling factor present 
in mitochrome and aged mitochondrial preparations which Hiilsmann et al. have 
recently deseribed.'? — [soéctane extracts of the R-factor preparations were 
found to stimulate respiration in the test system described. However, such extracts 
also caused substantially complete uncoupling of phosphorylation when ADP was 
present; these effects were found to be simulated by oleic acid.'4 The isodctane- 
extractable factor alone is therefore not responsible for the production of “loosely 
coupled” phosphorylating respiration. Although participation of the isodctane- 
soluble material in producing ‘“loose’’ coupling of respiration as here described is 
not excluded, other experiments indicate that a protein fraction is necessary. 
Separation and purification of the factors involved is under way. 

Freshly prepared rat liver mitochondria contain very large amounts of R-factor 
activity. For example, the amount of activity extracted by sonication from fresh 
mitochondria can produce complete loss of respiratory control in an amount of 
membrane fragments derived from about three times the equivalent amount of 
mitochondria. However, since freshly prepared mitochondria are ‘tightly’ 
coupled, the R factor must be present in a latent or precursor form. 

These experiments define three extreme poles of behavior of respiratory systems. 
In “tightly coupled” respiration, as in freshly prepared normal mitochondria, 
there is very little respiration in the absence of ADP, and the P:O ratio is maximum 
when ADP is present. In ‘loosely coupled” respiration the P:O ratio may be 
maximal or nearly maximal in the presence of a phosphate acceptor system, but 
the rate of respiration is as high in the absence of ADP as in its presence. The 
third extreme case is represented by fully “uncoupled” respirations, such as produced 
by DNP, in which ADP is not necessary for maximum respiration, nor is there 
any phosphorylation when ADP is present. Intermediate degrees of behavior are, 
of course, possible. It is important to note that both “loose ’ coupling and “un- 
coupling’? can be produced in mitochondrial fragments by adding substances 
derived from mitochondria, namely the R-factor preparations described here and 
the isooctane-extractable factors described by Hiilsmann ef al.,'* respectively. 
These mitochondrial factors may thus be involved in intracellular or intramito- 
chondrial control of respiration and phosphorylation. 

It is of significance that Hoch and Lipmann have found” that mitochondria 
isolated from the livers of thyrotoxie snimals have a normal P:O ratio when tested 
in the presence of a phosphate acceptor system. However, such mitochondria 
differ from those isolated from euthyroid animals in that they show a high and 
nearly maximum rate of respiration when ADP is omitted, a situation which cor- 


responds to the “loose-coupling” produced by addition of the mitrochoudrial R 
factor described in this paper. It is therefore suggested that the “unmasking” or 


release of its mitochondrial R factor may be subject to physiological or endocrine 
control and that it ean in turn determine the response of the respiratory chain to 
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the ambient ADP concentration in the cell. These relationships are being examined 
further. 

The following general reaction sequence, in which X designates a substance 
capable of forming a “high-energy” complex with a specific respiratory carrier 
during electron transport, has been postulated for the coupling mechanism: '!7 

(1) Carrier + X ‘lectron transport Carrier ~ X 

(2) Carrier~ X¥ + Pic ’ Carrier + P~ X 

(3) P~ X + ADP, > ATP + X 
It is generally assumed that uncoupling agents such as DNP act by causing the 
decomposition of Carrier ~ X or an equivalent high-energy intermediate early in 
the coupling sequence.*:'® Such an action could account for respiratory stimula- 
tion by DNP in the absence of ADP, as well as the uncoupling by DNP. On 
the other hand, the R factor stimulates respiration in the absence of ADP, without 
producing uncoupling when ADP is present. It may therefore be postulated that 
the R factor and ADP compete for the same intermediate in the coupling mechanism, 
such as P ~ X. The mechanism of action of the R factor is being investigated 
further by examination of its effect on the partial reactions of the coupling process. 

Summary.—A protein-containing factor (R factor), which is apparently “latent” 
in fresh rat liver mitochondria, is released in active soluble form by sonic oscillation. 
It was found capable of fully releasing respiration of mitochondrial subfragments 
from its dependence on ADP as phosphate acceptor, without lowering appreciably 
the P:O ratio when ADP is also added. The R factor thus converts “tightly 


coupled” into “loosely coupled’? phosphorylating respiration and is therefore 
qualitatively different in action from 2,4-dinitrophenol, which produces “uncoupled”’ 
respiration. The significance of these three poles of respiratory behavior, as well 


as the mode of action of the R factor, is discussed briefly. 

* Supported by grants from the U. 8. Public Health Service, the National Science Foundation, 
the Nutrition Foundation, Inc., and the Whitehall Foundation. 

* U.S. Public Health Service Research Fellow of the National Cancer Institute. On sabbatieal 
leave of absence from the Department of Agricultural Chemistry, Oregon State College, Corvallis. 
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GLYCOGEN SYNTHESIS FROM GLUCOSE, GLUCOSE-G-PHOSPHATE, 
AND URIDINE DIPHOSPHATE GLUCOSE IN MUSCLE PREPARATIONS* 
By Puitirps W. Roppins, Ropert R. Traut, AaNpd Fritz LipMANN 
ROCKEFELLER INSTITUTE, NEW YORK CITY 
Communicated November 7, 1958 
INTRODUCTION 


The function of insulin as a regulator of permeability of cell membranes is, at 
present, in the foreground of attention. Nevertheless. experiments such as those 
reported in the lucid paper by Kipnis and Cori! on the change of the temperature 
coefficient for pentose transfer into diaphragm with insulin addition indicate a 
possible shift from diffusion to a chemical and therefore probably enzymatic reac- 
tion. Such observations suggest that the search for a possible primary metabolic 
effect of insulin should not yet be fully abandoned. These considerations led us to 
investigate pigeon breast muscle, where Krebs had found effects of insulin on res- 
piration maintenance? which had been confirmed.* The fact that such an effeet 
appeared to be unique with pigeon muscle seemed, essentially, not to argue against 
its possibly reflecting an aspect of insulin activity. 

Although with the pigeon breast muscle system no insulin effect could be ob- 
tained when measuring glycogen synthesis, we continued because we became aware, 
during these studies, of some seemingly unresolved problems in the field of glycogen 
synthesis proper. Certain problematic aspects of the relation of phosphorylase to 
glycogen synthesis had been brought out already through work on the epinephrine 
effect on phosphorylase activity. This work, done mainly by Sutherland and his 
co-workers!~® and by Cori,’ showed clearly that phosphorylase activation invari- 
ably led to glycogen breakdown and never to glycogen synthesis. Since our ex- 
periments with breast muscle homogenates gave strong support to the proposition 
that phosphorylase was not involved in glycogen synthesis, we considered the possi- 
bility that the synthesis might go by way of UDPG. The demonstration by 
Leloir®’ that liver contains an enzyme which will form glycogen from UDPG, made 
it seem likely that the problem of muscle glycogen synthesis could be explained 
in terms of a UDPG pathway. In the experiments reported below, results with the 


pigeon breast muscle homogenate are described, and preliminary data are given on 


the rabbit and rat skeletal muscle system which synthesizes glycogen from UDP. 


MATERIALS AND METHODS 

C'Jabeled glucose was purchased from Volk Radiochemical Corporation. — It 
Was converted to glucose-6-phosphate with crystalline yeast hexokinase, the gift 
of Dr. M. Kunitz, and was converted to UDPG by a modification of the method of 
Glaser and Brown.’ Both glucose-6-phosphate and UDPG were purified by 
chromatography on Dowex-1 before being used. The UDPG was free of UDP- 
galactose. ATP, DPN, TPN, and CoA were purchased from the Pabst Labora- 
tories. Pigeon breast muscle kochsaft was prepared secording to Krebs and 
Kggleston.? 

Homogenate System for the Incorporation of Labeled Glucose into Glycogen. 
Pigeons were decapitated or anesthetized with Nembutal, and a section of breast 


6 
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muscle was removed rapidly and chilled. The muscle was minced and then ho- 
mogenized in 4 volumes of oxygenated 0.15 M KCl, 0.05 M glyelglycine, pH 7.6. 
Homogenization was carried out in a glass tube with a Teflon pestle. Incubations 
were in Warburg flasks in O» at 37°, with KOH in the center well. The flasks con- 
tained 0.6 ml. of homogenate, radioactive glucose, and other additions, such as 
kochsaft and Na aspartate, in a final volume of | ml. and were incubated for 60-75 
minutes. In most experiments the glucose concentration was 2 X 107% M, 100,000 
cpm/ml. ‘The reaction was stopped by removing the flask and mixing the contents 
with | ml. of 60 per cent KOH. The mixture was heated in a boiling-water bath 
for 30 minutes, and 0.25 ml. of saturated Na,SO, and 3 ml. of 95 per cent ethanol 
were added. After standing on ice, the precipitated glycogen was centrifuged, re- 
dissolved in 2 ml. of water, and reprecipitated with 3 ml. of ethanol. The final 
precipitate was dissolved in | ml. of water, plated, and counted. The glycogen 
isolation method is very similar to the method described by Hassid and Abraham." 
The loss due to the solubility of glycogen in the alcoholic solution was calculated 
to be about 50-100 ueg.!! 

Glucose Incorporation into Glycogen in Pigeon Breast Muscle Homogenates..—W hen 
experiments were carried out under conditions similar to those described by Krebs? 
for the insulin effect, it was found that glucose was not incorporated into glycogen. 
No incorporation was observed until the system was modified by the elimination 
of phosphate buffer. As shown in Table 1, the incorporation depended almost 
completely on the addition of either a muscle kochsaft or CoA. The kochsaft- 
stimulated system was found to be inhibited 90 per cent by 0.01 17 phosphate and 
100 per cent by 5 & 107° M 2,4-dinitrophenol. 


TABLE 1* 


STIMULATION OF GLUCOSE INCORPORATION INTO GLYCO- 
GEN BY KocusAFT AND CoA 

Addition to Incorporation 

Homogenate ¢.p.m 
None 200 
Kochsaft 12,000 
ATP 800 
DPN and TPN 900 
CoA 7.700 
CoA, ATP, DPN, and TPN 6,500 

* Incorporation experiments were carried out in Warburg 

flasks for 60 minutes in Oc, as described in the text rhe amount 
of kochsaft added was 0.2 ml./ml of incubation ATP concen 
tration was 5 X 1074 M, while 50 ug/ml each of DPN and TPN 
were added. The CoA concentration was 5 units/ml 


Since the energy supply for glucose incorporation in this system was derived 
mainly from oxidative phosphorylation, the effects of kochsaft and CoA became 
clear when it was found that at least one of these materials was required in order 
to maintain oxygen uptake throughout the incubation period. In later experi- 
ments, aspartic acid was also added to the incorporation system. It was found 
that aspartate increased the incorporation, probably by helping to maintain the 
adenine nucleotide concentration at a level adequate for normal oxidative phos- 


phorylation.’ Stimulation of incorporation by insulin was never observed. 


Epinephrine inhibited incorporation only at very high concentrations (5-50 ug/ml). 
Effect of Phosphorylase Kinase on Glucose Incorporation System.—We now pro- 
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ceeded to explore the relationship of phosphorylase to this incorporation of glucose 


into glycogen observed in the homogenate system. Although a direct answer to 
this question was not likely to be obtained in the complex homogenate system, it 
seemed possible to get at least a partial answer by studying the effect of phos- 
phorylase kinase. This enzyme has recently become available in a highly active 
state through the work of Krebs and Fischer. '* 

Muscle homogenates were supplemented with crystalline phosphorylase 6'* and 


— Control 
--- Phosphorylase 
kinase added 


’-Percent phosphorylase a 
4-Glucose incorporation ~~moles/ml. 


Time in minutes 

Fic. 1.—Phosphorylase kinase effect on 
glucose incorporation into glycogen. The in- 
cubations were carried out in Warburg flasks 
as described in the text with 0.2 ml. of koch- 
saft, 0.008 M Na aspartate, and 1 mg. of 
crystalline phosphorylase 6 (Fischer and 
Krebs, J. Biol. Chem., 231, 65, 1958) added 
to all flasks. Phosphorylase kinase was pre- 
pared according to Krebs and Fischer ( Bio- 
chim, et Biophys. Acta, 20, 150, 1956) and was 
added to a final concentration of 0.2 mg/ml. 
Since the enzyme is stored in a sucrose-phos- 
phate mixture, a comparable amount of su- 
crose-phosphate was added to controls. The 
reactions were stopped by adding 0.1 ml. of 
1 M NaF 0.1 M EDTA from the side arm. 
After small samples had been diluted for 
phosphorylase assay (G. T. Cori, B. Illing- 
worth, and P. G. Keller in 8. P. Colowick and 
N. O. Kaplan, Methods in Enzymology |New 
York: Academic Press, Ine., 1955], 1, 200), 
the rest of the sample was used for the assay 
of incorporation, 


any sizable incorporation. 


were incubated in pairs, one with and 
one without phosphorylase kinase added. 
As ean be seen in Figure 1, the addition 
of kinase leads to an appreciably higher 
steady-state phosphorylase a concen- 
tration than in the This 
higher phosphorylase a concentration, 


control. 


however, not only does not stimulate 
but, in fact, completely eliminates the 
incorporation of glucose into glycogen. 
A similar experiment is shown in Fig- 
Here labeled 
phate or labeled  glucose-6-phosphate 
and kinase were added after 8 minutes 
of incubation. On the left, it 
seen that the kinase lowers incorpora- 
tion. On the right, a corresponding 
increase in the rate of glycogen break- 


ure 2. glucose-6-phos- 


can be 


down may be noted. 
Fractionation of the homogenate 
system was carried out by differential 
centrifugation. The results of such a 
fractionation are Table 2. 
Since it had not been possible to sub- 


shown in 


stitute energy sources such as ATP, 
creatine P, or PEP for oxidative phos- 
phorylation in the homogenate system, 
it was necessary, of course, to include 
mitochondria in the incorporating sys- 
A surprising result, however, was 
the particle fraction 


tem. 
the finding that 
precipitating between 20,000 * g and 
100.000 X g seemed to be required for 


An explanation of this effect appeared subsequently 


when, on determining the intracellular distribution of the UDPG glycogen synthe- 
sizing system, activity was found chiefly in the particle fraction (Table 4). Al- 
though indirect, this experiment provides preliminary evidence that incorporation 
in the homogenate system was probably by way of the UDPG system. 


Synthesis of Glycogen from UDPG. 


When C'4-labeled UDPG was tested with 


muscle homogenates, it was found to be incorporated much more rapidly than 
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TABLE 2* 


DIFFERENTIAL CENTRIFUGATION OF PIGEON BREAST MuscLe System Tua? 
INCORPORATES GLUCOSE INTO GLYCOGEN 

Glucose 
Incorporation 

(e.p.m.) 
Mitochondria 100 
Mitochondria + 20,000 & g supernatant +, 800 
Mitochondria + 100,000 X g supernatant 600 
Mitochondria + 100,000 X g supernatant + glycogen 500 
20,000 X g. supernatant. < 100 
100,000 X g. supernatant < 100 


Fraction 


* Pigeon breast muscle homogenate was prepared as described in the text but was centri 
fuged for 30 minutes at 20,000 X g. Part of the supernatant fraction was centrifuged 
for 60 minutes at 100,000 X g and the precipitate discarded. Mitochondria were isolated 
by differential centrifugation (J. W. Harman and U. H. Osborne, J. Exptl. Med., 98, 81, 
1953) after resuspension and dilution of the 20,000 xX g precipitate. Incubations were 
carried out with 0.1 ml. of mitochondrial suspension and 0.6 ml. of the appropriate super- 
natant fraction in a final volume of 1 ml., and 0.2 ml. of kochsaft was added to stimulate 
incorporation. Glycogen was added to a final concentration of 0.2 per cent where indi- 
cated. After 1 hour at 37° in O2, glycogen was isolated and counted as described in the 
text. 


glucose or glucose-6-phosphate. In faet, as shown in Table 3, when incubations 
were carried out under anerobic conditions with a small amount of dilute ho- 
mogenate, incorporation was obtained only with UDPG. The maximum incorpora- 
tion rate was found at pH 7.5, with 0.01 17 Mg++. Addition of primer glycogen 
caused an increase in the rate of incorporation. The reaction is approximately 
Jinear with enzyme concentration and with time. Paper electrophoresis confirms 
that the only product formed from UDPG is glycogen. 

Homogenates or extracts prepared in the absence of EDTA had very little 


activity. The differential centrifugation of a rat muscle homogenate is shown in 
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Fic. 2.—Phosphorylase kinase effect on glucose-6-phosphate incorporation into glycogen. The 
experiment is similar to that described in Fig. 1, except that two large-scale incubations were 
carried out in 125-ml. Erlenmeyer flasks. The glucose-6-phosphate concentration was 5 X 10-5 M, 
100,000 epm/ml. Samples were withdrawn at the times indicated for phosphorylase assays and 
incorporation determinations. For the chemical glycogen determinations, samples of the same 
isolated glycogen which was counted were hydrolyzed and assayed for reducing sugar, as de- 
scribed by Hassid and Abraham (in S. P. Colowick and N. O. Kaplan, Methods in Enzymology 
|New York: Academie Press, Inc., 1957], 3, 34). 





BIOCHEMISTRY: ROBBINS ET AL. Proc. N. A. 


TABLE 3* 
IN« ORPORATION OF GLue OSE, Gut COSE-6-PHOSPH ATE, AND 
UDPG invo GiLycoGEN BY Rat SKELETAL MUSCLE 
HoOMOGENATE 


Incorporation 
Substrates ¢.p.m,) 


C'*-glucose, ATP, UTP < 50 
C'*-glucose-6-phosphate, UTP < 50 
C'UDPG ; 1,400 


* The incubation medium contained 0.05 M glycylglycine, pH 7.5 
0.01 M MgCh, 1 per cent rabbit liver glycogen, and enzyme in a final 
olume of 0.5 ml. The reaction was started by adding the labeled 
substrate after a brief pre-incubation at 30°. Incubations were 
carried out for 15 minutes at 30°. The reaction was stopped by add- 
ing 0.5 ml. of 60 per cent KOH, and the tubes were mixed, covered, 
and heated in a boiling-water bath for 10 minutes. The rest of the 
procedure was carried out as described for the pigeon breast muscle 
assay system except that the volumes were half as large. 

In this experiment, 0.02 ml. of a 1:20 homogenate of rat skeletal 
musele in 0.25 M sucrose 107-2 M EDTA was used as the source of 
enzyme. The labeled compounds were all present at a concentra 
tion of 4 & 1074 M and 10,000 ¢.p.m. Nucleotides were added to a 
final concentration of 10 MV 


Table 4. All the activity remains in the supernatant fraction after 10 minutes at 
20,000  g, while the 100,000 & g supernatant fraction is inactive. The particles 
alone have some activity, but, as shown in 
Figure 3, the addition of the dialyzed super- 
natant fraction and = inorganic phosphate 





causes substantial stimulation. Essentially, 
the same particle-supernatant system can be 
obtained from muscle minces by extracting 
at room temperature for 60 minutes with 
0.01 M EDTA. Mince extracts made in the 
cold or in the absence of EDTA have little, 
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if any, activity. 
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Many of the chemical and biological aspects 
01 — of glycogen synthesis from UDPG remain ob- 
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scure at the present. It seems likely, how- 
ever, that transfer of glucose to a glycogen 
primer leads to the formation of a 1,4-linked 
Fic. 3. Stimulation of UDPG in- polymer which undergoes branching under the 
corporation into glycogen by the di- i 
alyzed supernatant fraction from mus- 
cle. The particles sete for this ex- dase. Preliminary observations have shown 
periment were prepared from room- = ar . 
temperature extracts of rabbit muscle that when UDPG is incubated with the part- 
by acid precipitation, The Spinco jeles alone, a polysaccharide is formed which 
supernatant fraction was dialyzed over- : 
night against 0.01 17 EDTA and gave 
no stimulation in the absence of added — that only partial branching has taken place 
phosphate. It gave no incorporation , P a a 
bv itself in the presence or absence of with the particles and l DI G alone. 
phosphate. Incorporation experiments The reason that the particle reaction is 
were carried out as described in Table 3, ; et : 
except that 4X 10-4 M K.HPO, was stimulated by the addition of supernatant 
added to all tubes. The dialyzed protein is not yet clear. The undialyzed 


supernatant fraction was diluted with : ; ; ‘ 
cold water before being used. supernatant fraction stimulates, but after 


influence of amylo-(1.4 —  1,6)-transglucosi- 


gives a purple color with iodine, indicating 
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dialysis no stimulation was found unless inorganic phosphate was added. This 
observation suggested the need to check phosphorylase involvement in the reac- 
tion. It was found, however, that crystalline phosphorylase had little effeet on 
the incorporation of UDPG by the particles. To exclude a nonspecific effect, 
serum albumin was added and found inactive. Experiments are being continued 
at present to define the activity of the supernatant fraction. 


TABLE 4* 


DIFFERENTIAL CENTRIFUGATION OF UDPG-INCORPORATING 
SYSTEM OF SKELETAL MUSCLE 


Incorporation 
Fraction ¢.p.m.) 


Whole homogenate +, 500 
20,000 X g supernatant fluid +, 050 
100,000 X g supernatant fluid 150 
100,000 X g precipitate 1,200 

* Rat skeletal muscle was homogenized in 10 volumes 0.25 M 
sucrose 1072 M EDTA and was centrifuged successively for 10 minutes 
it 20.000 * g and 60 minutes at 100,000 * g. The 100,000 x g 
precipitate was resuspended in the original volume of sucrose EDTA 
Incubations were carried out with 0.025 ml. of the appropriate frac 
tion in a final volume of 0.5 ml Assays were carried out as de 
scribed in Table 3. 


Preliminary experiments have shown that the activity in the 20,000-100,000 x 
g particle fraction is firmly bound. Since this fraction contains a high concentra- 
tion of glycogen, the possibility that the enzyme is strongly adsorbed to glycogen 
must still be considered. The fact that almost half of the phosphorylase 6 of a 
muscle homogenate is also found in this fraction strengthens this possibility. 

The thermodynamic advantages of glycogen synthesis by way of UDPG are 
obvious. Calculations based on data for the sucrose-synthesizing enzyme, sucrose 
phosphorylase, and phosphorylase!® |? show that the transfer of a glucose unit from 
UDPG to glycogen should have a standard free-energy change of about —3,000 ea 
Conversion to glycogen, therefore, should be greater than 99 per cent. In con- 


firmation, transfers of this order of magnitude have already been observed. 


SUMMARY 


In a pigeon breast muscle homogenate, the synthesis of glycogen was studied by 


following the incorporation of radioactive glucose and glucose-6-phosphate. —In- 
corporation depended on active respiration, which was not replaced by ATP or 
energy-rich phosphate feeder systems. Maintenance of respiration in this prep- 
ation depended on the addition of muscle kochsaft, which could be almost com- 
pletely replaced by the addition of CoA and aspartie acid. 

Addition of phosphorylase kinase increased the concentration of phosphorylase a 
in the system and strongly inhibited the incorporation of glucose and glucose-6- 
phosphate. Insulin had no effeet; adrenalin inhibited in high concentrations. 

Glycogen synthesis from UDPG was found to be highly concentrated in particles 
which had microsomal sedimentation characteristics. These had to be supple- 
mented for optimal activity by a nondialyzable supernatant factor, Mg**. and 
phosphate. 

The probable function of the UDPG enzyme as the cellular system for glycogen 
synthesis is discussed. 
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* This investigation was supported in part by a research grant from the National Cancer 
Institute, National Institutes of Health, U.S. Public Health Service, and the National Science 
Foundation. 

The following abbreviations have been used: UDPG, uridine diphosphate glucose; ATP, 
adenosine triphosphate; DPN, diphosphopyridine nucleotide; TPN, triphosphopyridine 
nucleotide; CoA, coenzyme A; EDTA, ethylene diamine tetraacetic acid; UTP, uridine tri- 
phosphate; creatine P, creatine phosphate; PEP, phosphoenol pyruvate. 
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DISSOCIATION OF HUMAN CARBON MONOX YHEMOGLOBIN 
AT HIGH pH* 


By ULricn HAssERopDT AND JEROME VINOGRAD 


GATES AND CRELLIN LABORATORIES OF CHEMISTRY,T CALIFORNIA INSTITUTE OF TE¢ HNOLOGY, 


PASADENA, CALIFORNIA 
Communicated by Linus Pauling, February 20, 1958 


This communication presents a study of some of the hydrodynamic properties 
of the half-molecules of human carbonmonoxyhemoglobin at pH 11.0 and ionic 
strength 0.25, together with a demonstration that, at this pH, human carbonmon- 
oxynemoglobin is completely dissociated and the molecular weight of the product is 
one-half that of the whole molecule. The half-molecules appear to be stable when 
stored at pH 11 and 3° C., in contrast to their behavior in an acid environment, 
pH 4.5-6, in which denaturation occurs. The acid denaturation is so rapid at pH 
values low enough to induce complete dissociation that it has not been possible 
until now to study the properties of the half-molecules uncomplicated by the pres- 
ence of the whole molecules. Results of a comparison of some properties of un- 
treated hemoglobin and hemoglobin reconstituted from half-molecules by successive 
dialyses at pH 11 and pH 7 are also presented. 

Experimental.—C1 ystallized carbonmonoxyhemogiobin was prepared according 
to the method of Drabkin.! This material was used in most of the sedimentation 
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studies. Hemoglobin obtained prior to crystallization was also employed, and no 
significant difference was found in the sedimentation behavior at either pH 7.0 or 
pH 11.0. 

The following buffer solutions were used: Tris-HCl-NaCl, pH 8-9; glycine 
NaOH-NaCl, pH 9.6-11.6; and phosphate-NaCl, pH 11.0-11.6. The ionic 
strength was 0.25 in all cases. All hemoglobin solutions were dialyzed at 3° C. for 
12-15 hours before the sedimentation runs. 

The sedimentation and diffusion measurements were carried out in the Model E 
Spinco analytical ultracentrifuge, using a phase-plate schlieren diaphragm. A 
plastic commercial interference synthetic boundary cell? was used for the diffusion 
measurements, which were performed with an essentially static boundary at 8,000 
rpm. Sedimentation velocity measurements in the pH range 10.0-11.6 were carried 
out in a double-cell rotor, using a 1° positive wedge. The second cell contained 





| ! T | I 








| l | | | i | l 
4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 pH 


Fic. 1.—Sedimentation coefficients of human carbonmonoxyhemoglobin in the pH range 4.5 
11.6. @ from Field and O’Brien (Biochem. J., 60, 656, 1955): Hb concentration: 0.71 per cent. 
Ionic strength of the solution: 0.1. © this study. 








carbonmonoxyhemoglobin at pH 11.0 as a reference. Thus differences of 0.1 sedi- 
mentation units could be detected with considerable certainty. 

The sedimentation and diffusion coefficients are expressed in the text in ¢.g.s. 
units X 10'* and 10’, respectively. 

Results and Discussion.—Table 1 shows the average sedimentation coefficients 
obtained in basic solutions. The trend of the values indicates that the dissociation 
of the hemoglobin molecule in 


a | per cent solution begins just beyond pH 10. 
At pH 11.0 the value of sxo,, is 2.55. This is the value reported for horse globin® 


o. 
and also for human globin‘ at pH 7.0 and at concentrations of 1 per cent. Both 


these materials have half the molecular weight of the corresponding hemoglobins.*4 
Figure 1 shows a small plateau in the s versus pH curve around pH 11.0, indicating 
that in this region hemoglobin is completely dissociated into two subunits. At 
this pH the dependence on concentration is normal, as would be expected for a fully 
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dissociated hemoglobin, while at pH 10.5 the results suggest the presence of incom- 
plete dissociation (Table 1). 

Molecular weight determinations at pH 11.0, using the Archibald method of ap- 
proached sedimentation equilibrium,® gave an average value of 36,000 + 2,000.° 
Involved in this result is the assumption that the partial specific volume does not 
change significantly from 0.749 ml. g.7}. 

The spreading of the boundaries at pH 10.5 and 11.0 is in accord with the above 
conclusions. The sedimenting boundaries after 1 hour or longer at 56,100 rpm 
were skewed at pH 10.5, and the spreading was greater than could be accounted for 
by diffusion alone. At pH 11.0 the boundaries had an approximately Gaussian 
shape, with less spreading than at pH 10.5 or 11.6. Thus, in plots of the standard 
deviation of the schlieren patterns o versus +/t, curves with an upward trend were 
obtained for Hb at pH 10.5 and 11.6, and a straight line at pH 11.0. 

A more extended study at pH 11.0 gave the following results: the average diffu- 
sion coefficient (static) in phosphate as well as glycine-NaOH buffers is 7.6 + 0.4; 
the concentration dependence lies within the limits of error; s°29.. has the value 
2.7, and D°o9.., calculated from s and VW, has the value 7.4. 

Moving-boundary electrophoresis at pH 11.0, 7 = 0.1, phosphate buffer, showed 
a single symmetrical boundary in the ascending and descending limb. A solution 
first dialyzed to pH 11 and then to pH 6.75 gave electrophoretic patterns indistin- 
guishable from an equivalent sample dialyzed to pH 6.75 directly. A mixture of 
equal quantities of the foregoing materials again gave the same pattern. 

The frictional ratio calculated for the subunit at pH 11.0 is 1.27, whereas this 
value for hemoglobin at pH 7.1 is 1.21, if s°x9,. 1s taken to be 4.45 and the molecular 
weight 66,500.27 This slight increase in the molecular asymmetry, which corre- 


sponds, for example, to a change in the axial ratio from 3:1 to 4:1 in an ellipsoid, 


assuming 20 per cent hydration,’ suggests that cleavage is along the major axis of 
the molecule. 

The dissociation into subunits could be reversed upon neutralization. Samples 
of hemoglobin which had been at pH 11.0 for 12 and 36 hours, respectively, were 
dialyzed to pH 8.2, and the sedimentation coefficient was measured along with 
reference samples at the same pH. ‘The result was 4.25. 

Neither precipitation nor change in the spectrum (350-650 my) was observed on 
raising the pH to 11.0 and maintaining this pH for 36 hours. The region of the 
hemoglobin spectrum most sensitive to denaturation and shifting of the hemes is to 
be found at the maxima at 540 and 568 my. The ratios of the minimum at 376 
my to the two maxima were 1.26 and 1.29 at pH 7.1 and 1.25 and 1.28 for Hb which 
had been at pH 11.0 for 36 hours and was redialyzed to pH 7.1. 

On the other hand, the changes in the hemoglobin molecule at pH 11.6 are not 
reversible. The color of the solution changes, and the maxima at 540 and 568 my 
disappear slowly. Precipitation occurs on redialysis to pH 7.1, and the amount of 
this precipitate increases with the time of exposure to the alkali. 

From the data in Table 1 and the results of Field and O’Brien’ a pH. stability 
diagram may be formulated (Fig. 1). This curve is not symmetrical about the 
isoelectric point, indicating that dissociation is not entirely the result of repulsions 
generated by the net negative or positive charge on the protein. Neglecting ion 
binding, it may be estimated” that at pH 6, where acid dissociation begins, the 
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TABLE 1 
SEDIMENTATION AND DIFFUSION COEFFICIENTS OF HUMAN CARBONMONOXYHEMOGLOBIN * 


1 Per Cent 0.3 Per Cent 


' pH ie Dies 


9 


* All measurements at yp 0.25 in buffered solutions described in text 


t Runs in double-cell rotor with a solution at pH 11 present as a reference 


molecule has a net charge of +5, whereas at pH 10 the charge is about —30. These 
considerations neglect binding of negative ions, which may be expected to make 
the situation even more asymmetric. Thus it is not likely that general net charge 
repulsions are solely responsible for the dissociation reaction at low pH. 

An examination of the stability curve suggests that the amino acid residues 
critically involved in stabilizing the whole molecules are histidines, which are pro- 
tonated during acid dissociation, and lysines and tyrosines, which are neutralized 
or ionized upon alkaline dissociation. 

The temperature coefficients of the dissociation into subunits have been esti- 
mated from the temperature coefficients of the observed sedimentation coefficients 
at pH 5 and 10.5. The results corresponded to 0.0 + 0.5 keal./mole at pH 5 and 
—4.0 + 0.5 keal./mole at 10.5. These enthalpy changes are interpreted as the 
sums of the enthalpies of splitting the molecules plus the enthalpy changes involved 
in charging or discharging the critical groups. Since the number of the groups is 
not known, it is not possible at the present time to apportion these enthalpy changes. 


* This work was supported in part by a grant No. H-8394 from the United States Public Health 
Service. 
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RADIATION INDUCED REVERSE MUTATIONS IN: DROSOPHILA 
MELANOGASTER 


By M. M. GREEN 
DEPARTMENT OF GENETICS, UNIVERSITY OF CALIFORNIA, DAVIS 


Communicated by G. Ledyard Stebbins, November 7, 1958 


There can be little doubt that within the limits of current operational criteria 
spontaneous reverse (back) mutations have been recorded in Drosophila and other 


organisms.' Doubt has been raised whether the back-mutation frequency can be 
significantly increased following X-ray irradiation, despite early reports that it can.? 

The relatively high spontaneous back-mutation frequencies of two pseudoallelic, 
recessive, sex-linked forked bristle (f) mutants in Drosophila melanogaster suggested 
that they might provide critical information bearing on the question of induced re- 
verse mutation following X-ray or gamma treatments. These mutants, f*” and 
f', were found to back-mutate spontaneously in females at rates of 46.7/10® (14/ 
298,698) and 15.3/10® (7/457,316) chromosomes, respectively. ! 

By design, attached-X 9 9 of genotypes v f*"car and v f'B were treated with 
1,000 r X-rays or gamma, and reversals were sought among their female progeny. 
The use of attached-X 2 9 permits the scoring of two treated chromosomes for 
each 2 zygote. 

Adult females 0-48 hours old were irradiated with either 4000 r X-rays delivered 
from a General Electric ‘““Maxitron”’ machine operated at 30 m.a. and 100 KVP 
(rate ca. 1,000 r/min.) or 4,000 r gamma delivered from a Co® source (rate 265,000 
r/hour).* Treated 2 9 following mating to f! oo were allowed to oviposit for 
7 days and transferred to fresh media for a second 7-day interval. Partial and 
complete f reversals were sought among F; 9 9. All 2 9 showing phenotypic 
reversals were progeny-tested by crossing to o sibs. Progeny tests were obtained 
onall 2 2 discussed herein. 

The results of these experiments, listed in Table 1, represent summations of three 
separate X-ray treatments of v f*"car 9 Q and two of v f'B 2 2, together with 
two separate gamma treatments of v f*"car 9 9 and onev f'B 9 9. The reversals 
recorded were all complete and phenotypically inseparable from wild-type. All 
bred true on progeny-testing. Although a number of partial reversals were re- 
corded, none bred true. 

TABLE 1 


CompLeETE REVERSALS OF f3" AND f! IN ArracHED-X 92 2 FoLtLowing X-Ray or GAMMA 
IRRADIATION 


Radiation Genotype Total f Total Reversals / 10! 
4,000 r) Treated 9? § ' Chromosomes Chromosomes 


ace vf" car 38 115,774 328.2 

} vf'B | 68 , 268 14.6 

Gamma fo f®car 10 65, 206 153.3 
aus vf'B 0 21,550 0 


A comparison of the induced and spontaneous rates of reversals shows that 
1,000 r X-rays increased f*" reversals seven fold, whereas equivalent gamma irradi- 
ation somewhat more than tripled the reversal frequency. There is little doubt 
that radiation at the dose given significantly increased the reversal frequency. 


16 
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The experiments with f! show that no comparable significant increase in reversal 
frequency followed irradiation. Admittedly, the less extensive nature of these 
experiments permit only a tentative conclusion. 

Incidental to scoring f reversions, flies were also scored for possible reversions of 
v, car, and B. No reversals of v or car were found. At least one partial reversal 
of B was found, and this is not unexpected, since induced partial and full reversions 
of B have been recorded. 

By progeny-testing reversal 2 2, two relevant questions can be answered: (1) 
whether one or both treated f mutants reverted following irradiation and (2) whether 
gross chromosomal rearrangements are causally associated with reversion. The 
occurrence of crossing over in attached-X 9 9 provides answers to both questions. 
If but one f mutant reverted following irradiation, then the unaltered f should be 
recovered in homozygous condition as a consequence of nonreciprocal crossing over 
between the f locus and the centromere. If no gross chromosomal rearrangement 
occurred, the frequency of homozygosis of the unaltered f mutant should equal one- 
fourth the standard frequency of crossing over between f and the centromere. 

Without presenting the extensive data, the following conclusions can be made. 
First, homozygous f*" or f! 2 2 were recovered among the progeny of all reversals, 
proving that only one of the two treated f alleles reverted. Second, in each case 
the frequency of homozygosis was consistent with the interpretation that crossing 
over between the f locus and the centromere was unencumbered by any gross chro- 
mosomal rearrangement. 

Ample evidence has accumulated to warrant the conclusion that Drosophila cells 
in different stages of spermatogenesis differ significantly in their sensitivity to X- 
irradiation.4 A similar situation during oégenesis is indicated by studies of induced 
dominant and recessive lethals.° Following X-irradiation, the frequency of domi- 
nant and recessive lethal drops sharply after the seventh day of oviposition, pre- 
sumably because the supply of mature odcytes exposed at the time of treatment 
have been replaced by cells that were younger (o6gonia) at that time. Since 
dominant lethals are, in all probability, associated with chromosomal rearrange- 
ments, it follows that the maturer cells are more prone to induction of rearrange- 
ments than are younger cells. Hence a comparison between reversion rates during 
the two oviposition periods could demonstrate whether or not reversions are usso- 
‘ated with chromosomal rearrangements. 

A consideration of the f*" reversions induced by X-rays shows that the rates for 
the first and second weekly egg-laying intervals are essentially identical. For the 


first week, 8 reversions/29,452 chromosomes were recorded; for the second week, 


30/86,322. The rates computed per 10° chromosomes are 271.6 versus 347.5 and 
appear not to be significantly different. This comparison supports the conclusion 
stated earlier that the reversions appear to be independent of chromosome rearrange- 


ment. If such an association exists, a preponderance of reversions would be ex- 
pected among the progeny of the first 7-day egg-laying interval. The data sug- 
gest that perhaps most, if not all, of the reversions are induced in young oécytes 
or oégonia. Among D. melanogaster females 0-48 hours old, the majority of ovarian 
cells at the time of treatment will be at the oégonial or early od6cyte stage of develop- 
ment. Exposure to 4,000 r X-rays or gamma should induce dominant lethals in 


most mature odcytes, thereby leaving only the viable oégonia and young odcytes 
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eligible for reversions. If the latter cells are the primary source of reversions then, 
as observed, no difference is expected in rate between the consecutive egg-laying 
intervals. 

The question of whether the phenotypic reversion of recessive mutants to wild 
type or near wild type is in fact back mutation has been discussed in extenso else- 
where in connection with the problem of spontaneous back mutation.! The argu- 
ments militating for spontaneous back mutation as a specific change or alteration 
of a particular mutant gene apply equally well for the radiation-induced reversions 
of f*" reported here and need no further comment. In short, these results amplify 
results based primarily on partial reversions in Drosophila® and of reversions in 
Neurospora’ that X-rays and gamma rays can induce back mutations which, inso- 
far as current criteria permit, do not differ from spontaneous back mutations. 

Summary. Drosophila melanogaster 2° 2° homozygous for one or the other forked 
pseudoalleles, f°" and f', both known to back-mutate spontaneously, were irradi- 
ated with 4.000 r X-rays or gamma rays. 

Back mutations of f°" were significantly increased; those of f! apparently not. 

The back mutations are not associated with detectable chromosome alterations. 

M. M. Green, Genetica, 29, 1-88, 1957 
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National Laboratory. Thanks are due Drs. H. J. Curtis and R. C. King for making facilities 


available. 
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6H. J. Muller, and I. I. Oster, Adv. Radiobiology, pp. 407-414, 1957. 

7 N. H. Giles, F. J. deSerres, and C. W. H. Partridge, Ann. N.Y. Acad. Se7., 59, 536-552, 1955; 
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DEPENDENCE OF MUTATION FREQUENCY ON 
RADIATION DOSE RATE IN FEMALE MICE 
By W. L. Russeui, LIANE Braucu Russeti, AND Mary B. Cupp 
BIOLOGY DIVISION, OAK RIDGE NATIONAL LABORATORY, * OAK RIDGE, TENNESSEE 
Communicated by Alexander Hollaender, November 17, 1958 


In studies of specifie locus mutation rates in mouse spermatogonia, a much lower 
mutation rate was obtained from chronic gamma than from acute X-irradiation. !~! 
Differences in both intensity (dose rate) and quality of radiation were involved in 
the comparison, but it was pointed out? that, although gamma rays frequently 
show a biological effectiveness somewhat lower than that of X-rays, the observed 
difference between the mutation rates appeared to be greater than could be attrib- 
uted to this factor. This point of view has since been confirmed by three separate 
experiments conducted to test it.4 In these experiments, no appreciable differences 
were found in the effectiveness of acute gamma rays (from Co”), on the one hand, 
and acute X-rays, on the other, in inducing dominant lethals in spermatozoa, specific 
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locus mutations in spermatozoa and other postspermatogonial stages, or specific 
locus mutations in spermatogonia. It can be concluded that most of the difference 
between the chronic gamma and acute X-ray mutation rate results must be due to 
dose rate. 

However, the finding of a dose-rate effect in spermatogonia did not necessarily 
prove that this effect was on the mutation process itself. It might, instead, have 
been due to secondary causes, such as cell selection resulting from spermatogonial 
killing or from other interference with the dynamies of the eyele of the seminiferous 
epithelium." 2 To determine whether the mutation process itself could be affeeted 
by dose rate, it became important to carry out a dose-rate comparison in females 
where there seems little possibility that mutation-rate measurements can be 
interfered with by such secondary processes as elimination of sensitive cells. The 
adult ovary contains no oégonia, and the population of germ cells consists solely of 
odcytes, Which are non-mitotic and which appear to be remarkably uniform.® 
Furthermore, the continued fertility of females in chronie radiation experiments® 7 
provides no evidence of extensive killing, selective or otherwise, of the oGeytes. 

The first part of the dose-rate comparison for mutation rates induced in oéeytes 
has already been reported. The specific locus mutation rate in adult females 
exposed to continuous chronic Cs!*7 gamma. irradiation at a dose rate of 86 r 
per week proved to be considerably lower than that produced by ncute X-rays in 
spermatogonia.’ Carter’ also observed a low mutation rate in females exposed to 


) 


interrupted chronic Co® gamma irradiation of somewhat higher intensity than that 
used by us. The mutation rates from chronic irradiation of females, as far as they 
go, are not significantly different from the rates obtained from chronie gamma 
irradiation of spermatogonia. 4 

The completion of the dose-rate comparison for mutation rates induced in 
odcytes, namely, the determination of mutation rates in females exposed to acute 
irradiation, promised a clear-cut answer of broad significance. For, if this rate 
proved to. be no higher than that from chronic irradiation, then it would show 
that odcytes are much less sensitive than spermatogonia. On the other hand, if 
the rate from acute irradiation of females turned out to be similar to that from 
acute irradiation of spermatogonia, namely, much higher than from chronic ir- 
radiation of both sexes, then it would indicate that there was a dose-rate effect 
directly on the mutation process. 

It was, therefore, most important to obtain a reliable mutation rate from acutely 
irradiated females, even though this task is beset by the practical difficulty that 


the number of offspring that can be obtained before sterility ensues is very small. 


The experiment reported here Wiis accordingly started to extend the limited data 


available from an earlier pilot experiment.’ \ brief statement of the new findings 
has already been presented. In the present account, the results to date are 
described in more detail. Although the collection of data is still proceeding, the 
results are already sufficient to give a clear answer to the questions posed. 
Materials and Method.Young adult F, hybrid females from a cross of 101 X 
C3H strains were exposed to 400 r of whole-body X-radiation (250 Kwp; 15 ma.; 
inherent filtration 3mm. Al; H.V.L. 0.4 mm. Cu; dose rate 92 r min; animals in 
lucite containers on turntable rotating above masonite seatter block). The 


females were mated for 2 weeks, beginning on the day following irradiation, to 
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males of our T stock, which is homozygous for seven recessive marker genes. 
Young subsequently born were examined for mutations at the seven loci. Details 
of the experimental procedure have been described earlier.® 

The data obtained in the present, 1958, experiment are additive with those of the 
1950-1951 series,® since the stocks used were the same and the methods were 
identical except for one detail: in 1950-1951, when fertility effects of radiation were 
still being explored, females were kept with their mates indefinitely instead of 
being separated after 2 weeks, as in the 1958 series. However, the bulk of the 
young (97 per cent) were conceived within 2 weeks after irradiation, since most 
females were sterile by the time they could have had a second mating, and all were 
sterile 5 weeks after irradiation. 

In 1950-1951, 379 females were irradiated, and 319 of these gave 1,729 young 
surviving to the age at which mutations are scored (i.e. approximately 3 weeks). 
In 1958, 6,130 young have, so far, been scored for mutations. These were off- 
spring of 1,088 irradiated females. The results of the two acute X-ray experiments 
are presented in Table 1, together with data from chronic irradiation of females 
and control data. The data from the chronic gamma experiment have been 


brought up-to-date since the last report.® 


TABLE 1 
Muvuvations At Specrric Locr in Mouské O6cytTEs 


Dost INTENSITY No. OF No. or Murations at Locus 
RADIATION r r/min) OFrrsPRING b ri d p 8 8 d s¢ 


* 0 5,845 
 &. 100 O6 1,729 
X 100 q2 6, 130 


0 10,918 
Gamma (Cs!*7) 258 0.0085 26,468 1 


* 1950-1951 experiment (Russell, Russell, Gower, and Maddus, these Procrrpinas, 44, 901-905, 1958) 


Conclusions.-On combining the 1950-1951 and the 1958 series on acute. ir- 
radiation of oédcytes, the mutation rate is 381.8 X 10~%/r/locus (95 per cent con- 
fidence interval: 12.8-65.5 K 107). The mutation rate calculated from the up- 
to-date tabulation of the chronic gamma data for females is 4.2 K i0-°/r/locus 
(95 per cent confidence interval: 0.51-15.1 & 107%). A test of the significance of 
the difference between the mutation rates per roentgen in the two experiments gave 
P = 0.006 for a one-tailed test. Thus a significantly greater frequency of specific 
locus mutations is induced in oécytes by acute than by chronic irradiation. 

It might be argued that the oédcyte stages irradiated were not strictly comparable 
in the chronic and acute radiation experiments. Thus, in the former, offspring 
were obtained for almost a year, representing, therefore, odcytes irradiated in 
follicles of a variety of stages of development, from youngest to most mature. In 
the acute irradiation experiment, on the other hand, all offspring resulted from 
odcytes irradiated in relatively mature follicles. However, as already noted, the 
state of the odcyte nucleus in the adult ovary appears quite uniform (dictyate) 
throughout the development of the follicle, except for the last few hours before 
ovulation (a stage avoided in the present experiment by not mating females until 
1 day after irradiation). Development of the follicle seems to involve only growth 
in odcyte size and multiplication of follicle cells. Moreover, a difference between 
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chronic and acute irradiation can be demonstrated even if one uses only that part 
of the chronic radiation data which pertains to oécyte stages more comparable to 
those involved in the acute radiation experiment, i.e., only litters conceived within 
2 weeks following accumulation of the dose. In Carter’s experiment,® in which all 
litters were of this classification (because of subsequent sterility from the dose 
distribution used), one mutation occurred among 10,117 offspring. Of the 26,468 
offspring from chronically irradiated females in our experiments, 4,303 were con- 
ceived in the first 2 weeks, and one of the two mutations was among these. These 
two sets of results from chronically irradiated later odcyte stages, taken together, 
give a significantly lower rate than that now found to result from acutely irradiated 
females (P = 0.005, for one-tailed test), indicating that, as suspected, the dif- 
ference between the acute and the total chronic data is not due to odcyte stage. 

The finding of an effect of dose rate on the frequency of mutations recovered 
from irradiated females—-where, as explained, secondary factors, such as cell 
selection, cannot be easily invoked—argues strongly that the mutation process 
itself must be affected by dose rate. Furthermore, although secondary mechanisms 
cannot yet be definitely ruled out as the cause of the intensity effect demonstrated 
earlier for spermatogonia, it is simpler to assume that the explanation for the 
results in o6cytes also applies to spermatogonia. 

That the frequency of radiation-induced mutations is independent of dose 
rate has been one of the basic tenets of radiation genetics. The finding that dose 
rate affects the frequency of specific locus mutations observed following irradiation 
of mouse spermatogonia!~* did not. by itself violate this tenet, since the result 
might have been due to secondary factors. However, the present demonstration 
of a dose-rate effect on mutation rate in odcytes shows that the tenet does not have 
universal validity. The classic findings of intensity independence were made in 
Drosophila spermatozoa.’” Recent experiments on mouse spermatozoa’ are in 


complete agreement with the Drosophila results. Therefore, the explanation for 


the new phenomenon of intensity dependence resides in gametogenic stage. Hy- 
potheses concerning the mechanism for this effect have been discussed elsewhere,’ 
and it was shown that a threshold response for all mutations in spermatogonia and 
odcytes is not a necessary consequence of the findings. 
The specific locus mutation rates in odcytes, 4.2 X 10/r/loeus and 31.8 & 10 

r locus for chronic and acute irradiation, respectively, may be compared with the 
rates in spermatogonia, namely (as calculated from straight lines fitted in the range 
of 0-861 r*), 5.2 K 10-%/r/loeus and 21.3 & 10°%/r/locus for chronic and acute 
irradiation, respectively. It will be noted that, at each dose rate tested, there is, 
at present, no evidence of marked difference between odcytes and spermatogonia. 
Before the mutation rates from acute irradiation of females were known, it was 
argued,® on the basis of the results then available, that the low mutation rate 
obtained with chronic irradiation was due either (a) to the same factor in males 
and females, namely, the low dose rate, or (6) to different factors in the two sexes, 
namely, to secondary causes in males and to a low sensitivity of the sex or of the 
oocyte stage in females. The new results, showing a high mutation rate from 
aeute irradiation of females, obviously disprove the second hypothesis. The 
conclusion,® that “the mutation rate in women exposed to chronic gamma radiation 
may be less than that estimated from our earlier yardstick, namely, the mutation 
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rate obtained from acute X-irradiation of male mice,” can now be extended to 
chronically exposed men also, since the new results have made secondary factors 
unlikely. Carter? thought it most likely that the low mutation rate found 
following chronic irradiation of o6cytes was attributable to sex. As has just been 
pointed out, this interpretation is not upheld by the results reported in this paper. 
Carter’s emphasis on the consequence of his interpretation, namely, that only a 
small part of the genetic hazard from medical irradiation would come from ex- 
posure of females, now appears to have been misleading. 

Of course, more information is needed on the relative mutational sensitivities of 
males and females under various conditions of irradiation—intensity, fractionation 
of dose, ete.—-but there is no justification at present for assuming less genetic 
hazard from irradiation of females. In faet, the present point estimate from the 
mutation rate induced by acute radiation in o6cytes is actually slightly bigher than 
the spermatogonia rate. Although the difference is not significant in the present 
data, the possibility that the mutation rate may be somewhat greater in oécytes 
than in spe:matogonia is not out of line with the higher mutation rates found from 
irradiation of postspermatogonial stages,!' among which are spermatocytes. In 


Drosoniila, specific locus mutation rates are about half to three-fourths as high in 


oocytes as in spermatozoa,!? but the spermatozoa rate is about four times the 


spermatogonia rate.!® Thus the relative mutation rates of o6cytes and spermato- 
gonia may be similar in the two species. 

It may be noted from Table | that no mutations have, to date, occurred among 
16,763 offspring of our control females. Control males have, so far, yielded 17 
mutants among 288,616 young! at our laboratory and 4 mutants among 41,262 
young at Harwell. Thus there is some indication that the spontaneous mutation 
rate in females may be lower than in males in the mouse. <A similar situation 
exists in Drosophila." 

Summary. The specific locus mutation rate in females is significantly higher 
following acute than following chronic irradiation. Since conditions in the adult 
mouse ovary are such that complications from secondary factors are unlikely, it can 
be concluded that in odeytes there is a dose-rate effect on the mutation process 
itself. The influence of dose rate on the frequency of mutations following irradia- 
tion of spermatogonia, which was demonstrated earlier and which could have been 
due to secondary factors, is now, in the light of the odcyte results, also considered to 
be due to a dose-rate effect on the mutation process. 

The most signifieant features of the results, from the point of view of human 
hazards, are, first, that there is, at present, no evidence that females are less sensi- 
tive than males and, second, that the finding of a dependence of mutation frequency 
on radiation dose rate, now extended from males to females, applies to those germ 
cell stages that are important in human hazards. 


The authors are grateful to Miss E. M. Kelly and Mrs. KX. F. Stelzner for their 
extensive assistance in this experiment. 
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MODIFICATION OF TRRADIATION EFFECTS IN THE PIGEON, 
COLUMBA LIVIA, 1. EFFECTS OF HOMOLOGOUS AND 
HETEROLOGOUS BONE MARROW TRANSPLANTS* 


By DonaLp H. SHaw aANnp HALVoR VERMUND 


DEPARTMENT OF GENETICS AND DEPARTMENT OF RADIOLOGY, 
UNIVERSITY OF WISCONSIN 


Communicated by M. R. Irwin, November 
INTRODUCTION 


Previous investigations have demonstrated that the lethal effeets of total-body 
roentgen irradiation on mammals (mice, rats, rabbits, and guinea pigs) can be 
modified by isologous, homologous, or heterologous bone marrow transplanta- 
tion.!°!! These modifications appear to depend upon (1) the establishment and 
proliferation of the hematopoietic elements of the transplant and (2) either partial or 


complete repopulation of the host’s depleted hematopoietic system by the cells pro- 
liferating from the graft.°:47 

The present report includes preliminary findings of experiments in which homol- 
OZOUS and heterologous bone marrow Injections were administered to pigeons pre- 
viously subjeeted to 2,500 rads (approximately the LD joo) total-body X-irradia- 
tion. The evidence indicates strongly that partial and eventually total repopula- 
tion of the erythrocytes results from the homoglous and heterologous transplants. 


Also, host-species-specifie hemagglutinins (hemagglutinins specific for the erythro- 


cytes of host species) appeared in the circulation ot irradiated pigeons which re- 


ceived heterologous bone marrow. 
MATERIALS AND METHODS 
Dosage response studies have shown that the LD5o 39 for the pigeon is approxi- 
mately 2,050 rads and that the effeetive LD oo 39 Was 2,500 rads. These studies also 
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indicated that the death pattern appears to be comparable to that observed in mice, 
rats, ete. This information will appear in detail elsewhere. 

Non-inbred sexually mature pigeons, 6-10 months of age, were individually ex- 
posed to 2,500 rads total-body roentgen irradiation. A single dose was adminis- 
tered by a General Electric Maxitron resonant transformer X-ray machine operated 
at 1,000-KV peak. The target served as a filter, and no additional filtration was 
provided. The half-value layer was 3.4 mm. of lead. The birds were treated in- 
dividually and rotated at 5 revolutions per minute in a lucite cylinder measuring 7.5 
by 20 em., with the longitudinal axis of the animal as the axis of rotation. The 
radiation field covered the entire cylinder with the bird in the middle. The dis- 
tance between the target of the X-ray tube and the axis of rotation was 70 em. 
The dose rate was 92 r/min, as determined with a Victoreen roentgen meter using 
an ionization chamber placed within the cylinder with tissue-equivalent scattering 
material around it to secure electron equilibrium. The chamber had been eali- 
brated by the National Bureau of Standards. The conversion factor from roent- 
gens to rads was 0.96. 

A total of 80 pigeons was used: 60 received 2,500 rads, and 20 served as unir- 
radiated controls. The animals were separated into 4 experimental groups. 

(iroup | consisted of 20 irradiated birds, 10 of which were treated in July, 1958 
(Series A), while the remaining were treated in September (Series B). These birds 
were each given a single intravenous injection of a saline suspension of fresh ho- 
mologous bone marrow within 3 hours after irradiation. The volume was adjusted 
so that each bird received 2 ml. of a suspension containing approximately 250 10! 
nucleated cells. 

The homologous bone marrow was obtained from the femora and tibiotarsi of 10 
backeross hybrids which were derived from mating to dévza of hybrids and selected 
backeross hybrids from crosses between C. guinea and livia. These birds have 
been sufficiently backcrossed to be considered genetically homologous with /zv7za, 
except for the genetic material producing one or both of two erythrocytic antigens, 
B and ©, derived from guinea.® ° Such markers on the donor’s erythrocytes al- 
lowed detection in the host’s circulation of blood cells derived from the bone mar- 
row transplant because neither of these factors, B or C, is present in /zvza. In addi- 


tion, a separation of the erythrocytes of the host from those of the donors was pos- 


sible beeause the erythrocytes of each host gave positive agglutination with saline 
extracts of lima beans Sieva (Phaseolus limensis), while those of each donor were 
not agglutinated by this reagent.!? Each donor supplied bone marrow for two ir- 
radiated hosts, and at no time were marrow samples from two donors pooled. 

Gioup 2 consisted of 20 irradiated birds (10 in Series A and 10 in Series B), each 
of which received 2 ml. of a saline suspension of fresh heterologous bone marrow 
within 3 hours after irradiation. Hach host received approximately 250 X= 105 
nucleated cells, which were obtained from femora and tibiotarsi of a single ring 
dove (Streptopelia risoria). 

The 20 irradiated pigeons in Group 3 (10 each in Series A and B) received a single 
intravenous injection of 2 ml. of the sterile saline. Group 4 consisted of 20 unir- 
radiated animals, 10 each in Series A and B, which were injected as in Group 3. 


Two birds from each of Groups 1-4 were caged together following the irradiation 
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in an effort to reduce the influence 
of differential environmental fac- 
tors to a minimum. 

Blood samples were drawn from 
the wing vein of each animal on 
days 0, 4, 6, 8, 10, 12, 14, 16, 20, 
24, 28, and 30 or until the animal 
died. Hematocrit readings Ol 
these samples determined the ap- 
proximate volumes of red and 
white cells. The samples from 
Group 1-A and B were tested with 
antisera specific for the B and C 
antigenic factors, to reveal the pres- 
ence or absence of the donor’s 
erythrocytes, and with Sieva ex 
tract to detect the host’s red cells. 

The samples from Group 2-A and 
B were tested with rabbit anti-livza 
serum (absorbed with cells of rZso- 
ria, rendering it specific for livia 
erythrocytes in contrast to those 
of risoria) to detect pigeon (host) 
red cells. These samples were also 
tested for the presence or absence 
of donor’s red cells with an anti- 
serum (rabbit anti-rzsoria serum 
absorbed with cells of lévza) which 
Was reactive with risorta erythro- 
cytes, not with those of lévza. 

In Series B, serum samples were 
obtained at each bleeding from the 
individuals in Groups 1+. These 
sera were tested for possible hemag- 
glutinins against astandard panel 
consisting of erythrocyte suspen- 
sions from l2via, résoria, and back 
cross hybrids carrying the guinea 
factors Band C.  Inaddition, each 
serum was tested against cells of 
the bone marrow donor and its 
own erythrocytes obtained on day 
QO. 

All birds were weighed at the 
intervals above for blood sampling 
and those which died were weighed 


as soon as possible after death 
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Autopsies were not conducted, 
since most of the birds were 
not discovered until several 


TYPES OF 


hours (3-8) after death. 


RESULTS 


DIFFERENT 


Group 1. The results for 
this group of 20 animals, which 


received homologous bone 


ACH OF 


marrow after being subjected 
to 2,500 rads total-body irra- 


rlutination, 


diation, appeared to fall into 


‘ 


two distinct categories: 
Birds Which Died within 10 
Days after Trradiation.—Eleven 


, SHOWING 


e ol ag 


> 


Groupe 2 


of the 20 birds receiving homol- 
ogous bone marrow died be- 
tween days 5and 10. Of these, 


ra increase the de 


7 were from Series A and 4 


MARROW 


from Series B. These pigeons 
showed a rapid decrease in 
packed cell volume of both 


individuals in the 


erythrocytes and white cells. 


This was accompanied by se- 


on all the 
tled to the 


vere weight loss, and, in most 


ESULTS OBTAINED 


cases, both nausea and diar- 


R 
data 


rhea preceded death! There 
Was no indication of erythro- 
cyte repopulation by the mar- 
row transplant in any of these 
animals (for typical example 
see bird No. 3095.20, Table 1). 

Birds Which Survived Longer 
than 15 Days.—Ot the 3. birds 
in Series A which survived 10 


with 


7. 
I 
=< 
= 
7. 
/ 
Z. 


serum 
glutination 


pb Prat 
yf 
In VIVO: 


( 


days, | died on day 16, 1 on 
day 21, and 1 on day 30. The 
6 birds of Series B which sur- 


Isotonic saline 


gelutination”’ 
activity 


IRRADIATE 
livia 


risoria 


vived beyond the tenth day 


in Vivo agglutination 
ictivity 
i 


“in vivo : 


min 


died subsequently on days 16, 
21, 26, 29, 50, and 78. Inall 
birds which survived more 
than 10 days there was defi- 
nite evidence of erythrocyte 


rlutinin 
rrees of agglutination san 


C activity wi 


DATA FOR 3 OF 


] Jegree ol 
itination in these samples is due to ‘ 


repopulation by the donormar 
row. That is, the red blood 
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cells of the 9 survivors were 
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agglutinated by antisera specific for the donor’s blood type, anti-B and/or anti-C 
(see birds Nos. 3090.40, 3095.13, and 3095.29 in Table 1). The earliest evid- 
ence of repopulation appeared 8 days after transplantation (see No. 3095.13, Table 
1), though in some birds repopulation was not demonstrated until day 12. 

In the test with host-specifie agglutinin (Sieva extract), a definite reduction in 
activity occurred, paralleled by an il creased activity with donor-specific antisera 
beginning with the sample obtained on day 8 (see birds Nos. 3090.40, 3095.13, and 
3095.29 in Table 1). The erythrocytes of all 4 birds which survived 28 days or 
more were ho longer agglutinated by host-specific agglutinin but were, however, 
completely reactive with donor-specific antisera, indicating total replacement of 
erythropoiesis by the transplant marrow (see birds Nos. 3090.40 and 3095.29, 
Table 1). 

Group 2. The 20 animals which received 2,500 rads total-body irradiation fol- 
lowed by heterologous bone marrow transplantation can be separated into 2 classes. 
The serological data for 3 of the birds in this group which show each of the different 
types of results obtained are shown in Table 2. 

Birds Which Died Prior to Sampling on Day 6.—Ot the 12 birds which died be- 
tween days ! and 6, 4 were from Series A and 8 from Series B. Hematocrit read- 
ings on day 4+ revealed that 11 of these birds had subnormal erythrocyte and white 
cell volumes. The other bird died prior to sampling on day 4+. These 12 pigeons 
manifested weight loss, ind in most cases symptoms of nausea and diarrhea were 
noted prior to death. 

The blood sample obtained on day 6 from bird No. 3090.27 (Series A) appeared to 
be agglutinated as it dropped from the wing vein and ran down the side of the col- 
lection tube (see Table 2). Second and third samples, taken from this pigeon, | 
and 2 hours later, confirmed that the erythrocytes were definitely agglutinated. 
This evidence indicates that such agglutination occurred in vivo: hence the term 
“in vivo agglutination” would apply. 

The day-4 samples from birds in Series A were not examined for “in vivo agglu- 
tination,”’ since such a reaction had not been anticipated. All 7 birds of Series B, 
which died prior to the day-6 bleeding, yielded blood samples on day 4 showing ‘‘in 
vivo agglutination.”? Serum samples from these birds on day 0 and 4 (with one 


exception) contained no demonstrable circulating agglutinin in excess of that im 
the agglutinated cell mass. Bird No. 3095.38 (see Table 2), whose day-4 sample 


Was strongly agglutinated “in vivo,” yielded serum on that day which showed defi- 
nite agglutination with blood samples from 9 pigeons, although it was negative 
with 7 dove samples. Erythroeyte suspensions of the day-O0 samples of the host, 
pigeon number 3095.38, and the bone marrow donor, dove No. D.'T. 20, were in- 
cluded. 

Birds Which Died between Days 6 and 14.—Seven birds died during this period 
(4 on day 6, 1 on day 7, and 2 on day 8). All but | bird (which died before the 
phenomenon had been discovered) in Series A showed strong “in vivo agglutina- 
tion” at the time of day-6 bleeding. 

The sera of 2 birds in Series A contained an agglutinin which was specific for 
pigeon erythrocytes (host species) i contrast to dove red blood cells (donor species). 
“Tn vivo agglutination” was obvious on day 4 and appeared to be more pronounced 
on days 6 and 8 in one bird in Series B which died on day 8. On days 4, 6, and 8 





28 GENETICS; SHAW AND VERMUND Proc. N. A. 8S. 


circulating agglutinins specific for pigeon red cells (in contrast to those of doves) 
were present, while the day O serum samples revealed no agglutinin activity for 
either pigeon or dove cells. 

Bird No. 3095.37 of Series B (Table 2), whose day-O sample was entirely normal, 
survived until day 13 and showed “‘in vivo agglutination” beginning with the day-4 
sample. The sample on day 6 appeared to be more completely agglutinated, while 
the degree of “in vivo agglutination” seemed to be reduced in the day-8 sample and 
not detectable (even microscopically) on days 10 and 12. The day-4 sample was 
specifically reactive with pigeon red cells, even those obtained from this host on 
day 0. This host-species-specific hemagglutinin activity was stronger on day 6 and 
remained apparently constant for the 8-, 10-, and 12-day serum samples, in spite 
of the absence of any indications of “in vivo agglutination” on days 10 and 12. A 
test of the blood samples from this bird on days 10 and 12 revealed that the ery- 
throcytes were completely reactive with dove-specific antisera. Tests with pigeon- 
specific antisera were negative, even under microscopic examination, indicating 
total replacement of erythropoiesis by the transplanted dove bone marrow. 

Group 3. All 20 birds in this group died between days 5 and 17, with an 
average survival time of approximately 8 days. These birds manifested increasing 
anemia, weight loss, nausea, and diarrhea. None of the blood samples from birds 
in Series B showed “‘in vivo agglutination” or agglutinin activity with erythrocytes 
of livia or risoria. 

Group 4. The 20 birds in this group showed no anemia, emaciation, nausea, or 
diarrhea and are still apparently normal in all respects at the time of this report. 
As in Group 3, neither “in vivo agglutination” nor agglutinin activity was demon- 
strable. 


DISCUSSION 


The results of this preliminary work indicate that bone marrow transplantation 
alone may not be a satisfactory protective treatment for acute lethal irradiation 
damage in the species C. livia. In fact, those birds treated with heterologous mar- 
row died slightly earlier, on the average, than those lethally irradiated birds receiv- 
ing only saline. Though 6 of 20 of the irradiated pigeons which were injected with 
homologous marrow survived considerably longer than did irradiated controls, 
none has approached the survival time obtained with mice.*: * > 1% 

Lack of prolonged survival in livia was possibly due to transplantation of inade- 
quate amounts of bone marrow, though this is probably not the only contributing 
factor. Intestinal damage may prove to be more serious in C. livia than has been 
suggested for mice and other mammals.'! The fact that a large percentage of the 
irradiated pigeons which received homologous transplants died within the first 7 


days, showing symptoms of severe nausea accompanied by diarrhea, could be in- 


dicative of gastrointestinal damage. This, coupled with the observation that a sub- 
stantial number of birds died during a period after detectable erythrocyte repopu- 
lution and increasing red cell count, does not suggest anemia as the cause of death. 
The symptoms associated with these deaths appear to be more closely allied with 
those reported for rabbits, which Porter and Murray"! attributed to bacterial and 
viral infection. These investigators have demonstrated that treatment with anti- 
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bioties in addition to the bone marrow transplant increases survival time for 
lethally irradiated rabbits. 

Modification of the dove bone marrow transplantation technique sufficiently to 
obtain therapeutic effects is complicated by the hemagglutinin produced in pigeons 
so treated. This host-species-specific agglutinin appears very rapidly and is suffi- 
ciently potent to result in a very high degree of “in vivo agglutination,” and this, 
superimposed on the cessation of host erythropoiesis, may prove to be therapeu- 
tically insurmountable. Characterization of this hemagglutinin from the present 
data partially elucidates the mechanism involved. Since no agglutinin could be 
demonstrated in either the host’s or the donor’s circulation prior to irradiation and 
transplantation, it seems logical to assume that this agglutinin was either produced 
or released after treatment. lIurthermore, the fact that such hemagglutinins 
were not detectable in either the saline-treated or the homologously treated birds 
tends to rule out the possibility of viral or bacterial hemagglutinins. For the same 
reason, it is highly unlikely that these agglutinins are radiation-induced. The most 
plausible interpretation is that these agglutinins were elaborated by the cells of the 
heterotransplants. 

Whether this agglutinin is an antibody or a complex of antibodies elicited by anti- 
genic components of the host remains unanswered. ‘The present data show that 
the specificity of the “agglutinin” is comparable to the specificity anticipated if 
risoria were immunized with the cells of /ivza. As far as has been investigated, 
these hemagglutinins meet all requirements of an antibody produced by the foreign 
graft against an antigen or antigens of the host. 

The fact that this reaction occurs very early following the transplantation and is 
of such magnitude is possibly only a reflection of the high metabolic rate in avian 


species. Here we may also be dealing with bone marrow which contains a relatively 


high proportion of immunologically active cells, since lymph nodes, comparable to 
those present in mammals, are presumably absent in avian species. 


SUMMARY 


1. Sexually mature pigeons (C. livia) were subjected to 2,500 rads total-body 
irradiation and subsequently injected intravenously with: (1) homologous bone mar- 
row (backcross hybrids possessing erythrocytic antigens, Band or C, derived from 
C’. guinea); (2) heterologous bone marrow (ring dove, S. risorza); or (3) sterile saline 
(irradiated controls. ) 

10 days showed definite evidence of erythrocyte repopulation (initial appearance 
of donor-derived red cells was between days 8 and 12). ‘Total replacement of 


2. The 9 out of 20 birds of the homotransplant group which survived more than 


erythropoiesis by the homotransplant was demonstrable on day 28 in 4 out of 20 
of the birds so treated. These birds died subsequently without any indication of 
chimera reversal, but prolonged survival comparable to that observed in mice was 
not obtained, 

3. Heterotransplantation resulted in the rapid appearance (4-6 days) of host- 
species-specific hemagglutinins sufficiently potent to produce a high degree of “in 
vivo agglutination.”” The most plausible interpretation is that such agglutinins 

I 


were elaborated by the cells of the heterotransplant. 





GENETICS: SZILARD Proc. N. A. 8. 


* From the Department of Genetics (paper No. 712) and Department of Radiology, University 
of Wisconsin. This project was supported in part by the Research Committee of the Graduate 
School from funds supplied by the Wisconsin Alumni Research Foundation and in part by a re- 
search grant and C-4051(R1) to H. Vermund) (G-3884 to M.R. Irwin) grants from the Depart 
ment of Health, Education, and Welfare, of the National Institutes of Health, Publie Health 
Service 

'C.C. Congdon, Delta Uphoff, and Egon Lorenz, J. Natl. Cancer Inst., 13, 73-107, 1952. 

71D. W. H. Barnes, C. FE. Ford, P. L. T. Ilbery, and J. F. Loutit, Transpl. Bull., 5, 101-106, 
1958 

1d. W. van Bekkum and O. Vos, J. Cellular Comp. Physiol., 50, 139-156, 1957. 
C. kk. Ford, J. L. Hammerton, D. W. H. Barnes, and J. F. Loutit, Nature, 177, 452-454, 


N. Gengozian, T. Makinodan, C. C. Congdon, and R. D. Owen, these PROCEEDINGS, 44, 

565, 1958 

M.. R. Irwin, L. J. Cole, and C. ID. Gordon, J Exptl. Zoél., 73, 285-308, 1936. 

1). L. Lindsley, T. T. Odell, and F. G. Tausche, Proc. Soc. Exptl. Biol. Med., 90, 512-515, 1955 
’'T. Makinodan, Proc. Soc Exptl Biol. Med., 92, 174-179, 1956. 

’W. J. Miller and C. R. Bryan, these PROCEEDINGS, 39, 407-412, 1953. 

T. T. Odell, F. G. Tauseche, D. L. Lindsley, and R. D. Owen, Ann. N.Y. Acad. Se7., 64, 811 
, 1957 

K. A. Porter, and J. I. Murray, J. Natl. Cancer Inst., 20, 189-206, 1958. 

> 1D. H. Shaw, and J. H. Stimpfling, unpublished material. 

J.J. Trentin, Proc. Soc. Exptl. Biol. Med., 96, 139-144, 1957. 


ON THE NATURE OF THE AGING PROCESS 
By Leo Szitarp* 
ENRICO FERMI INSTITUTE OF NUCLEAR STUDIES, UNIVERSITY OF CHICAGO 


Communicated by Theodore Shedlovsky, November 24, 1958 


Introduction. This paper represents an attempt to describe a hypothetical bio- 
logical process that could conceivably account for the phenomenon of aging. Aging 
manifests itself in much the same general manner in all mammals, and we are in a 
position to learn enough about the aging of mammals to be able to test the validity 
of a theory that leads to predictions of a quantitative kind——-as does the theory 
here presented. 

We know that a gene can be responsible for the synthesis of a specific protein 
molecule, which in many cases has a known enzymatic activity. When we speak 
later of a mutant, or “incompetent,” form of a gene, we mean an altered form of the 
gene, Which cannot synthesize the specific protein molecule in its chemically active 
form. 

Our theory assumes that the elementary step in the process of aging Is an “aging 


hit,’ which “destroys” a chromosome of the somatic cell, in the sense that it renders 


all genes carried by that chromosome inactive. The “hit”? need not destroy the 


chromosome in 2 physical sense. (See note | added in proot.) 

We assume that the “aging hits” are random events and that the probability 
that a chromosome of a somatic cell suffers such a ‘hit’ per unit time remains con- 
stant throughout life. We further assume that the rate at which chromosomes of 
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a somatic cell suffer such “hits” is a characteristic of the species and does not vary 
appreciably from individual to individual. 

As a result of an aging process of this nature, the number of the somatie cells of 
an individual organism which have “survived” up to a given age (in the sense of 
remaining able to fulfil their funetion in the organism) decreases with age. On the 
basis of our assumptions, spelled out below, the “surviving” fraction of the somatic 
cells decreases with age at an accelerating rate. 

Our theory postulates that when f, the surviving fraction of the somatic cells of an 
individual, approaches a certain critical vaiue, f*, then the probability that. that 
individual may die within a period of one year will come close to 1. On this 
basis, the theory establishes a relationship between the surviving fraction of the 
somatic cells and the age of death of the individual. 

Because the young mammalian organism may be assumed to have a large func- 
tional reserve, we shall assume that the surviving fraction of the somatic cells of an 
individual may fall substantially before the organism loses its capacity to live, per- 
haps toa value somewhere between |; and !/ jp. 

The precise meaning of the term ‘critical value,” f*, will shift as we go from the 
crudest form of the theory, which we shall discuss first, to a less crude form of the 
theory, which we shall discuss thereafter. In the crudest form of the theory, we 
shall assume that an adult does not die of natural causes until the surviving fraction 
of his somatic cells comes very close to the critical fraction /* and that he dies at 
the critical age, i.e., within the year in which this surviving fraction reaches the 
critical fraction f*. Thus, in its erudest form, the theory postulates that the age 
at death is uniquely determined by the genetic makeup of the individual. 

Clearly, this cannot be strictly true, for, if it were true, identical twins would die 
within one vear of each other. In fact, the mean difference between the ages at 
death of female identical twins ean be estimated to be about 3.5 years. The 
discrepancy arises from the failure of the crude theory to take into account that in a 
cohort of identical individuals the number of deaths per year may be expected to 


rise as a continuous function with advancing age and that an appreciable number 


of deaths may be expected to occur at ages lower than the “critical age. 

If not otherwise stated, our discussion here relates to man and, in particular, to 
the female of the species. In the case of man, the somatie cells of the female con- 
tain m 23 pairs of homologous chromosomes. There may be in man perhaps 
15,000 genes. There may be a larger number of specific DNA molecules which 
are inherited from generation to generation, but we designate as ‘‘genes’”’ here 
only those DNA molecules which would handicap the individual if present in a mu- 
tant form. An individual who is a heterozygote for a mutant gene might not 
necessarily be handicapped under the conditions prevailing at present in the United 
States, where essentially no adult dies for lack of food or shelter and no adult has 
a reduced propensity to procreate because of his inability to proy ide food or shelter 
for his offspring. But such a heterozygote would have been handicapped (nccord- 
ing to our definition of the term “gene’’) under conditions which were prevalent in 
the past-——up to recent times. The present abundance of mutant forms of genes 
in the population may not correspond to the steady state under present conditions. 

We may assuine that the “genes” somehow affect differentiation and morpho- 
genesis during the embryonic development of the individual and that 2 mutant form 
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of a gene may cause, with a certain probability——appreciable even in the hetero- 
zygote—a developmental abnormality of the individual. 

We assume that among the 15,000 genes, there are perhaps 3,000 genes which 
are important for the functioning of the somatic cells of the adult. We shall call 
these genes “vegetative genes,” and a mutant form of such a gene we shall desig- 
nate as a “fault.” Of the remainder of the genes, we shall assume that they are 
irrelevant for the functioning of the somatic cells of the adult organism. 

We postulate that, in the course of aging, a somatic cell remains functional as 
long as, out of ‘each pair of homologous vegetative genes, at least one of the two 
genes is competent and active and that the cell ceases to be functional when both 
genes are out of action. Accordingly when a chromosome suffers an aging hit, the 
cell will cease to be functional if the homologous chromosome has either previously 
suffered an aging hit or if it carries a fault. 

According to the views here adopted, the main reason why some adults live 
shorter lives and others live longer is the difference in the number of faults they 
have inherited. If we assume that faults are distributed in the population at ran- 
dom, then we can compute the distribution of the faults, from the mean value of 
faults per person (which we shall designate by n). From the observed distribution 
of the ages at death, between seventy and ninety years of age, we shall be led to 
conclude that we haven > 2. For n = 2 we would obtain from the crude theory for 
the critical surviving fraction of the somatie cells f* @ !'/y. For n t we would 
obtain f* ® 3/4. On this basis we shall be led to conclude that we have n < 4. 

We shall, for the purposes of our discussion, adopt, as a reasonable value, n 
and then we obtain f* & '/., which would seem to be a reasonable value. 

The “Surviving” Fraction of the Somatic Cells.—We shall now proceed to com- 


‘ 


pute the “surviving” fraction of the somatic cells of a female who has inherited r 
faults, asa funetion of her age. 
We designate by é the average number of “aging hits” that have been suffered by 


2 chromosome of a somatic cell, and we may write 


1 age (1) 
2m 
where 7 is the average time interval between two subsequent aging hits suffered 
in toto by the m pairs of homologous chromosomes contained in a somatic cell. 
We may call this average time interval 7 the “basic time interval of the aging 
process.” 
Let us now consider a female who has inherited r faults. If none of the pairs of 
homologous chromosomes contain more than one fault—a condition likely to be 
fulfilled if r is small compared to m—then we may write for the “surviving” frac- 


tion of her somatic cells at a given age 


i {1 (1 


In f = (m — r) In [I 


lor & < | we may write, from equation (3), neglecting mé! and r&*, ete., 
o | ? tan] s s 


' 
In — m(& 
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Writing 


we may write 
l - 7 
In ~ = m[— + p]-[1 — ( 
/ 


a 


provided r < « < 2m (i.e., p < 


| = 
In — = m7?(1 — 7), 
i 


where yn = &+ p. In place of equation (7) we may write, in our approximation, 


We may also write inversely 


or, expanding, 


\, In 4 2 (10) 
” J é . 


According to the assumption of the crude theory, f, the surviving fraction of the 
somatic cells, reaches the critical value f* at the age of death, which we designate 
by t,. Further, we designate by x, the average number of aging hits suffered in 
toto, up to the age of death, by the m pairs of homologous chromosomes of the 
somatie cells. ‘Thus we have, at the age of death, 

fr 2mé 
e 


it 
(12} 
2m 


Accordingly, we may write at the age of death, where we have f = /*, from equation 


(8), 
l (vr, + r)? uetr 
In i ; 
Ds tm 2m 


Similarly, we may write at the age of death, from equation (10), 


| 


| | 
utr f4m In — + In — 
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Bons. l l e 

+r tm In — + In —, (15) 
T * a 

For the genetically perfect female, for whom we have r = 0, we shall designate the 
age at death by tf. We shall call ty the “life-span”’ of the species. 

From equation (15) we may write for, the life-span, , 


i. 


(17) 


(18) 


As may be seen from equation (17), the addition of one fault to the genetic makeup 
of an individual shortens the life of that individual by At = 7, so that we may write 


At per fault = r. 19) 


This expresses one of the basic results of our theory. According to equation (19), 
an individual whose genetic makeup contains one fault more than another individual 
has a life-expectancy which is shorter by 7, the basic time interval of the aging proe- 
ess. This holds true within the crude theory for individuals who have inherited a 
small number of faults. 

Concerning the life-span, tf), we may write, from equations (11), (13), and (16), 


l l /to\? 1 to 
In — = ( y( = ), (20) 
7 tm \r 2m 7 


and, from equations (15) and (16), we may write 


lo ben lh Sa | 
: . ee 
T \ ‘i \ i 


The Distribution of the Ages at Death.—The above equations hold within the 
framework of the crude form of the theory. In this form of the theory, members 
of one cohort would die only in certain years—at the critical ages, t-—-and thus 
the years in which death occurs within one cohort would be separated from each 


(21) 


other by time intervals of 7 years; no deaths would occur in the intervening yeurs. 
Further, if the distribution of the faults in the population is random, then the 
?,, occurring at each age, is given by the Poisson distribution: 


number of deaths, / 


ne n on 
P; (22) 
7! 


— ft, : 
and where n stands for the 


where, according to equation (18), we have r 


average number of faults per individual. 
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The distribution of the ages at death in the population is actually a continuous 
function of the age. Even though the probability that an individual may die 
Within a year may increase rather steeply as the surviving fraction of his somatic 
cells approaches the critical value f*, genetically identical individuals do not all die 
at the same age. The observed mean age difference at death of identical twins 
may be regarded as a measure of the scattering of the ages at death, which is left out 
of account by the crude form of the theory and to which we shall refer as the “‘non- 
genetic scattering.” 

For the time being, we shall continue to leave this non-genetic scattering out of 
account; yet, for the sake of convenience, we shall henceforth describe the distribu- 
tion of the ages at death by P(r), a continuous function of r, in place of the diseon- 
tinuous “‘Poisson”’ values, P;. For P(r) we may write 


? 
i] 


P(r) 


(r+1) 


where [ represents the gamma function (which for integral values of r + 1 assumes 
the values of r!) and where we have 


(24) 


For the number of deaths occurring within a cohort per unit time, we may then 


write, according to our theory, 
d(theor. ) . (25) 


From equation (24) we obtain 


Thus we may write, from equation (25), 


d(theor. ) " per year, 


where r is given by equation (24) and where 7 is expressed in years. 
The approximation used throughout this paper holds for small values of r, which 


correspond to high ages at death. We may say that at high ages of death the dis- 


tribution of ages at death in the population is represented by a reversed Poisson 
distribution (27), where small values of r correspond to high ages at death. 

Lower Limit for n.—-We shall now proceed to compare the distribution of the 
ages at death, as given by our formula (27), with the actually observed distribution 
of the ages at death, as given by the U.S. Life Tables, based on the 1949-50 Census. 

For the purposes of this comparison, we shall use Table 6 for white females, which 
lists the number of deaths per year, in yearly intervals, as a function of age. Ac- 
cording to this table, the maximal number of deaths occurs between the eightieth 
and eighty-first year; the corresponding maximal number of deaths per year is 
0).0344 per person. 
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The distribution of the ages at death is not symmetrical around the age at the 
maximum, /* = 80.5 years of age; the nuimber of deaths per year fall faster toward 
higher ages than toward lower ages. Thus the table lists, for the number of deaths 
per year, 0.0230 per person between the ages seventy and seventy-one and 0.0179 
per person between the ages of ninety and ninety-one. 

We may derive from this table a “normalized” distribution of the ages at death by 
forming /(obs.), the ratio of the number of deaths per year and the maximal number 
of deaths per year, 0.0344. Thus we obtain R(obs.) = 0.667 at 70.5 years of age; 
Robs.) 1 at 80.5 years of age; R(obs.) = 0.520 at 90.5 vears of age. 

We may similarly obtain from the number of deaths per year, given as a function 
of age by the theory, a “normalized” distribution of the ages at death, by forming 
R(theor.), the ratio of the number of deaths per year given by equation (27) and 


the maximal number of deaths per year given by d(theor.) jax. : 


{ 
d(theor.) max . (28) 
1 ; 


. 
(y+1)/ max 


If we designate by r* the value of r for which this expression becomes maximal, 
we may write, forn 2 2, 
r an — 0.5. (29) 
Accordingly, we may write 


(30) 


d(theor.) n Dieak 


R(theor.) 


d(theor.) max mas 14) 

We may now ask for what value of n would the normalized Poisson distribution, 
R(theor.), fit R(obs.), both above and below 80.5 years of age so that we have for a 
suitably chosen value of Ar, for r= r* + Ar, R(theor.) = 0.667 (the value of R(obs.) 
at 70.5 years of age) and that we also have, for r = r* — Ar, R(theor.) 0.520 
(the value of R(obs.) at 90.5 years of age). 

It turns out that such a fit is possible only for a value of n which is very close to 
n 2. For the corresponding value of Ar we obtain Ar 1.4. For the corre- 


sponding value of 7 we may write 


10 
years. 


Ar 


Forn 2 and with Ar 1.4, we obtain 7 = 7.15 years. 

For values of n which are substantially larger than 2, it is not possible to fit the 
normalized Poisson distribution /(theor.) to R(obs.) in this manner. If R(theor.) 
is made equal to 0.520 (the value of R(obs.) at 90.5 years of age) for r = r* — Ar, 
then, for r r* + Ar, we have R(theor.) < 0.667 (the value of R(obs.) at 70.5 
years of age). 

Because there is reason to believe that, below 80.5 years of age, the crude theory 
gives too low values for R(theor.), we cannot exclude the possibility that we have 
n> 2. Therefore, from the fact that R(theor.) derived from the crude theory fits 
Robs.) for n 2 between the ages of 70.5 and 90.5, we may not conclude that we 


actually have n = 2, and we may only conclude that we have 
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n22 (32) 


Approximation of the Poisson Distribution by a Gaussian Relationship between r 
and n.—R (theor.) given by equation (30) goes over into a Gaussian for n > 1. 
For a Gaussian, the value of R(theor.) = 0.667 (the value of R(obs.) at 70.5 years 
of age) corresponds to a distance from the maximum of 0.9 o, where a is the stand- 
ard deviation of the Gaussian. Similarly, R(theor.) = 0.520 (the value of R(obs.) 
at 90.5 years of age) corresponds to a distance from the maximum of 1.140. Thus 
the time interval of 20 years around the maximum corresponds to 2.040, and hence 


we have 
20 
years 9S years 
2.04 * ‘ 
Because the variance of a Poisson distribution is given by its mean, n, we may 


write (forn > 1) 
(34) 


9.8 
and thus we obtain 7 years. 
V 


nM 

While equation (34) holds, strictly speaking, only for large values of n, the error 
is small even forn = 2. 

For n = 2 from equation (35) we obtain 7 = 6.82 in place of the previously given 
value of 7 = 7.15. 

For n = 2.5 from (55) we obtain 7 = 6.2 years in place of the “correct”’ value of 7 

6.3 years, which we find by fitting the “normalized” Poisson distribution, as 
well as possible, to R(obs.). 

Thus, for most of our purposes we may use equation (35) for values of n for n 2 2. 
From equation (35) we obtain 

NT oV Nn, (36) 
where nz represents the average life-shortening caused by the “load of faults,” 
n. Since the value of o is empirically fixed, the higher we assume n to be, the 
higher is the lite-shortening effect which we must attribute to it. In this sense, the 
life-shortening effect of the mutation load increases with Yn. 

Correction of + for the Non-genetic Scattering of the Ages at Death.—Because the 
non-genetic scattering of the ages at death has so far not been taken into account 
by us, the observed distribution of the ages at death may be expected to be actually 
somewhat broader, and, accordingly, the actual value of 7 may be expected to be 
somewhat lower than the values given above. 

The mean age difference at death between female identical twins has been re- 
ported by Franz J. Kallman to be about 2.6 years for twins dying above the age 
of 60. From this value we may estimate, on the basis of the Life Tables, the mean 
age difference at death of female identical twins who die as adults above the age 
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of 40 to be 3.4 years. If the distribution of the ages at death of genetically identi- 
eal individuals resembled a Gaussian, then the variance of the distribution of the 
ages at death in the population would be equal to the sum of the variance of this 
Gaussian and that of our theoretical distribution of the ages at death. By making 
such an assumption, for the purposes of this computation, we may then correct 
the values of 7, given above, as follows: 

From the fact that the mean age difference at death of female identical twins may 
be taken to be about 3.4 years, it follows that the standard deviation of the distribu- 
tion of their ages at death is about 3 years. Using this value, we find that the non- 
genetic scattering increases the variance of the distribution of the ages at death by 
a factor of about 1.1 and, accordingly, the previously given values of 7 must be 
reduced by 5 per cent. 

Thus we may now write, for the corrected values of 7, for n = 2, 6.8 vears; 
and for n = 2.5, 7 = 6 years. 

We may also write on this basis--within the limits of the approximation—for 
n> Zz, 

9.3 
— years. 
Vn 
Substituting this value of 7 in equation (16), we obtain 


to t, 
Wn =a 7. 38) 


T 9.5 

The Value of the Critical Surviving Fraction of the Somatic Cells f*—-Upper Limit 
for n.—-In order to compute the critical surviving fraction of the somatic cells, {*, 
we shall now make use of the fact that (for white females) the maximal number of 
deaths per year occurs at 80.5 years of age. Our theory demands (29) that the 
maximal number of deaths per year should occur for individuals for whom we have 
r= n— 0.5. Accordingly, we may write t, = 80.5 and r = n — 0.5. We thus 
obtain, from equation (38), 


and, from equation (20), 


l Ll #805: .; NA 1 /80.5 J 
In : VYntn— 0.0 | — Vnt+n— 0.9 ; (40) 
i fi tm \ 9.3 iF 2m \ 9.8 


From this equation we may now compute for a given value of x the corresponding 
value of f*. Thus we obtain for n = 2, f* & '/,:forn = 2.5, f* & '/,:forn = 4, 


f*¥Q} 


| 


On this basis we may then write, by assuming f* > !/j», 


(ie eae (41) 


A value of f* & '/5 would seem to be rather reasonable and, therefore, we shall 


adopt, as a reasonable value for n, the value of 


n A 2.5, (42) 
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and, as a reasonable value for 7, the corresponding value of 

tT = 6 years. 3) 
We shall in the remainder of the paper base all our discussions on these values of 
nand rv. We cannot exclude, of course, the possibility that n might be somewhat 
larger and that 7 might be somewhat lower. 

The ‘Physiological Age.’’—The general physiological age may be defined for a 
given population on the basis of its age-specific death rate; according to our theory, 
it may be defined as the age of the genetically perfect female who has the same 
surviving fraction of the somatic cells, f. Accordingly, we may say that two fe- 
males, whose genetic makeup differs by A faults, differ from each other in their 
physiological age by 7A years at sufficiently high ages as demanded by the approxi- 
mation used. 

Changing the Load of Faults.—If, as a result of living under ‘‘modern” conditions, 
our load of faults should, in time, be doubled, then the average adult woman would 
live nr years less than she does today. 

For n = 2.5 we have nr = 15 years. Thus the physiological age of the average 
female at 65 would be the same as that of the average 80-year-old woman today. 


~ 


If we were to assume that n > 2.5, then nz would amount to more than 15 years 


because nz increases, according to equation (36), with /n. 
5 | 


A doubling of our load of faults might conceivably occur, in time, through the 
exposure of the population, generation after generation, to ionizing radiation, in an 
intensity that doubles the mutation rate. 

Such an increase in our load of faults might perhaps occur also as a result of the 
current practice of controlling the family size. As spelled out below this practice 
might conceivably eliminate one of the selection pressures which have tended to 
keep our load of faults low. 

We may, on this occasion, also ask how much advantage the genetically perfect 
(faultless) female would have over the average female of today. Assuming n = 
2.5, we may say, on the basis of considerations similar to those presented above, 
that the genetically perfect female would at 50 years of age have the same phys- 
iological age as the average female of 35 today. Her most probable age at death 
would be 92 instead of 80. If n were larger than 2.5, the advantage of the geneti- 
cally perfect female would be greater. 

Life-Shortening Effect of Ionizing Radiation on the Adult Offspring of the Exposed 
Population.—Experiments of W. L. Russell have shown that the offspring of mice 
which have been exposed to a dose of fast neutrons have a reduced life-expectancy. 
This has generally been interpreted to mean that exposure of the parents to ionizing 
radiation induces mutations in the germ cells of the gonads and thus “reduces the 
viability” of the offspring. 

From the point of view of our theory, however, we have to distinguish between 
that reduction of the life-expectancy of the offspring which is due to an increased 
mortality of young animals and that reduction which is due to a decrease in the life- 
expectancy of the adults. All the mutations induced by ionizing radiation may 
contribute to the former, but only the “faults” contribute to the latter. 

11 the case of man, at least, it should be possible to make a fairly clean separation 
between these two categories of life-shortening. In the case of man the U.S. Life 
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Tables show that the number of deaths per year falls, from an initial high value in 
the first year of life, to about 40 per 100,000 per year at the age of 10. Moreover, 
of the few deaths occurring at 10 years of age, a substantial fraction is due to acci- 
dents. Thus we are led to believe that, in the heterozygous individual, mutant 
(incompetent) forms of genes may cause the death of the embryo, or of the im‘ant 
below 10 years of age, while they do not cause death with an appreciable probability 
after the tenth year of age, unless they represent “faults.”’ Faults increase the age- 
specific death rate above 10 years of age only in conjuction with aging hits, and they 
increase it appreciably only above 40. 

If we observe the life-shortening of the adult animal in the offspring of an ir- 
radiated population, resulting from the induction of faults by ionizing radiation in 
one species, we may be able to predict, on the basis of our theory, the life-shortening 
for another species. This may be seen as follows: 

We obtain from equations (19) and (21) for the relationship between the life- 
shortening, At per fault and for the life-span of the species, fo, 


\” In re + In pe ) (44) 


The right-hand side of the equation contains only the chromosome number m 


At per fault r ( 
to a Se 


and the critical value f*. Therefore, if two species of mammals have the same 
value of f* and the same chromosome number m, their life-shortening per fault 
caused by exposure of their parents to ionizing radiation amounts to the same frac- 
tion of their life-span. We may call the ratio the “specific life-shortening” of a 
fault. 

If the two species of mammals may be assumed to have also the same number of 
vegetative genes and if the sensitivity of their genes to the ionizing radiation em- 
ployed may be assumed to be the same, then the number of faults produced by a 
given dose of radiation will be the same for the two species. Thus, according to 
equation (44), the radiation exposure will shorten the lives of the two species by the 
sume fraction of the life-span. 


If m, the number of their chromosome pairs, is different for the two species, then 


the “specific life-shortening”’ will be larger for the species which has the smaller 
chromosome number. According to (44), the specific life-shortening increases 
about inversely with Ym. The number of chromosome pairs is 11 for the Chinese 
hamster and 39 for the dog. Therefore, according to (44), the specific life-shortening 
per fault induced may be expected for the Chinese hamster to be higher than for 
the dog by a factor of about 2. 

The mouse has m = 20 pairs of homologous chromosomes and we may therefore 
estimate the life-shortening of man from the life-shortening of the mouse and vice 
versa, by postulating that the life-shortening per rep in man and in the mouse 
umount to about the same fraction of their life-span. (Some authors believe that 
man is about twice as sensitive to X-rays as the mouse and, if they are correct, then 
our estimated value for the life-shortening of man would be low by a factor of about 
two.) 

W. L. Russell found that an X-ray dose of 300, as well as a dose of 600 rep, in- 
duces 25 & 10~° gene mutations per rep per locus in the spermatogonia of mice, — As- 
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suming 15,000 loci, we may conclude that an X-ray dose of 667 rep would induce 5 
mutations in the diploid offspring of exposed mice. If we assume that one-fifth 
of the genes are vegetative genes, then 667 rep of X-rays induce on the average one 
fault in the offspring. Assuming 7 = 6 years for man, we may thus expect in man 
a life-shortening of the adult offspring of 6* = 3.3 days per rep. 

The number of faults induced in the offspring per rep depends on the nature of 
the ionizing radiation, and may be assumed higher for fast neutrons than for X-rays. 
Also, in the case of X-rays, the number of faults induced might conceivably depend 
not only on the total dose but also on the dose rate and be lower for lower dose rates. 

The actual value of 6* ought to be determined experimentally, for the different 
kinds of ionizing radiations which are of interest, by direct observation of the life- 
shortening of the adult offspring. Experimental data so far available are inade- 
quate. 

If the average mutation rate per gene per generation is 1/ 60,000 and if we assume 
for V,, the number of vegetative genes, .V, 3000, and tor V, the total number of 
genes, V, = 15,000 then we obtain for 4;, the spontaneous mutation rate of the hap- 
loid set of vegetative genes, yu, 0.05, and for we, the total spontaneous mutation 
rate of all genes, uw. = 0.25. We shall use these values for the purposes of our dis- 
cussion below. 

The Life-Shortening Effect of Ionizing Radiation on the Exposed Population. 
We may expect ionizing radiation to produce gene mutations in the chromosomes 
of the somatic cells of an exposed individual and we shall assume that the sensitiv- 
itv of the genes of the somatic cells is about the same as that of the genes in the 
spermatogonia and the oogonia. Because a certain fraction of these mutations, per- 
haps one-fifth of them, uffect vegetative genes, faults are induced in the chromo- 
somes of the somatic cells of the exposed individual. It can be shown that an ex- 
posure to ionizing radiation which induces on the average one fault per somatic 
cell must reduce—on this score along—the life-expectancy of the exposed individual 
by about 7 years. If exposure to ionizing radiation had no other life-shortening 
effect, the life-shortening, 6, of the individual in the exposed population would be 
equal to 6*, the life-shortening of the adult offspring of the exposed population (see 
note 2 added in proof). 

Maternal Selection Pressure against Faults?—-It is conceivable that a woman who 
carries a particular fault in her genetic makeup ceases to be capable of bearing 
children 7 6 vears earlier, on the average, than her counterpart who lacks that 


particular fault. This is what one would expect on the assumption that the ter- 


mination of a woman’s reproductive period is determined by her physiological age 
if all other factors are equal. If this assumption is correct, then a powerful selection 
has operated in the past that has tended to keep the load of faults low. 

In the past, infant mortality was high, the birth rate was high, and women kept 
on having children until the end of their childbearing period. Clearly, the ma- 
ternal selection mentioned above is switched off when women have two or three 
children between the ages of 18 and 25 and then avoid having further children. 

If such a “maternal” selection was the predominating selection of the past, then 
we may expect that, when this selection ceases to operate, our load of faults may 
double, in time. As discussed before, senescence would then set in about 15 years 


earlier. 
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However, even if we assume the worst in this respect, our load of faults still can- 
not increase by more than 24; = 0.1 per generation. This means that—at the 
worst—it would take 25 generations for our load of faults to double. 

The effect of the “maternal” selection here discussed might be estimated as 
follows: 

Let us single out one vegetative gene. If a woman carries, as a heterozygote, 


this particular gene in a mutant form, her physiological age is 7 = 6 years higher 


than that of another women who does not carry this particular “fault”? but who 
is otherwise identical in her genetic makeup. We now assume that the “phys- 
iological age”’ sets the termination of the reproductive period, and we take for the 
“most probable duration” of the reproductive period 30 years. Thus the fault 
singled out “most probably” shortens the reproductive period by one-fifth of its 
length. The fertility of younger women is higher than that of older women. Near 
the end of the reproductive period, the (average) time interval between two suc- 
cessive pregnancies might be by a factor k > 1 (of perhaps 2 or 3) longer, than the 
interval between two successive pregnancies, averaged over the whole reproductive 
period. 

If we postulate that this “maternal” selection constitutes essentially the sole 
selection pressure against faults, then we may write for the mutation equilibrium 
of the fault, singled out: 


n 
9 Mi 


“Ni, 2X 5KN,’ - 


and, for the value of yu; 0.05 per generation, we thus obtain k = n. For n 


we have 
k meek (46) 


Under “natural” conditions, for young women the average time interval between 
two successive pregnancies might be about 1 year; toward the end of the repro- 
ductive period it might be about 5 years; and, averaged over the whole reproductive 
period of about 20 years, it might be about 2 years. These values correspond to 
k} = 2.5. This coincidence might, of course, turn out to be fortuitous. 

Refinements of the Theory—Specialized Vegetative Genes.—It appears likely that 
there exist genes which are not essential for the functioning of most of the somatic 
cells of the adult but each of which is essential for the functioning of one particular 
kind of specialized somatic cell. We shall call these genes 
genes, and mutant forms of such genes we shall call “specialized faults.” 


“specialized vegetative” 

We shall now single out, for the purposes of our discussion, specialized cells which 
synthesize a gene product, a particular enzyme, for instance, that serves not the 
needs of these cells themselves but rather those of the organism as a whole. Such 
specialized cells might in some cases fulfill their function in the organism by releasing 
the enzyme into the circulation. 

We may in general assume that the normal young person has a considerable re- 
serve of such enzymes, and we shall specifically assume, for the purposes of this 
particular discussion, that the maximal output of a normal young person is higher 
than the need of the organism by a factor of about 6. 

The maximal output of enzyme by such a specialized cell may be assumed to be 
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lower by a factor of '/. in the heterozygous individual, who has inherited a mutant 
form of the specialized vegetative gene. Further, it may be shown that, for an 
individual who has inherited n = 2.5 faults, the surviving fraction of the somatic 
cells is about one-third at 54 years of age. Since, in the heterozygote, the special- 
ized cells under discussion carry one additional fault, the surviving fraction of the 
specialized cells will reach one-third about 6 years earlier. Thus, at 48 years of 
age, the maximal output of the enzyme of the surviving specialized cells of the 
heterozygote will be lower, by a factor of about '/s, than for a normal young person. 

On this basis we may then expect. that around 50 years of age there may become 
manifest, in such heterozygotes, symptoms of disease due to the insufficiency of the 
output of one kind of specialized cell. The inheritance of diseases of this class may 
be expected to show a marked degree of dominance. 

Speaking more generally, we may expect to see in certain heterozygotes, late in 
life, narrowly circumscribed degenerative phenomena which are caused by special- 
ized faults they have inherited. 

The Number of “Segments” per Chromosome.—Instead of assuming that a whole 
chromosome is ‘‘destroyed” in one aging hit, we might choose to assume that the 
elementary step in the process of aging consists in the random destruction of one- 
half of a chromosome. The formulas given above then remain unchanged, except 
that we have to write 2m in place of m. As one may see from equation (40), we 
then obtain, for the same value of n, a higher value for f*. Thus for n 2.5 we 
obtain f* @!/;. Apart from this, the general character of the theory remains un- 
changed. 

However, one might ask at this point whether one could not generalize the theory, 
presented above, by assuming that each chromosome consists of g segments and 
that the elementary step in the process of aging consists in the random destruction 
of such segments, independently of each other. By choosing the value of g larger 
and larger, we might then graduaily change the character of the theory and might 
end up with a theory which postulates that the aging process consists in a sequence 
of gene mutations of the chromosomes of the somatie cell. 

A theory of this kind would, however, come up against difficulties, which are as 
follows: 

As may be seen from equation (40) (where we now have to write gm in place of 
m), for a fixed value of f*, n goes up roughly parallel with increasing g. A very 
large value of n might, however, be incompatible with the known fertility of con- 
sanguinous matings. 

Further, as we increase g, we would also increase the difference of the life-expee- 
tancy of the female and the male. The male of the species has only one X chromo- 
some, While the female has two. Let us disregard here the possibility that a sub- 
stantial piece of the X chromosome might be covered, in the male, by genes con- 
tained in the Y chromosome. Let us also assume, for the sake of argument, that 
f* has the same value for the male as it has for the female. On the basis of these 
assumptions, we may then identify the male, from the point of view of his life- 


expectancy, with a female who has suffered g aging hits, prior to birth. Aecord- 


ingly, we may expect the adult male to live a shorter time, by gr years, than the 
adult female. 
Actually, according to the 1949-50 Census, the maximal number of deaths for the 
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white male occurs between the ages of 77 and 78, i.e., three years earlier than for 
the white female. This difference is three years less than what we would expect on 
the basis of our theory, which assumes g 1 and which gives an estimate for 7 of 
r = 6 years. This discrepancy indicates that perhaps the value of {* is somewhat 
larger for the male than tor the female. 

Because of the possibility that this might be the case, we conceivably have 
" 2 (9 2 would mean that the elementary process of aging consists in the 
“destruction” of one-half of a chromosome rather than a whole chromosome). 

However, there is no reason to believe that {* may be very much larger for the 
male than for the female. Therefore, the observed small difference between the 
life-expectancy of the female and that of the male may rule out a modification of 
the theory that assumes g > 1. 

Experimental Test of the Theory.-The most stringent experimental test of the 
validity of our theory is likely to come from experiments in which one observes a 
reduction in the life-expectancy of the adult offspring of, say, an irradiated mouse 


population. Experiments of this sort are needed in order to determine the value 


of 6*. Experiments of this sort will also show whether among the different phe- 
nomena which generally accompany senescence, such as the graying of the hair, 
the loss of accommodation of the eye, etc., there are any which are determined by 
the general physiological age, defined on the basis of the age-specific death rate. 


Arrangements for experiments along these lines are now under discussion. 


Iam grateful to Dr. Gertrud Weiss and to Dr. N. Conant Webb, Jr., of the Medical 
School, University of Colorado, for detailed criticism of some of the computations 
involved, and to Samuel W. Greenhouse, of the National Institute of Mental Health, 
for eritical evaluation of some of the approximations used. For carrying out 
auxiliary computations I am also grateful to Mrs. Dorothy Lathrop, Office of 
Mathematical Research, National Institutes of Health, and to Mrs. Norma French, 
Section on Theoretical Statistics and Mathematics, Biometries Branch, National 
Institute of Mental Health. 


* The author wrote this paper while serving as consultant to the Basic Research Program, Na- 
tional Institute of Mental Health, National Institutes of Health, U.S. Public Health Service, 
Department of Health, Mducation, and Welfare 

NOTATIONS 

fis the fraction of the somatic cells which “survive,” in the sense of remaining functional up to : 
given age 

f* isthe “surviving” fraction of the somatic cells at the age of death. 

ris the number of inherited faults. 

p isthe number of inherited faults per chromosome. 

n is the average number of faults per person in the population. 

r is the basic time interval of the aging process, defined as the average time interval between two 
successive aging hits suffered by the chromosomes of the somatic cell. 

is the number of pairs of homologous chromosomes of the female of the species. 

is the number of aging hits suffered, on the average, by the chromosomes of the somatic cells up 

to the age of death, by an individual who has inherited r faults. 

is the age, at death, of an individual who has inherited r faults. 
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f, is the life-span of the species, defined as the age at death of the genetically perfect female who 


did not inherit any faults 

¢(theor.) is the number of deaths per vear that will oceur in a cohort, as given by the theory. 

d/(theor.)max. is the highest number of deaths per veal that will occur in a cohort, as given by the 
theory 

R(theor.) is the ratio of the number of deaths per vear and the maximal number of deaths per year 
in a cohort, as given by the t heory 

Robs.) is the ratio of the number of deaths per year and the maximal number of deaths per year 
in a cohort, as given by the U.S. Life Tables for white females. 

o is the standard deviation of the Gaussian that approximates the observed distribution of the 

number of deaths per year, between the ages of 70.5 and 90.5 . 

r/lois the specific life-shortening per fault of the species. 

6 is the life-shortening per rep for a population that has been exposed to ionizing radiation. 

5* is the life-shortening per rep of the adults in the offspring of a population that has been exposed 
to ionizing radiation 

wu; is the spontaneous mutation rate of the haploid set of vegetative genes per generation 

uris the tetal spontaneous mutation rate of all genes in the haploid set 

Nis the haploid number of vegetative genes of the species 

\, is the haploid number of all genes of the species. 


gis the postulated number of ‘‘segments”’ per chromosome 


NovTes ADDED IN PROO! 


1. When we say that an aging hit ‘‘destroys’’ a chromosome of the somatic cell, we mean that 
that chromosome has been rendered inactive as far as its vegetative functions are concerned, 1.e., 
the genes which the chromosome contains will fail to produce the corresponding gene products. 
The question whether the chromosome is inactivated in any other sense is left open for the present. 
Thus it is left open whether, if a cell containing an inactivated chromosome were to duplicate, the 
inactivated chromosome would or would not duplicate and whether or not it would remain inac- 
tive after such a duplication. One might, for instance, imagine that the chromosomes of the 
somatic cell contain NNA strands which fulfil a vegetative function in the somatic cell by produe- 
ing the specifie gene products but do not duplicate when the cell duplicates. Aging hits would 
then inactivate these vegetative ‘‘copies”’ rather than the genetic copies. The latter would dupli- 
cate when the cell duplicates and would then produce fresh vegetative copies. This is yust one of 
several assumptions which one may make concerning the nature of the aging hits. For the pres- 
ent, we are free to choose among several such ad hoc assumptions 

2. In the case of exposed animals it is conceivable that their life is shortened, not only through 
the induction of gene mutations in the chromosomes of their somatic cells by the ionizing radia- 
tion, but perhaps also through some other effects of the ionizing radiation on their somatie cells, 
which may involve the chromosomes or some other components of the cell. Among such effects 
might be the breakage of chromosomes which may lead to the loss of a chromosome. However, 
the theory here presented does not cover the life-shortening effeet of ionizing radiation which is 
due to causes other than the induction of gene mutations in the somatic cells of the chromosomes. 
Disregarding such other effects, the “surviving” fraction of the somatie cells of an exposed female 
may be computed on the basis of the faults induced in the chromosomes of her somatic cells by the 
lOnIZINE radiation For a genetically perfect female who is exposed to a dose of ionizing radiation 
which induces, on the average, p faults per somatic cell, we may write for the “surviving” fraction 


of somatic cells: 


for p/2m< 1 and p/r< »mav write, in analogy to 


at ae " C ») [1 


where ¢, is the age of death of a genetically perfect female exposed to a dose of ionizing radiation 
that induced an average ot p faults in the chromosomes of her somatic cells 





ON IDENTICAL VANISHING OF HOLOMORPHIC) FUNCTIONS IN 
SEVERAL COMPLEX VARIABLES* 


By 8. BocuNER 
MATHEMATICS DEPARTMENT, PRINCETON UNIVERSITY 
Communicated November 7, 1958 


We will make some comment on the following recent result of J. J. Kohn.! In the 
complex space C,,:(2:, .. . ; Za) Uf; f, are holomorphic in the unit ball and 
have continuous boundary values on the boundary, and if on the boundary the 
linear combination Z:f; + ... + 2,f, is 0, then the functions fi, ...,f, are identically 
(). 

1. We denote by D any bounded circular domain in C,, by B its boundary, and 
by By an open neighborhood in B. We do not assume that D contains the origin in 
(,, but only that Bo has a positive distance from it. For any integer m, m = n, 
, fm Which are holomorphic in D and have continuous 


we take functions fi, 
boundary values on By. We also take m polynomials P, P,, in n symbols 
C1, ..., €,, homogeneous and of a common degree and we denote the degree by r. 


Finally, we introduce the function 


m 


$(f;2) = D> Pfs, (2) (1) 
u=l 


in the 2n variables ¢, Gni Zi... 2, If for z ¢ D we make the specialization 
c, z,, v | n, then we obtain a function (2; z) which in D is analytie in 
the 2n real variables 71, ..., %y3 4 Ynys 2 = xX, + ty,, and has continuous 
boundary values on Bo. 

THEeoreM 1. (i) Jf the function ®(2; z) is 0 on Bo, then the function ®(¢; z) zs 
identically O in (€; 2). 

(ii) Hence, if the polynomials P; P,,, are linearly independent over constants, 
then f; f, are all identically 0. 

Proof: In the case considered by Kohn, P; = (1, ..., P, = €,. These poly- 
nomials are independent and his conclusion follows. F 
In our general case, let a = (a, a,) be a fixed point of D, not the origin in 

If we put z, = fa,, v = n, then we have 


m 


(ia; ta) = U > P,(a)f,(ta), (2) 
w=! 


and in the complex ¢-plane the latter sum is defined and holomorphic in a certain 
annulus with center at ¢ = 0. Because of our assumptions, there is a neighborhood 
Dy in D such that for each fixed a ¢ D, the sum in equation (2) is 0 on an are of the 
f-annulus. Hence the sum is identically 0 in ¢, and hence for z € Do, ®(Z; 2) isO0.  Bya 
uniqueness theorem of E. Cartan, the function ®(¢; 2) is identically 0 as claimed. 
2. In the case of the unit ball a statement can be made even if the polynomials 
P,, are arbitrary non-homogeneous, and the reader will readily verify that 
this statement subsumes Theorem 1, for such a D. 
THeoreM 2. Jf D is the unit ball, and Pi, ... , Pm are any polynomials in $1, 
. €n, then the vanishing of ®(z; z) on a neighborhood Bo of the boundary of D implies 


a representation 
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@(¢05z) = 1 — 141 — 2... — Sn2n) WE; 2) (3) 
in which W(¢; z) ts a polynomial in (¢,) with coefficients which are holomorphic in D. 
Proof. If we put 
u=l—ae-—...— 32; (4) 


then D is defined by u > O and B by 0. If now we substitute 


in &(2; z) we obtain an expression A (1, 32, ... , Zn} 21, ..., 2,) Which is a polynomial 
in U, Z,..., 2, With coefficients in 2, ..., 2, which are holomorphic in a neighbor- 
hood Dy in D bordering on an open part of By. For «w = 0 the expression has value 


0. Hence we can put 
ANU Bas ks y Bg 2) = UB Ss, sss, Baek 


and if herein we replace u by the sum (4), we obtain a representation (3) to the 
following extent. The object (3) is a polynomial in (¢,) with coefficients which are 
holomorphic in Do, and the equality (3) is valid for z « Dy and ¢, = 2,. However, we 
can form the quotient 


P(Z; 2 


V(Z; 2) 


1— 32, —... (eA 
in the entire unit ball D, and it is analytic in the real variables (z,; y,). Therefore 
¥(¢; z) has a holomorphic continuation from Do into D, which completes the proof 
of the theorem. 


* This research was supported by the United States Air Force through the Office of Scientific 
Research of the Air Force Research and Development Command 
J.J. Kohn, “A Boundary Condition for the Vanishing of n Holomorphic Functions in Complex 


n-Space,”’ Proc. Am. Math. Soc., 9, 175-177, 1958. 


LINEAR AND ALGEBRAIC DEPENDENCE OF FUNCTIONS ON COMPACT 
COMPLEX SPACES WITH SINGULARITIES* 
By 8S. BocHNER 
MATHEMATICS DEPARTMENT, PRINCETON UNIVERSITY 
Communicated November 20, 1958 
We will axiomatize the following situation. Take a (non-compact) complex 


space W”, entirely regular, and form the quotient space V W"/R relative to 
some equivalence relation , say relative to the action of a group of homeomor- 


phisms. It is the space V we are interested in. We assume that it is compact, 


but about its “regularity” (or what may be left of it) we only assume as follows. 
For some complex dimension n,n < m, there is in W” a finite number of separate 
complex n-cells, each holomorphically immersed, such that the union of their pro- 
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jections covers V. In the axiomatization we eliminate V altogether, and only 
retain the separate n-cells and a “weak” structure to bind them, as suggested by R. 

In the space of n complex variables z = (z',..., 2") we take a finite number of 
neighborhoods l’,, a 1, each represented by the unit ball 2.) < 1, and, 
for every a, some subneighborhood 7',, as represented by |/2q) < ra, fa < 1. With 
each z, in each U,, we associate permanently some point z, of some 7), (b # a or 


a 


bh = a), insymbols 
$(2,). (1) 


Phat is all. This is our (compact) pseudo-space, and we denote it by S. 
1. A scalar function on S isa set f = (fi, ...,f4) such that f,(z,) is an ordinary 


point function on V, and that equation (1) implies 
Ha(23) =F dlZe)).. (2) 


Our first statement is as follows. Jf each f,(z,) ts holomorphic, then each is a constant 
C, (maximum modulus theorem). 
2. More generally, an automorphic (or density) function arises if for every 


2, € Us, a 1,..., A, there is given a fired number \(z,), with a joint bound, 


\(2,)| <M, and if we replace equation (2) by 
Fala) = (Za) fr(O(Za)). (3) 


The automorphic functions f f f4) constitute, in an obvious manner, a 
vector field over complex constants and the next statement is as follows:! For 


holomorphic automorphic functions this vector space has a finite basis (linear depend- 


ence of automorphic functions). 

3. Even more generally we may assume that the 
not a number but a N by N matrix of numbers, and that each f,(z,) has NV com- 
ponents, so that equation (3) is an affine transformation between such components. 


‘ 


‘automorphy factor’ \(z,) is 


The object (fi, ... , f4) is then representative of a general tensor (linear vector 
bundle over S), and the following statement still remains:? Any vector field of 
holomorphic tensors of a given type has a finite basis. 

!. Finally, we consider meromorphic functions f = (f,). Each f,(z,) is a quo- 
tient Po(2,) :Ga(Z.) and there is a system of factors A(z,) which depends on the given 


function, again with | A(z,)| < .W, such that 


DalZa) (Za) Po(b(Za)), = alZa) N(2a)9n(P(Za)). 


Then the following theorem holds:’ To any n + 1 meromorphic functions f', ..., 


f"*" there is a polynomial P(¢! c"*!) # 0 such that P(f' ft’) =0 (alge- 
braic dependence of meromorphic functions). 

There are versions of this theorem on other types of complex spaces with singu- 
lurities.* 

5). Finally, we note that in all preceding statements the ordinary regular neigh- 
borhoods l’, can be replaced by more general local varieties which are themselves 
already singular. All that is needed is that, for holomorphic functions suitably 
defined on them, a suitable version of the classical Schwarz’s Lemma shall continue 


tohold. Such generalizations are indeed available.® 
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* This work was supported by the United States Air Foree through the Office of Scientific 
Research of the Air Research and Development Command. 

'S. Bochner, “Algebraic and Linear Dependence of Automorphic Functions in Several Vari- 
ables,’ J. Indian Math. Soc., 16, 1-6, 1952. 

2S. Bochner, ‘Tensor Fields with Finite Basis,’ Ann. Math., 53, 402-411, 1951. 

3J. P. Serre, “Fonctions automorphes,” Séminaire Ecole Norm. Sup., 2, (Paris 1953/4); C. L. 
Siegel, ‘‘“Meromorphe Funktionen auf kompakten analytischen Mannigfaltigkeiten,’’ Nachr, 
Akad. Wiss. Gottingen, 2a, 71-77, 1955. 

*R. Remmert, ‘‘Meromorphe Funktionen in kompakten komplexen Riiumen,”’ Math. Ann., 132, 
277-288, 1956. 
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LINE-ELEMENT FIELDS ON THE TORUS* 
By Bruce L. REINHART 
UNIVERSITY OF MICHIGAN 


Communicated by Solomon Lefschetz, November 10, 1958 


1. A line-element field on the two-dimensional torus 7? is a cross-section in the 
unit tangent bundle (or the projective space bundle if we consider unoriented line 
elements). Since either of these bundles is the product bundle 7? & S', we may 
consider such a field as a mapping Ff’: 7? + S', which we shall assume to be of class 
('. We shall parametrize S!' as follows: Consider T? as the quotient of the plane 
k? by the points with integer coordinates. At each point, let (7) be the angle 
between the vector v and the positive x axis. Then 6/27 will be taken as parameter 
on the fiber of the tangent sphere bundle and 6/7 on the fiber of the tangent pro- 
jective space bundle. 

Let Fy:H,(T*) — H,(S') be the induced homology mapping. Then for any 


closed curve C on 7? lying in the homology class C*, F’'4(C*) = uwé, where & is the 
field F’. 


The C' homotopy classes of mappings F are in one-one correspondence with the 


generator of H,(S') and yw is the winding number of C with respect to the 


induced maps F’x.'_ In what follows, we shall be interested in the relation between 
the integral curves of F’, considered as defining a differential equation on 7T?, and 
the map F’x. 

LemMa |. Jf C ts a closed integral curve, the homology class C* lies in the kernel 
of Fy and, in fact, generates it if F's ts not identically zero. 

2. All the closed integral curves of F belong to the same homology class, say 
C*. If * is any other homology class containing a simple closed curve, and T ¢« T'* 
is a simple closed curve, it follows from elementary properties of the intersection 
ring of the torus and the Poinearé-Bendixson theorem? that all integral curves of 
F meet L. 


LEMMA 2. There exists a C' simple closed curve T ¢ T* 


which has everywhere a non- 
zero tangent vector and is tangent to F at only a finite number of points. 

Henceforth, let us denote by Ta curve in ['* having the fewest possible points of 
tangency. Let u(T) be the winding number of [ with respect to F. 

THEOREM. u(T) 0 af and only if Fs 0. In particular, if there is no closed 


integral curve, these conditions hold. 
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Proof: Jf there exists a closed integral curve, then both ['* and C* are in the 
kernel of Fy, so Fy =0. If there is no closed integral curve, let '’ be a cycle with- 
out contact.? Any integral curve of / meets I’ infinitely often; let P' be a point 
of accumulation of these intersections. We may then construct a C! closed curve 
(’ which differs from the integral curve only in an arbitrarily small neighborhood 
of P. C’ and T” always cross in the same direction; hence their intersection number 


is non-zero, so that they neither bound nor are homologous. The winding number 
of T’ is zero because it is without contact, while that of C’ is zero because it differs 
from an integral curve by a small amount. Hence ys = 0 as before. 

Coro.iary 1. If there exists on T? any curve whose winding number with respect 
to F is non-zero, then there exists a closed integral curve of F. 

3. Making use of the cycle T, we may generalize the rotation number A(T), 
which is classically defined for the case that there exists a cycle without contact. * 
Let [ and [, form a basis for H,(T?), and let P4(l) = 7, Px(T) = 7€. 

Proposition. yw(l) = tand X(T) = —j/1, the latter holding only in case u(T) ¥ 0. 

By considering the number of points of tangeney on the cycle T', we may classify 
qualitatively the various kinds of integral curve families of non-oriented line- 
element fields.4 From this classification, the following corollaries are immediate: 

Coro.uuary 2. The non-ortented line element field F is orientable if and only if 7 
and j are both even 

Coro.iary 3. The number of closed integral curves ts at least the greatest common 
divisor of tand j. 

* This research was supported by the Office of Naval Research. 

1N.E. Steenrod, The Topology of Fibre Bundles (Princeton: Princeton University Press, 1951), 


§34. 
2. A. Coddington and N. Levinson, Theory of Ordinary Differential Equations (New York: 
McGraw-Hill Book Co., 1955), chaps. xvi-xvii. 

C. L. Siegel, “Note on Differential Equations on the Torus,” Ann. Math., 46, 423-428, 1945. 


‘H. Kneser, “Reguliire Kurvenseharen auf den Ringfliichen,” Math. Ann., 91, 185-154, 1924. 


ON A PROBLEM OF MAZURKIEWICZ CONCERNING THE BOUNDARY 
OF A COVERING SURFACE 
By PETER SEIBERT 

4 MEMBER OF RIAS, 7212 BELLONA AVE., BALTIMORE, MD. 

Communicated by Solomon Lefschetz, November 13, 1958 
1. In his note! Mazurkiewicz showed that, by introducing a certain metrie on 
a Riemann covering surface and completing the obtained metric space, a concept 
of “boundary elements”? can be defined which corresponds precisely to the usual 
notion of transcendental singularities of meromorphic functions defined on the 
surface. This boundary, considered by itself, is a complete separable metric 
space. Moreover, Mazurkiewiez proved in the above-mentioned note that the 
dimension (in the sense of Menger-Urysohn) of the boundary never exceeds the 
value 2 and raised the question whether covering surfaces with two-dimensional] 
boundaries exist. To our knowledge an answer to this question has, so far, not 


been given. 
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In the present note we give a simple example of a covering surface with a two- 
dimensional (and, moreover, compact) boundary. Furthermore, we show that 
the boundary of a simply connected covering surface is, at most, one-dimensional 
if it is locally compact. 

2. Let F be a two-dimensional manifold, S the Riemann sphere and s an in- 
terior mapping? of F' into S. In this sense we call F a covering surface of S and s 
the projection mapping of F. 

Now consider two points p, g, € F and denote by C the set of all Jordan ares C on 
F connecting p and q. Then, if the diameter of a set in S with respect to the spheri- 
cal metric is denoted by diam, the distance between p and q will be defined by 


p(p, g) = inf diam s(C). (1) 
CeC 


It has been shown in a previous paper’ that this definition yields a metric which is 


compatible with the original topology of F, so that a sequence } p,} tends to a point 


p if and only if p(p,, p) ~ 0.4 The completion of the surface F' will be called the 
extended surface and denoted by F*. Then the boundary Z of F is defined by 


Z=F*—F. 


Therefore, the boundary elements (which we will usually denote by ¢) correspond 
to the equivalence classes of those fundamental sequences of /’ which possess no 
accumulation point on F, 

We now extend the projection s continuously onto the boundary, thus obtaining 
a mapping of the extended surface F’* into S. 

3. In the present section we shall assume F’ to be simply connected. In our paper 
(RS) it was shown that in this case a cyclic order can be defined on the boundary. 
Subsequently, we established a relation between the two topologies on the boundary 
defined by the Mazurkiewiez metric and by the cyclic order, respectively. As a 
corollary of this relation the following result was obtained 

LemMa. The boundary of a simply connected covering surface always permits a 
one-one continuous and order-preserving map into the circular line. (Continuity here 
is meant with respect to the metric on Z.) 

THeoreM 1. [f the boundary Z of a simply connected covering surface F is locally 


compact, the condition 


dimZ S 1 (2) 
holds. 
Proof: Z is a separable metric space.® If, in addition, it is locally compact, it 
can be represented in the form 


: Z uU ®, (3) 


l 
where each of the sets ®, is compact.?. The mapping mentioned in the lemma maps 
each ®, homeomorphically onto a subset of the circle, so that 
dim ®, < 1 (v | ee (4) 


follows. Finally, formulas (3) and (4), together with a theorem of Hurewiez,* yield 
the required inequality (2). 





MATHEMATICS: P. SEIBERT Proc. N. A. S. 


Remark; If the surface F is generated by an entire function f, its boundary is 
totally disconnected and therefore of dimension 0. This statement is true despite 
the fact that some entire functions assume every complex number as an asymptotic 
value.’ Indeed, owing to a well-known theorem of Lindeléf, two non-equivalent 
asymptotic paths of f with finite asymptotic values are always separated by paths 
along which f tends to infinity. In our notation’ this means: Given two bound- 
ary points 1, f € Z with s(f1), s(f2) + ©, each of the sets (1, 2), (€, ¢:) contains 
at least one boundary point over ©. Now let A be a connected subset of Z, for 
which s(A) does not contain the point at infinity. Then ¢(A) is connected, because 
¢ is continuous. Since ¢ is also order-preserving, any two points of g(A) are 
separated by images of points over ©, so that g(A), and therefore A, reduces to a 
single point. The set of boundary points over ©, so that ¢(A), and therefore A, 
reduces to a single point. The set of boundary points over ©, on the other hand, 
is known to be a discontinuum.!' Therefore, we have the following proposition: 
A simply connected Riemann covering surface which does not cover a point a é€S is 
either a hyperbolic type or its boundary ts a discontinuum. 

+. THEOREM 2. There exists a covering surface F of the sphere S whose projection 
mapping s maps the boundary Z of F homeomorphically (and, moreover, isometrically) 
onto S. More generally, given a closed subset A of S, there exists a covering surface 
whose boundary ts mapped by s isometrically onto A. 

Proof: We first consider the case A = S. Denote by d(a, b) the spherical 
distance between two points a,beS. Let }L,{6 bea sequence of sets each of which 
consists of a finite number of segments L,, (a = 1, 2,..., m,) with endpoints 


” . . . ° °,° 
An, yy Gn, , and satisfies the following conditions: 


” ” 


with (n’, un’) ¥ (n”, pw”). 


” 


(i) Ly ye O Lae ye = bo forall n’, nr”, uw’, uw 


yp’ 


ee / " sé ’ . 
(ii) A(ay, yy Qn.) S 2° “wforalln, w. 


(iii) Denoting by V.(B) the eneighborhood }a eS: inf d(a, b) < e! of 
beB 


Bc SS, the relation 


holds for all n. 
Example: Denote by (6, ¢) the point on the sphere of longitude @ and latitude 
<9<2 —< *), and consider the seg I ‘the parallel 
¢(0 S0< 27, —-5 8 ¢ S +5), and consider the segments L,,,,,, of the paralle 


circles ¢ = constant bounded by the points 


” / 
a (Se, 2 eR) aa eS ae re (2 ee), 


where 


jO forn = 0 
(+1 42 ...,+ (2"—1)forn>O, 


Then choose as set L,, the set | L,.,.,} with «, A varying as indicated above. 

We then construct the surface F in the following manner: By removing from 
S the sets of segments!? L,-; and L,, we obtain a sequence of “sheets” F,,(n 
0,1,2,...) which, by joining 7, with F,,,; along both edges of all slits of L,, vields 


9 
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a surface of the required type. (Here the projection s is, of course, the mapping 
which assigns to every point of F the point of S with the same coordinate.) 

The proof will be carried out in two steps by showing that: 

(a) Every point of S possesses an inverse image with respect to sin Z: s(Z) = S. 


(8) Given two elements ¢’, ¢” ¢ Z, the relation 


pte’. &*) d(s(¢’), s(€”)) 


? 


holds. 

Given a point a eS, there exists, according to (ii) and (iii), a sequence | p,} of 
branch points of / such that every p, is an end-point of a segment of L, and s(p,) 
tends to a for n > ©. We have to show that }p,{ is convergent. The points 
P,—-1 and p, both lie on the boundary of the sheet F,.. Now denote by S, the seg- 
ment with end-points s(p,—1), 8(p,) and consider the subset 

D, = \peF,: 8(p) € Va-n,(S,) 
of F,. We conclude from (ii) that none of the cuts L,,, separates p,—1 and p, 
within the set D,, and from (1) that the same is true for any union of such segments. 
Therefore, the two points can be connected by an are in D, and, according to the 
definition of the distance p, the inequality 


P(Pn—1, Pn) S diam s(D,,) = d(s(pp-1), 8(Pn)) + 227" 9 


+: ( { ° . 
holds. Since } s(p,){ converges, so does } p,, and assertion (a) follows. 
To prove (8), consider two points ¢’, ¢” « Z and two sequences of branch points 


‘ y ° Pv a: . ‘ P - 
1p, | and }p,"{ with ¢ and ¢" as their respective limits!’ and assume s(p,’), s(p»”) 
to be end-points of segments L,,,. Then for every n the following pair of inequali- 


ties holds: 


d(s( Pr’ ee Pr”)) = pl Dri Pn ") < d(s( Dr. P s(pr”) ) + Q!-"x, 


The first is obvious. The second is obtained in the same way as (5). Since 


lim d(s(p,’), 8(p.”) d(s(¢’), s(¢")), 

the assertion follows. 

In the case A ¢ S we denote by A,, the set Vs—n,(4A) and omit in L,, all segments 
which are disjoint to A,.'* It is then easy to see that s(Z) is precisely the set A. 

5. Finally, we mention two open questions concerning the dimension of the 
boundary, one of which is topological, the other function theoretic: 

I. Does relation (2) hold for simply connected covering surfaces in general? 

II. Do there exist (concrete) Riemann surfaces with null-boundary (in the sense 
of Nevanlinna) whose boundary (in the sense of Mazurkiewicz) is of dimension 2? 


1S. Mazurkiewiez, ‘Sur les points singuliers d’une fonction analytique,’’ Fund. Math., 17, 26 
29, 1931. 
2 Le., a continuous mapping which transforms open sets into open sets and under which no 
non-degenerate connected set has a single point as its image 
Cf. Seibert, “Uber die Randstukturen von Uberlagerungsflichen,” to appear in Mathematische 
Vachrichten, 19. | This paper will be referred to as “RS” hereafter 
‘For the pertinent arguments in the paper RS, it apparently makes no difference whether the 
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metric introduced there or the one defined by (1) is considered. The distance (1) is called the 
distance with respect to the mapping (s, F) by L. Cesari, who introduces it in compact metric spaces 
(cf. “Surface Area,”’ Ann. Math. Studies, 35, 159 (1956). 

5 Cf. RS, §7(f). 

6 Cf, ibid., §3(c). 

7 Cf. C. Kuratowski, Topologie (Warzawa-Wroctaw, 1950), I, 51, théor. 4. 

8 Cf. ibid., 1, 176, théor. 1; ef. also Hurewiez-Wallman, Dimension Theory (Princeton: Prince- 
ton University-Press, 1948), p. 30, theor. IIT 2. 

* Cf. W. Gross, “Eine ganze Funktion, fiir die jede komplexe Zahl Konvergenzwert ist,’’ Math. 
Ann., 79, 201-208, 1918. 

” Cf. RS, §4 for the relationship between asymptotic paths and boundary points and §5 for the 
notation used here 

1 Cf. ibid. §3(d) 

3 f, , =¢. 

'§ Apparently every boundary point is the limit of a sequence of branch points. Indeed, if 
¢ez, either every neighborhood of ¢ contains parts of infinitely many sheets, or all sufficiently small 
neighborhoods are contained in a single sheet F,. In the first case ¢ is a cluster point of cuts Ln, 
and therefore, according to (ii), also of branch points. The second case cannot occur because 
every boundary point of a sheet F,, is an interior point of F. 

14 We can, of course, assume A to be non-empty. 


THE EFFECTS OF STRONTIUM-90 ON MICk* 
By BarcLay KAmB AND Linus PAULING 
CALIFORNIA INSTITUTE OF TECHNOLOGY, PASADENA, CALIFORNIA 


Communicated November 17, 1958 


On Sept. 19, 1958 there was published in Science a paper by Dr. Miriam P. 
Finkel of Argonne National Laboratory in which she communicated her observa- 
tions on the effects of strontium-90 injected into mice on life expectancy and on 


incidence of tumors of bone and blood-forming tissues.! She discussed the question 
of whether or not the effects are proportional to the amount of injected strontium- 
90 at low doses, and reached the conclusion that it is likely that there is a threshold 
with value for man between 5 and 15 ue. (as compared with the present average 
value from fallout, about 0.0002 ue., and the predicted steady-state value from 
fallout for testing of nuclear weapons at the average rate for the past five years, 
about 0.02 yue.). Her paper ends with the sentence ‘In any case, the present 
contamination with strontium-90 from fallout is so very much lower than any of 
these levels that it is extremely unlikely to induce even one bone tumor or one case 
of leukemia.” 

On the same day, Sept. 19, 1958, newspapers throughout the United States 
published accounts of this work. For example, the Pasadena (Calif.) Star-News 
contained an article with the headline ‘Tests on Mice Show Fallout Safe’? and the 
first sentence, “A woman researcher says tests on mice show that the present fallout 
from nuclear weapons tests will not produce a single case of bone cancer or leukemia 
in humans.” The New York Times published accounts of the work on both Sep- 
tember 19 and September 28. 

We have made an analysis of Dr. Finkel’s data that shows that she had no justi- 
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fication whatever for her concluding statement. All of her data are compatible 
with a zero threshold for strontium-90. Moreover, the statistical analysis shows 
that in order for Dr. Finkel to have been justified with 90 per cent confidence 
(10 per cent type-II error) in making her concluding statement on the basis of her 
data she would have to have used over 1,000,000,000,000 mice in each of her groups, 
instead of the 150 or less that were used. It is hard for us to understand how 
such a serious error could be made in Dr. Finkel’s argument, leading her to 
publish her seriously misleading statement about this matter of great impor- 
tance. 

The Mice Experiments.—-In the studies described by Dr. Finkel young adult 
female mice (strain CF No. 1, about 70 days old) were given a single intravenous 
injection of an isotonic equilibrium mixture of strontium-90 chloride and yttrium-90 
chloride. There were twelve injected groups, ranging in size from 150 mice for the 
group receiving the smallest amount (1.3 we/kg body weight) to 15 for that receiv- 
ing the largest amount (9330 we/kg, an amount that caused death of about 50 per 
cent within 30 days). The control group contained 150 mice. The author states 
that there is 11 per cent retention (at 600 days) of the injected radioactive material. 
Report was made of the fractional decrease in average survival time, the incidence 
of animals with osteogenie sarcomas (among 150-day survivors), and the fractional 
decrease in time to a 20 per cent incidence of reticular tissue tumors compared with 
the 20 per cent incidence time of the controls. 

Studies of this sort may be of great value in providing information about the 
probable amount of damage done to human beings by exposure to high-energy 
radiation, such as that from strontium-90 produced by nuclear weapons. It is 
important that the analysis of the experimental results be carried out correctly. 
We have found that in the treatment of problems of this sort the assumption that 
the probability of damage is strictly proportional to the amount of radiation ex- 
posure does not in general require that a response such as decrease in life expectancy 
be linear, except over a very small range. Moreover. we have found that this 


assumption together with the theory of statistics can be applied in a reasonably 


straightforward way in the discussion of data such as those obtained by Dr. Finkel, 
as shown in the following sections. 

Analysis of the Experimental Data on Life Shortening.—Our analysis proceeds 
from the hypothesis, induced by Lewis? as a result of his study of the incidence of 
leukemia, that exposure of the bone marrow of an animal to radiation results in an 
increase in the probability per unit time that the animal wil! die at any time there- 
after, the increase being proportional to the quantity of radiation absorbed. We 
shall suppose that this hypothesis applies to all of the radiation-induced effects in 
Dr. Finkel’s experiments with mice. 

Let No be the number of animals at the beginning of a given experiment, 
taken to be at ¢ 0, and let N(t) be the expected number (average for many 
experiments of the same kind) at the later time ¢. Further, let N°(t) be the 
expected number in a “control” experiment in which no strontium-90 is injected, so 
that 


1 dN°® 
No dt 
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is the natural specific death-rate function. We denote by a the quantity of stron- 
tium-90, in we/kg body weight, that is retained in the animals. Then our hypothesis 
vields the equation 


dN ' N ; i 
di —NBat — Noy NI (1) 
where 8 is a constant of proportionality relating the quantity of strontium-90 re- 
tained to the increased probability per unit time that the animals will die, this 
probability of course increasing linearly with time owing to the nearly constant 
irradiation by the decaying strontium-90. 
Let n°(t) N°(t)/No, so that g(t) —n'(t). Then on integrating equation (1) 


, : | t ‘n%(t) 
N No exp: | = Bat? — | dt 
2 » wt) 


ro 2) Bal? ’ 2) Bat? . 
Nyn(te—" a” = N%(t)e~" i (2) 


we obtain 


To compare this result with the experimental data we calculate A, the fractional 
decrease in life expectancy (fractional decrease in average survival time after in- 


jection), 


where (, is the life expectancy for a retained quantity @ of strontium-90, 


| — 
be No | N(t)dt 


The equation for A is 


(4a) 


(4b) 


where for simplicity we put 7 , 23. The result in equation (45) is obtained by 


wn integration by parts, and the error function is defined as 
2 ee : 
Erf x | e” dy 
VT Jo 
If normally (for a = 0) all animals lived to the age f and then died, so that g(t) 
were a delta function 6(¢ — to), then we would have simply 
Var Ert lyoW¥ 


A ; 
2 bV y 
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However. the actual lifetimes scatter with sizable dispersion about t. The 
extent of this dispersion can be estimated from the acceptance region A < 0.07 
quoted by Dr. Finkel (Figs. 3 and 4) as appropriate to a test of the hypothesis of 
no difference between the responses of the control population and of a population 
injected with a given dose of strontium-90. If we assume (1) that the test of no 
difference applies to Curve A, the life-shortening data, (2) that the test was one- 
tailed, and (3) that the test accounted for the uncertainty in the mean lifetime of 
the control population and for the uncertainty in the mean lifetime of her group 10, 
the highest-dosed population to fall within the acceptance region (except for the 
“peculiar result’? for group 8), then we find that the estimated standard deviation 
for g(t) is ¢ = 258 days. 

These assumptions are somewhat uncertain, as explained later, but they are the 
best that can be made from the information given in Dr. Finkel’s paper. The 
uncertainty in drawing any conclusions about g(¢) from Dr. Finkel’s data lead us to 
take a more general approach. Gompertz discovered that for animal populations 
the logarithm of the ‘‘age-specific death rate”’ is closely a linearly increasing fune- 
tion of time. For man the age-specific death-rate doubling time is about 8 years. 
Jones*® has pointed out that the doubling times for different animal species are ap- 
proximately proportional to the mean life spans for the species. We shall use this 
information to derive a hypothetical death-rate function g(t) for the mouse popu- 
lation used in Dr. Finkel’s experiments. 

The Gompertz Jaw is 


In g(t) C + Bt 
which yields 
0 . Bi | = 
n(t) exp [—A(e™ — 1)] (7) 
where A (= e°/B) is a constant and where B is related to the doubling time rp by 
In 2 


TD 


(S) 


A is to be chosen so as to give the correct mean life span: 


. 1 
exp [—A(e® — 1) dt = 5 [-Bi(—A)] 


The exponential integral Ei(v) is defined by Jahnke and Emde.‘ If ¢ Q is 
tuken to be a time shortly after birth, but long enough after birth to exelude infant 
deaths (which are omitted in Gompertz’ treatment), then f is 7’, the mean life 
span from birth to death. If then 7/rp is a constant for all animal species, we 
find from equations (8) and (9) that A is a constant, independent of species, given 


by the solution of the equation 


(10) 


Assuming 7’ 60 vears for man, with 7p 8 yveurs, we find A 0.0032. The 


solution of equation (9) is obtained with the help of the expansion! 
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— Ei(—-2x) —InT'z + 2 — — ee (11) 


where I 1.781. 
The death-rate function is 


g(t) = AB exp [Bt — A(e® — 1)] 


The dispersion of life spans is measured by the standard deviation o of g(t): 


a” [ tq(t)dt — to? = 2f in(t)dt — ty” 


‘The second expression results from an integration by parts. A numerical integra- 
tion is required to obtain o*, which can most easily be carried out with n°(t) values 
from equation (7). In this way we find from equations (9) and (10), with t) = 7, 
that 


a 0.247 


For Dr. Finkel’s mice, reported to have T = 670 days, we have ¢ = 161 days, in 
rather poor agreement with the value o = 258 days inferred above from her paper. 

In the calculations that follow we have used equations (7) and (10), with the as- 
sumption 7’ = ft) = 600 days, although actually the mice were about 70 days old 
at the beginning of the experiment. Thus we have used a doubling time rp of 80 
days, and our g(t) has standard deviation ¢ = 141 days. The assumption 7p = 
80 days agrees with the value quoted by Jones’ for mice. Values of n°(¢) for these 
parameters are given in Table 1. The difference between assuming 7 = 600 days 
and assuming 7’ 670 days is not great; in fact, survival curves calculated from 
equation (5), which assumes ¢ = 0, do not differ greatly from curves obtained by 
the more refined procedure that we have used. 


TABLE |! 


t 
Days) u(t) 
0 OOO 
SO Q97 
160 990 
240 978 
320 953 
100 908 
ISO SIS 
560 666 
640 145 
720 195 
SOU 038 
SSO OOO 


We proceed now to compare equation (4a), evaluated with the help of equations 
(7) and (10), with the experimental life-shortening data. We assume with Dr. 
Finkel that a 0.1 la* where a* is the injected dose of strontium-90 in we/kg body 
weight, and we attempt to choose the available parameters so as best to reproduce 
the observed life-shortening data A(a*). There are two parameters: the constant 
8, and the no-dose life-shortening Ao, the latter arising from the fact that we cannot 
give great weight to Dr. Finkel’s zero point because of the statistical uncertainty 





Von. 45, 1959 PATHOLOGY: KAMB.AND PAULING 59 


in the observed mean life span of the control population. Thus the theoretical 
curve to be fitted to the data is 


“ 
A(a*) = Ao +1 —- 4 Re re a (12) 
0 
The inclusion of Ap simply as a constant in equation (12) is not strictly correct: the 
additive term arising from an adjustment of the zero point should be written 
Ao(a*), where Ao(a*) is a somewhat complicated decreasing function of a* that 
tends to zero as a* — ©. For simplicity, however, we ignore this complication, 
which proves to be unimportant for the ower radiation levels (a* >= 1000 ye/kg), 
and which in any case does not much change the results obtained, because Ao is 
small. 
For very small values of a* equation (12) reduces to 


| eo 
A(a*) Ao + a* E (O.11)¢ | nnneat| 


dA 
Ay + a*® (13) 


Q@ * 
da* * 


a 0 


This is the linear response region. We can therefore choose preliminary values of 
Ay and 6 by estimating a linear fit to the experimental points at low values of a*. 
The integration in equation (13) is performed numerically, with use of equation (7). 

We have calculated theoretical curves from equation (12) in three steps: (1) For 
a* < 50 ue/kg equation (13) applies; (2) for selected values of a* in the range 50 < 
a* < 1000 we carry out the integration in equation (12) numerically, using time 
intervals At = 80 days; (3) for a* > 1000 it is found that the asymptotie form of 


equation (12) is valid: 


I T 
A(a*) Ao +1 —- (14) 
7 ty Vo.11g0* 


Because we wish to examine the result statistically, we adjust the parameters by 
a weighted least squares procedure. We calculate two theoretical curves y fo + 
f(x, By) and y fo + f(x, Bz) (here y is fractional life shortening, A; x is injected 
dose, a*; and fy is the constant Ao) for two nearly correct values 8; and 8» of the 
parameter 8. We ask for values f fo + Afo and 8 B, + (B. — B,)6 of the 
parameters such that the weighted sum of the squares of the differences between 
the experimental values y,; and the theoretical values y(a,) is a minimum: 


dwily; — fri, B) — fe)? min (15) 
t 


Since B. — 8, is small we can assume that 
fi bt, B) fi Lt, B,) + Lf Li, (2) ~ ie fl Vi; B,) |6 
f Bis By) io Af 6 


The parameter adjustments Afy and 6 are then given by 
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(dow Ay,) | S w,Af;?) — ( yw Ay; Af;) (> w,Af;) (17) 
] 
(ow, (> w,Af;?) — (> wiAf;)? 


(dow, Ay, Af;) (dow, — (>ow,Ay,) (dow, Af;) 
(dw, (wf?) — (Cw,df)? 


where Ay, Yi — fo — f(x, Bi). 


The weights w, appearing in equations (15) and (17) should be inversely propor- 
tional to the a priori variances of the experimental values y;. We take the vari- 
ances to be inversely proportional to the number of animals in each experimental 
group. This ignores the effect of radiation in changing the dispersion of life spans, 
but a detailed examination shows it to be a not unreasonable procedure. 

The theoretical curve obtained in the above manner is shown in Figure 1, with 


Fig. 1.—Percentage decrease in life expectancy, A, as a fune- 
tion of injected dose a* of strontium-90. Solid curve is the theore- 
tical curve calculated from equation (12). Solid circles are the 
experimental values reported by Dr. Finkel. 


the experimental points for comparison. The two points at highest radiation levels 


lie well above the curve, doubtless because the mechanism of life shortening at the 
high radiation levels departs from what we have assumed, owing to the importance 
of subacute and acute irradiation disease, which Dr. Finkel reports to be the pri- 
mary cause of death at injected doses above 2200 ue/kg. These points make little 
contribution to the least squares parameter adjustment, owing to their low weights, 
and can be omitted without sensibly changing the result. The parameters ob- 
tained are 8 = 1.8 X 1077 day? (ue/kg retained) ~!, and Ay = 2.5 per cent. 

The least-squares-fitted curve can be used to estimate the death-rate standard 
deviation o: for experimental A values of unit weight (taken here to be for the con- 
trol group and “group 12”), the estimated variance of the experimental A values is 

1 ; 
> wily: — fo — f(x; B)/? (18) 


_) — 


m Fear 


where m is the number of experimental points. Equation (18) takes into account 
the two-parameter adjustment. If J/ is the number of animals in groups having 


unit weight, then 
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oe? = MG,? (19) 


From equations (18) and (19) we find ¢ = 191 days in case the two highest points 
mentioned above are omitted (@ = 222 days in case they are included). Com- 
paring the value 191 days with the value « = 161 days based on the Gompertz 
death-rate curve (670-day life span) and the value ¢ = 258 days inferred from Dr. 


Finkel’s data, we see that the theoretical curve fits the experimental values about 
as Well as would be expected from the Gompertz curve, and somewhat better than 
would have been expected on the basis of Dr. Finkel’s acceptance region. 

In Figure 2 there is shown the portion of the theoretical curve for the lower radi- 


€ 


Fic. 2.—Percentage decrease in life expectancy in the low-dose 
region of Dr, Finkel’s experiments. Theoretical curve and experi- 
mental points as in Fig. 1. 
ation levels. The curvature is pronounced, and the linear response region is re- 
stricted to injected doses less than about 50 we/kg. Most of Dr. Finkel’s experi- 
ments were carried out in the nonlinear portion of the curve. 

It is interesting to compare the above analysis with an alternative one based on 
the approach developed by Jones,* in which the effect of a given exposure of an 
animal to radiation is regarded as equivalent to an increase in physiological age of 
the animal by an amount proportional to the amount of radiation received. In 
terms of the Gompertz formulation of the natural death rate, this results in the 
case we are considering to the addition of a linear term in ¢ to equation (6): 


1 dN ' - 
In| — N dl = (‘+ Bt+ Brat (20) 


The constant By in this treatment plays a role analogous to the constant 8 used 
previously. 


From equations (3) and (20) we obtain 


Ki ( - 
1 + na A 


e 
Bl + na) 
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Cc » ° ° ° ° 
where A ¢ /B as betore. In practice A is so small that the exponential integral 


can adequately be approximated by the logarithmic term in equation (11): 
7 A | 
—Ki { — In + In (1 + na) 
l +90 PA 


Recognizing from equations (21) and (22) that 


and making use of equation (8), we obtain 


1 + «In (1 + na) 
1 + na 


A 


where « tp/ (ty In 2) = 1/5. 


By choosing 7 0.014 (uwe/kg retained)~', and by adjusting the zero point 
slightly as done previously, we calculate from equation (23) a theoretical curve 
that matches closely the curve calculated from equation (12), which is the curve 
shown in Figures | and 2. The discrepancy in A between the two curves is less 
than 0.02 over the range 0 < a* < 3000 we/kg, and increases to 0.04 at a* 
9000 uwe/kg. The two curves fit the experimental data equally well, as shown by the 
estimates ¢ 191 days for the curve from equation (12) and ¢ 189 days for 
equation (23), calculated by the weighted sum-of-squares procedure described pre- 
viously. Life-shortening data, at least of the accuracy involved here, are therefore 
unable to discriminate between the two analytical approaches. 

Analysis of Incidence of Leukemia and Related Diseases.—The experimental data 
on the incidence of diseases of the blood and blood-forming tissues can be analyzed 
in the framework of the above treatment. However, because of the peculiar form 
in which the experimental results are presented (“Curve C: percentage decrease in 
time to a 20 per cent incidence of reticular tissue tumors compared with the 20 per 
cent incidence time of the controls’’), the analysis is subject to greater uncertainties 
and difficulties and the data cannot so readily be evaluated statistically as those 
for the decreased life expectancy. We therefore content ourselves with a somewhat 
sketchy treatment, which should suffice to indicate the general nature of the 
problem. 

Let A(t, a) be the expected number of deaths due to these diseases that have 
occurred by the time fin a population having retained body burden a of strontium- 
90. We may then expect to find a death-rate probability parameter 8, for these 
diseases such that the death rate is 

dr ; dX N(t) 
N(t)Biat + ; 
dt dt Nt) 
Where A(f) is the number of deaths due to these diseases expected in the control 
population. To carry the analysis further we need to know the function do(t), but 
unfortunately Dr. Finkel presents no data that enable us to determine it. Of the 
various assumptions that could be made, we have chosen to assume that the nat- 
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ural deaths due to leukemia are distributed as though they were radiation-induced 
according to the same model as the deaths due to radiation from strontium-90. 
The natural leukemia death rate will then be equivalent to a “background” body 
burden a of strontium-90, and equation (20) becomes 

dx 


= N(t)B,(a0 + a)l 
dt 


Obtaining V(¢) from equation (2), we have 


X(t) re 1/2) Bat 
\ = B,(ao + a) | n° (bye ee had 


The expected 20 per cent incidence time 7 is then the implicit solution of 


0.118,(ao* + a*) | n(t)e~ OM/2BaNt® + dt 


a | 


and the expected no-dose 20 per cent incidence time 7» is given by 


0.118 ,a* | n'(t)t dt = (25) 
7 as given by equation (25) is not necessarily the same as the 20 per cent incidence 
time 7)’ = 465 days observed for the control population. 

To compare the theory with the experimental data we calculate from equation 
(24) the fractional decrease function | — r(a*)/7’. An adequate approximate 
calculation for values of 7 less than about 450 days (1 — 7/70’ > 0.20) can be made 
by approximating °(f) by a Gaussian e-“, as shown in Table 1. In this case 


equation (24) becomes 


a (0.11)Ba* 


0.118,(ao* + a*) | 
i = exp (=r( t 30.118") 


To evaluate the parameters 8, and ay* we have fitted a smooth curve, by eye, to the 


experimental values of 1 — 7/70’, and used this curve to pick pairs of values (a*, 
t(a*)) from which the quantity 0.118,(ao* + a@*) was calculated from equation 
(26). The quantity 0.11 B,ao* was caleulated from equation (25) by numerical 
integration, with the assumption 7» rt). When plotted against a*, the values 
of 0.11 B(ao* + a*) calculated in this way lie nicely along a straight line, as re- 
quired by the theory, for values of a* in the range 0 < a* < 1000 we/kg. Above 
1000 we/kg the linear relation breaks down, reflecting the fact that the one experi- 
mental value in this higher range, at 2200 ue/kg, lies rather far from the theoretical 
curve. Ignoring this highest value we obtain in this way the parameters 8, 

0.7 X 1077 day~? (we/kg)~! and ay* = 200 we/kg, from which the theoretical curve 
shown in Figure 3 is calculated. In addition to the point 7(0), and the points 
t(a*) calculated from equation (24) over the range of validity of the Gaussian ap- 
proximation, we have calculated the slope of the theoretical curve at a* = 0 from 


the following formula, which can be derived from equation (24): 
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ka 
3 


Fic. 3.—Percentage decrease in time to a 20 per cent incidence 
of bone tumors, as a function of injected dose of strontium-90. Solid 
curve is calculated from equation (26). Solid circles are experi- 
mental values reported by Dr. Finkel. 


| n(t)t°dt 
Bo So 


: | " 
L 108, sss 2 
| n(tjidt 
d ( 4 Jo 


da* 


"9 (Oo04 )B, 


T 1O-yo7797"( 70) 


where yo = '/2(0.11)8,a0*. The ratio of integrals appearing in the second term of 
the numerator in this equation can be shown to have a value close to unity (actually 
1.06). 

A comparison of the parameters 8, = 0.7 X 10-7 and 8 = 1.8 X 1077 suggests 
that of the radiation-induced deaths the fraction due to leukemia and related dis- 
eases in Dr. Finkel’s experiments on mice is rather larger than has been estimated 
for man. <A particular sensitivity to these diseases on the part of this strain of 
mice is suggested also by the large “background dose level’’ ao*, reflecting the rela- 
tively large number of deaths due to these diseases in the control population. 

Statistical Examination of Dr. Finkel’s Conclusions.—In searching for evidence 
for the existence of a threshold body burden of strontium-90, below which no harm- 
ful effects are caused, Dr. Finkel uses two methods: (1) statistical analysis of the 
experimental data, and (2) extrapolation of experimental curves. We now con- 
sider these two methods. 

The statistical analysis consists of a ¢-test of the hypothesis of no difference in 
response between the control population and a population dosed with strontium-90. 
Dr. Finkel accepts the null hypothesis at the 10 per cent significance level (“10 
probability level or higher’) for the three lowest-dosed experimental groups, and 
considers that this acceptance constitutes “evidence that there might be a thresh- 
old” or that “a threshold ... may lie between 4.5 and 44 we/kg.”’ 

It constitutes nothing of the kind. It is clear that the width of the acceptance 
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region for the null hypothesis (shaded region in Figs. 3 and 4 of Dr. Finkel’s paper) 
should vary inversely as the square root of the number of animals in the experimental 


groups, assuming approximate normality of the death-rate curve g(t), as Dr. Finkel 


must have done in applying the ¢-test. The threshold for which she finds “‘evidence”’ 
in the experiments is thus no threshold at all but simply a reflection of the statistical 
uncertainty of her information. It is clear that she could have found ‘‘evidence” 
of this sort for a threshold at any arbitrarily large radiation level (perhaps short 
of what would produce acute radiation sickness) by simply using few enough 
animals in her experiments. 

The fallacy in Dr. Finkel’s statistical argument is a failure to control the prob- 
ability of type-II error of her test. T'ype-LI error’ is acceptance of the null hypothe- 
sis When it is in fact false. Consideration of the type-II error requires consideration 
of the alternative to the null hypothesis, which in this case is the theoretically likely 
linear response at low doses. If we use for the slope dA/da*),«—) of the life-short- 
ening response at low doses the value obtained above (eq. [13]) from a study of 
Dr. Finkel’s results, namely, dA/da* 0.14% (ue/kg) 
Dr. Finkel’s é-test acceptance region is appropriate to a one-tailed test at a* = 8.9 


and if we assume that 


uc/kg, the highest experimental value for which the null hypothesis was accepted, 
then we can calculate the probability of type-II error. It is 85 per cent. This 
means that if there exists in fact no threshold at 8.9 we/kg, Dr. Finkel’s test would 
nevertheless have produced ‘evidence’ for one in 85 experiments out of every 100 
experiments performed. On the other hand, if there were in fact a threshold, the 
test would deny it in only 10 per cent of the experiments. Evidently the test is 
worthless as a proof of the existence of a threshold at this dose level (or lower, for 
which the probability of type-I] error approaches the maximum that is possible, 
90 per cent. fora 10 per cent probability of type-I error). 

It is incumbent upon those who would extrapolate their threshold conclusions 
from 150 mice to 3 X 10% human beings that they demonstrate the existence of a 
significant experimental departure from the theoretically likely linear response, 
because although the existing burdens of strontium-90 are low, the number of 
individuals involved is very large, and the harmful consequences of proceeding 
on an unfounded assumption of a threshold are great.6 As we have shown above, 
Dr. Finkel’s results are in complete harmony with a linear law; in fact, the agree- 
ment between the linear law and the experimental results is better than could have 
been expected on the basis of the width of her null-hypothesis acceptance region. 

As an alternative to the statistical tests, Dr. Finkel determines a threshold by 
extrapolating the experimental life-shortening curve. She states: “‘Since the [life- 
shortening | values for 1.3, 4.5, and 8.9 we/kg do lie along a straight line when plotted 
semilogarithmically, it may be argued that they represent true departures from 
the control value. An extension of this straight line crosses the control value at 
0.4 we/kg.”’ It is difficult to see why the semilog plot rather than some other should 
be used for the extrapolation. But in facet an extrapolation of any kind is ground- 
less. The three response values lie within less than half the range of probable 
error (within —'/: P.E. to +'/. P.E.) of the difference d between experimental 
values of A, as determined from the width of the null-hypothesis acceptance region 
(og = 5.5 per cent). If it is not obvious that no non-zero regression slope deter- 


mined from these points can have any statistical significance, one can show? that 
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the standard deviation of the regression slope estimator derived from the semi- 
logarithmic plot is 2.2 times the estimated slope itself. If the semilogarithmically 
linear relation of the three points can be ascribed to anything but chance, then all 
of Dr. Finkel’s statistical arguments are false. It is hard to imagine how two such 
mutually contradictory “proofs” could be advanced at one time. 

There are other statistical points in Dr. Finkel’s paper that merit scrutiny. In 
our discussion of her results we have had to rely on the correctness of the null- 
hypothesis acceptance region that she presents, but there are serious reasons for 
doubting its correctness. The width of the acceptance region corresponds to the 
estimate ¢ = 284 days for the standard deviation of the death-rate function g(t), 
if it applies to a one-tailed test on the difference between the life-shortening values 
obtained from two experimental groups of 150 animals each. On the other hand, 
Dr. Finkel’s statement’ that groups of 1393 animals would have been required to 
establish as significant at the 1 per cent level the difference observed (2.5 per cent) 
at the lowest dose corresponds to ¢ = 171 days, a gross discrepancy. The latter 
value, we note, agrees reasonably with the values ¢ = 161 days from the Gompertz 
relation or ¢ = 191 days from the agreement between experimental data and our 
theoretical curve. 

It is clear that since the number of animals differs from one experimental group 
to another in Dr. Finkel’s experiments, the null-hypothesis acceptance region can- 
not have width independent of injected dose a*, as shown in her figures. From the 
information given there is no way to tell to which experimental groups the test 
appropriately applies. 

More serious is the evident fact that Dr. Finkel applies the same acceptance re- 
gion indiscriminately to the three very different sets of experimental data repre- 
sented by her curves A, B, and C. It seems likely that the test was designed to 
handle the life-shortening data (curve A), because a ¢-test would not be inappro- 
priate to life-span data, since the death-rate function g(t) is (rather crudely) Gaus- 
sian. A statistical analysis of the curve-C data would be difficult, because the 


experimental statistic 7 (20 per cent incidence time) is cumbersome to handle mathe- 


matically, as is evident in our discussion. But it is easy to show that Dr. Finkel’s 
acceptance region is entirely inapplicable to the curve-B data (“proportion of ani- 
mals that survived the latent period of 150 days and then died with osteogenic 
sarcomas’’). 

The number of bone-cancer deaths in populations of a given size during a given 
time interval will be Poisson-distributed, if we neglect variations in population size 
due to deaths during the first 150 days, which is legitimate, as can be seen from 
Table 1 or from numbers given by Dr. Finkel, which show that the control group 
still contained close to 150 animals at ¢ = 150 days. Whatever the low-dose re- 
gression function for curve B, it is clear from Dr. Finkel’s Figure 4 that the ex- 
pected number & of bone-cancer deaths is close to 3 for groups of 150 animals not 
dosed with strontium-90. To find acceptance regions for the null hypothesis of no 
significant difference in the number of such deaths between the control population 
and a dosed population of equal size we therefore find the value of the difference 6, 


such that the probability of type-I error (one-tailed test) is P: 
6p + « en 


1-—P - ¢ 
; 2 be aus a 


(27) 
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The results of a numerical evaluation of equation (27), in case = 3, are P = 0.10, 
bp = 2.6; P = 0.01, 6p = 5; P = 0.001, 6p = 8.5. Using the fact that the observed 
number of control deaths was 3 (2 per cent of 150), we find that the upper limit of 
the acceptance region at 10 per cent significance level is 3.7 per cent of 150. From 
Dr. Finkel’s Figure 4 we therefore see that the highest strontium-90 dose that pro- 
duced a “statistically non-significant increase’ in the number of bone-cancer 
deaths is 8.9 uwe/kg, not 200 we/kg as stated by her. Her acceptance region repre- 
sents for curve B a test having 0.1 per cent probability of type-I error. Her entire 
discussion of the statistical significance of the curve-B data is erroneous. 

The Proper Testing of Evidence for a Threshold.—¥From the above discussion it is 
clear that a valid statistical test of the null hypothesis of “no response” at a given 
radiation level or a given dose of strontium-90 must use the type-II error as the 
basic parameter, rather than the type-I error, as is employed in the standard “cook 
book”’ tests, which are designed basically for application to the manufacture of 
goods. Alternatively stated, the null hypothesis that must be tested in the stand- 
ard way is the hypothesis that the observed response values are in accord with a 
linear response curve at low doses. 

Our analysis of Dr. Finkel’s data for mice enables us to estimate reliably the 
linear decrease in life expectancy for low doses of strontium-90, and it is therefore 
possible for us to determine how many animals would have to be used in an experi- 
ment in which the mean lifetime of a control group is compared with the mean 
lifetime of a dosed group in order to establish the existence of a threshold at or 
above the dose used. We may consider two types of test: (A) the “minimal” test, 
that is, the test that requires as few animals as possible; (B) the “most powerful” 
test, Which minimizes the probability both of type-I and of type-LI error. 

The null-hypothesis “‘no response”’ for test A is to be accepted if the dosed group 
exhibits no decrease in life expectancy, or an actual increase, When compared to the 
control group. Clearly this acceptance region makes the test minimal, because the 
probability of type-I error is 50 per cent, so that the test gives a neutral decision in 
case a threshold actually exists. If the number of animals used is greater than re- 
quired for test A, then a decision as to the existence of a threshold will be more 
often right than wrong, in case that the threshold does actually exist. At the same 
time we can protect ourselves adequately against the serious alternative possibility 
by suitably choosing the type-II error. 

Since the expected decrease in life expectancy for the dosed group is a-dA/da \o, 
the type-I and type-II errors are simultaneously minimized, and made equal, by 
choosing as the upper limit of the acceptance region for test B a decrease in life ex- 
pectancy of '/sa-dA/da |, 0.63° >:a, Where @ is given in uwe/kg retained in the 
body. 

Since the expected decrease is proportional to a, the number of animals M re- 
quired for the control group, if an equal number is used for the dosed group, is given 
by 

‘s 


M (28) 


a 


where v is a constant that depends on the type of test (A or B), on the probability 


of type-II error, and on the standard deviation o of the natural death-rate function 
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g(t). Because M proves in all cases of interest to be large, it is adequate to use the 
normal distribution in computing the constants v in equation (28), owing to the 
Central Limit Theorem. 


TABLE 2 
VALUES OF THE CONSTANT vy IN EQUATION (28) 
Type or Test 
A 3 
(‘Minimal’) (‘‘Most Powerful’) 
o (days) 284 170 284 170 
10 per cent probability of type-II error 1560 1630 18240 5520 
1 per cent probability of type-II error 14950 5360 59800 21440 

Using all of these principles, we have computed the coefficients v for the various 
circumstances shown in Table 2. In particular we compare the results for ¢ = 
284 days, derived from Dr. Finkel’s acceptance region, with the results for ¢ = 170 
days, which seems most reasonable on the basis of the previous discussion. 

The numbers vy given in Table 2 are equal to the number of animals in the control 
and in the dosed groups required to establish the existence of a threshold at a* = 
9.1 pe/kg, just above the highest injected dose (8.9 we/kg) for which Dr. Finkel 
accepted the hypothesis that a threshold exists. The numbers of animals used in 
her experiments were too small by factors of 10 to 400, for the conclusion that she 
reached. By solving equation (28) for a we may compute very simply the lowest 
threshold a7* that could have been recognized with statistical significance in her 


experiments, assuming that 150 animals were used both in the control group and in 


the dosed groups, which in general was not the case (fewer were used). These val- 
ues of ar* are 91 and 54 for test A and 181 and 109 for test B (in each ease for 
ao = 284 days and 170 days, respectively). It is clear from the experimental data 
that no threshold exists at any of these levels, and accordingly we are required to 
conclude that Dr. Finkel’s data show that there is no threshold large enough to 
have been recognized with statistical significance from her data. 

Conclusions about Effects on Man of Strontium-90 from Fallout.-We now turn to 
the discussion of Dr. Finkel’s conclusion that the present contamination with stron- 
tium-90 from fallout is so very much lower than the “threshold” levels that it is 
extremely unlikely to induce even one bone tumor or one case of leukemia. 

This statement by Dr. Finkel is shown by the argument given above to have no 
justification whatever from her experimental results, obtained with 150 mice or 
fewer in her control group and injected groups. We may ask how many mice 
would be needed in each group in order to permit Dr. Finkel’s statement to be 
made with statistical significance (or to be shown to be false). 

The present average body burden of strontium-90 in the world’s population is 
about 0.0002 ue. per person. This corresponds, with Dr. Finkel’s conversion factor 
(5 to 10 we. per 70-kg man equivalent to | we. retained per kg for mice) to a retained 
dose a = 0.00002 to 0.00004 we/kg in mice. Hence in order to justify Dr. Finkel’s 
statement evidence would be needed that the mouse threshold is as great as about 
0.00004 pe/kg; that is, we must place @ in equation (28) equal to 0.00004 uwe/kg. 
From the values of the constant v in Table 2 (we use the values for o = 170 days, 
which we believe to be better than those for ¢ = 284 days) we find MW = 1 Xx 10” 
for the “minimal” test and 3.4 X 10! for the “most powerful” test with 10 per cent 
type-II error, and 2.3 X 10’ and 13.5 & 10', respectively, with | per cent type-II] 
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error. We hence conclude that a study like that made by Dr. Finkel would have 
to use a much greater number of mice than the number of people in the world, in 
order to provide evidence that would justify her extreme statement to be made 
with statistical significance. 

This conclusion is, of course, not at all unexpected. The difficulty of detecting 
by statistical methods an effect that causes a small increase in the annual number 
of deaths among the world’s population is well known. For example, let us as- 
sume that the average number of deaths per year is 50 million. The statistical 
fluctuations from this average from year to year are measured roughly by the 


square root of this number, 7000; and accordingly a study of a larger population 


would be needed to show with statistical significance the existence of an effect 
resulting in an additional 1000 deaths per year (the rough estimate of the world- 
wide effect of the present body burden of strontium-90 from fallout, if there is no 
threshold). The same number of mice would be needed to test the equivalent 
effect in mice. 

Summary.—We have developed methods of theoretical analysis of the results of 
experimental studies of the eifects of injection of radioactive substances into ani- 
ma!s on their life expectancy and on the incidence of tumors. These methods have 
been applied to the data reported for mice by Dr. Miriam P. Finkel, and it has been 
shown that her conclusion from these data that it is extremely unlikely that the 
strontium-90 from the fallout from nuclear weapons tests will induce even one bone 
tumor or one case of leukemia in human beings is completely unjustified. 

* This paper is a contribution from the Division of the Geological Sciences (No. 908) and the 
Division of Chemistry and Chemical Engineering (No. 2421) of the Institute. A brief account of 
the work has been published by us (Letter to the Editor, The New York Times, Nov. 16, 1958), 
and a reply has been made by Dr. Finkel (cbid., Nov., 30 1958). 
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GRAVITATIONAL PROPERTIES OF ANTIMATTER* 
By L. I. ScHIFF 


DEPARTMENT OF PHYSICS AND INSTITUTE OF THEORETICAL PHYSICS, 
STANFORD UNIVERSITY, STANFORD, CALIFORNIA 


Read before The Academy, November 7, 1958 
Kétvés Experiments and Equivalence Principle... A very precise series of experi- 
ments involving gravity was performed: by E6tvés and collaborators between 1890 
and 1922.' These experiments sought for possible variations in the ratio of gravi- 
tational to inertial mass from one substance to another. By “gravitational mass”’ 
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we mean here what Bondi? refers to as the ‘passive gravitational mass,”’ that is, a 
quantity proportional to the force exerted on a body by an external gravitational 
field, say of the earth or the sun. In contrast, one could also work with the active 
gravitational mass, which is proportional to the force exerted by a gravitating body; 
we shall return to a discussion of this quantity in the last section. By ‘inertial 
mass’? we mean the ratio of some applied mass-independent force (electromagnetic, 
for example) to the acceleration of a body. We choose units such that th ratio 
measured by Eétv6s is equal to 1 for normal matter. 

The experiments were carried out for a number of substances, with varying de- 
grees of accuracy. Three of these substances are especially convenient for our 
purpose, since they are nearly elemental and spread throughout the periodic table, 
and they are among those for which the ratio was measured with greatest precision: 
magnalium (90 per cent aluminum, 10 per cent magnesium), copper, and platinum. 
Comparison of magnalium and platinum showed that the difference in the ratio of 
gravitational to inertial mass was (0.4 + 0.1) X 10~-%; this ratio difference was (0.4 
+ ().2) & 10-8 for copper and platinum. The authors did not regard these small 
positive effects as significant, and we shall follow this interpretation of their work. 
We can therefore say that for these three substances there is no difference in the 
ratio to within about '/. part in 108. 

The assumption that there is no difference whatever in the ratio of gravitational 
to inertial mass from one substance to another leads to the equivalence principle and 
is basic to Kinstein’s general relativistic theory of gravitation. However, it is im- 
portant to realize that the experiments that have been performed, while of very great 
accuracy, are not perfectly accurate, so that it cannot be said that the equivalence 
principle has been verified experimentally. This point has been emphasized by 
Dicke,* who noted that the energies involved in beta decay and in the self-gravita- 
tional interaction of a body of laboratory size would, even if they were to violate the 
equivalence principle, be far too small to affect the results obtained by Eétvés. In 
any event, we shall accept here the validity of the special theory of relativity, in ac- 
cordance with which inertial mass and energy are interchangeable. 

The Eétvés experiments are sufficiently accurate to show that the main factors 
that contribute to the inertial mass of a body also contribute equally or nearly 


equally to its gravitational mass.*: * To this end, we first quote formulas for the 


amount 1/64 by which the inertial mass would exceed the gravitational mass of an 
atom if one or another of these factors were to make no contribution to the gravita- 
tional mass. In this discussion, Z and A are the charge and mass numbers of the 
atom under consideration, which in practice is electrically neutral; 6A is the dif- 
ference between inertial and gravitational mass numbers; e and m are the electron 
charge and rest mass; J/ is the atomic mass unit (931 * 106 electron-volt /c?); PR is 
the nuclear radius; and ¢ and &# are the speed of light and Planck’s constant divided 


by 2x. The pertinent formulas are as follows: 
Rest masses of the atomic electrons: 
M6A mZ. (1) 


The kinetic energy of the electrons combines with the electromagnetic energy of 
interaction between electrons and nucleus to give the atomic binding energy; for Z 


not too small, § 
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Mé6A ~ —15.6 Z’ electron-volt/c?. 
Electrostatic energy of the nucleus:? 
M6A = 0.6Z(Z — 1)e?/R. 


The kinetic energy of the nucleons combines with the mesonic interaction energies 
and the electrostatic energy (3) to give the nuclear binding energy :§ 


MéA = — B(Z, A). (4) 


In making a comparison between these formulas and experiment, we obtain suffi- 
cient accuracy by using the Z and A values for the most abundant isotope present in 
each of the substances mentioned above: magnalium, copper, and platinum. Table 


| presents the values of |64/A], in units of 10-5, computed from equations (1)—(4) 


for each of these isotopes. The nuclear radius 2 is taken to be® 
R = 1.24A"* & 10-' cm. (5) 


TABLE 1 


VALUES OF Z, A, AND 6A/A COMPUTED FROM EQUATIONS (1)-(4), FoR Most ABUNDANT IsoToPE 
PRESENT IN EACH OF THE SUBSTANCES INDICATED 


{|< 105 from Equation 
SUBSTANCE Z d (1) 3) (4) 


Magnalium : 27 26,400 25 144,000 890 , 000 
Copper 2 63 25,200 38 242,000 923 , 000 
Platinum 8 195 21,900 223 397 , 000 838 , 000 


It is apparent from Table I that the magnalium-platinum and copper-platinum 
differences in each column are much greater than the experimental limit of !/». 
Indeed, if we ignore the very unlikely possibility of fortuitous cancellations, we can 
conclude that the gravitational mass of an electron is the same as its inertial mass 
within a part in 104; that the atomic binding energy makes the same contribution 
to both masses within '/; per cent; that the electrostatic repulsion between the 
nuclear protons has the weight expected from the equivalence principle to within 
two parts in 10°; and that the nuclear binding energy decreases the weight of an 
atom by the expected amount to within a part in 10°. 

The same kind of argument can be used to show that positrons have the same or 
nearly the same gravitational mass as electrons and with the same sign. This 
demonstration is of some current interest because of suggestions by several physi- 
cists, mainly in connection with cosmology, that antimatter may have the opposite 


sign of gravitational mass from ordinary matter.'° By “antimatter” is meant a 
substance in which normal particles (electrons, protons, neutrons) are replaced by 
antiparticles (positrons, antiprotons, antineutrons). Real positrons are not present 
in normal matter, so that it is necessary to make use of the virtual positrons that are 
produced because of the polarization of the vacuum by the Coulomb field of an 
atom. Since the probability amplitude for the presence of a virtual electron-posi- 
tron pair is of order Ze?/he = Z/137 (for Z not too large), we expect the mass dif- 
ference to be 

Mé6A = m(Z/137)? (6) 


if positrons have negative gravitational rest mass. This mass difference is caleu- 
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lated quantitatively in the next section; it is shown there that equation (6) has the 
correct order of magnitude and that comparison with the K6tvés experiments indi- 
cates that positrons have normal gravitational mass.!! 

A similar argument can be made in the case of antinucleons.''! '? Again, virtual 
antiprotons and antineutrons should be produced because of the polarization of the 
vacuum by the mesonic field of the nucleus. The corresponding 176A, if antinu- 
cleons possess anomalous gravitational masses, would be very much larger than 
that arising from virtual positrons, but it rests on a much less secure experimental 
and theoretical basis. The effect is expected to be larger because, first, the prob- 
ability amplitude for the presence of a virtual nucleon-antinucleon pair should be of 
order unity rather than 7/137 and, second, the mass anomaly should be of order J1/ 
rather than m. On the other hand, the theoretical structure of meson-nucleon field 
theory is not nearly so well understood as that of quantum electrodynamics, so that 
a reasonably reliable calculation is beyond the scope of present-day theory. Fur- 
ther, there is no direct experimental evidence for the existence of nucleonic vacuum 
polarization that is comparable with the evidence for electron-positron vacuum 
polarization. The latter, as produced by the atomic Coulomb field, manifests it- 
self in the Lamb shift, where its contribution can be verified quantitatively.!* For 
these reasons, we shall not attempt to compute the effect of an anomalous anti- 
nucleon gravitational mass on the Kétvés experiments, although it seems likely that 
the effeet would be large; instead, we shall devote the next section to a quantitative 
discussion of the positron case. 

We conclude this section by noting that several physicists'*: have remarked 
that a violation of the equivalence principle leads, in general, to a violation of the 
principle of energy (or inertial mass) conservation. Suppose, for example, that a 
reversible reaction (nuclear or other) can take place in which the aggregate of reac- 
tants A on one side has, say, a smaller ratio of gravitational to inertial mass than 
the aggregate of reactants B on the other side. We further assume, tentatively, 
that energy is conserved; then the inertial masses of A and B must be the same, 
since the reaction can occur. Now, if we permit the reaction A — B to take place 
at a higher altitude in the earth’s gravitational field than the reaction B — A, then 
a eyelic process could occur in which B falls and A rises, and energy would be 
steadily withdrawn from the system. 

In the particular case in which A is a nucleon-antinucleon pair and B is a pair of 
photons, Morrison argues that energy conservation can be maintained by requir- 
ing that the inertial mass (or energy) of B be increased at the lower altitude before 
it is possible for the reaction B —~ A to take place. This increase in inertial mass 
is assumed to be of such a magnitude that it just uses up the energy that would 
otherwise be withdrawn from the system. The argument thus implies either that 
the inertial mass of A is greater at lower than at higher altitudes or that the inertial 


mass of B is greater at higher than at lower altitudes. Rather than accept such a 


dependence of inertial mass on gravitational potential, this dependence necessarily 


being different for different substances, we prefer simply to place the energy-con- 
servation principle on the same basis as the equivalence principle. Both are valid 
with very great accuracy, but there is no assurance at present in either case that 


this accuracy is perfect. The idea that energy may not be precisely conserved is not 
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a new one and is, indeed, central to the continuous-creation hypothesis in cos- 
mology. '® 

Gravitational Mass of a Positron.—The system under consideration consists of an 
electrically neutral atom in interaction with the electron-positron pair field and the 
earth’s gravitational field. In the absence of gravity, the Hamiltonian is 


Hy = Ha + Hp + He, (7) 


where H4, Hp,and He are the Hamiltonians for the atom by itself, the pair field 
by itself, and the (primarily Coulomb) interaction between the atom and the pair 
field, respectively. For our purpose, the gravitational potential of the earth can be 


regarded as static and very nearly independent of the space co-ordinates; the scale 


of distance variation is of the order of the radius of the earth. We call this poten- 
tial gand choose its units such that the total energy (or rest mass times c”) of the atom 
in its ground state is changed from the value Fy far from gravitating bodies to 
Eo(1 + ¢) in the terrestrial laboratory. With this choice, the gravitational and iner- 
tial energies (masses times c”) of our system in a particular state are simply the co- 
efficient of @ in the corresponding energy eigenvalue and the ¢-independent. part of 
this eigenvalue, respectively. 

We consider here three cases, which correspond to different assumptions con- 
cerning the effect of gravity on a positron: (a) A positron behaves in the same way 
asanelectron. (6) The total gravitational energy (rest energy and kinetic energy) 
of a positron is equal to and opposite in sign from that of an electron. (c) The grav- 
itational rest mass of a positron is equal and opposite that of an electron, but its 
kinetic energy is acted on normally by a gravitational field. 

Case a corresponds to the equivalence principle. In this case the Hamiltonian 
(7) is replaced by 


H, = Hi(i+ 94), (8) 


where we assume that the variation of @ over the dimensions of the atom can be 
neglected. Then the eigenvalues of equation (8) are clearly equal to (1 + @) times 
the eigenvalues of (7), and the gravitational and inertial masses of the atom are 
precisely equal. However, as is well known, the lowest eigenvalue EH» of H, is, 
shifted by a logarithmically divergent amount, owing to the presence of Hp + He 
and this must be included as part of the renormalization of the atomic mass. This 
renormalization is, of course, the same for the gravitational and inertial masses, so 
that they remain equal. 
The simplest Hamiltonian that corresponds to case b is 


H, = (Ha + He)(1 + 6) + He + oHp'; (9) 


here Hp’ is so chosen that it assigns a positive energy to each occupied positive 
energy (electron) state, and a negative energy to each unoccupied negative energy 
(positron) state. In the momentum representation, 


Hy, = ie K,a,.*a,. — >! E,a;,a,*, (10) 


where the two summations are over positive and negative energy states £;,, as indi- 
cated, and a, and a,* are electron-destruction and -creation operators, respectively. 


The simplest choice for Hp’ is 
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H,' = >! +) Wa,*a, + pw ) E.a,a,*. (11) 


We are interested in those deviations of the eigenvalue of H, from the eigenvalue 
Ky of Ha that are of second order in H¢ and first order in ¢. A. straightforward 
perturbation calculation shows that, to this approximation (Z not too large), 


Ek, ~ El + o) — (1 + 26)>0Ep (P| H-|0)\2, (12) 


where the summation is over all eleetron-positron pair states P, of total energy Fp, 
and the matrix element of /¢ is between the vacuum and pair states. 

As remarked above, this summation is logarithmically divergent. The ¢-inde- 
pendent part of (12) is, of course, the renormalized total energy (inertial mass times 
c*) of the atom, since 7, and H, differ only in their ¢-dependent parts. However, 
the difference between the inertial and gravitational energies (that, is, the @-inde- 
pendent part of / less the coefficient of the ¢-dependent part) is no longer zero, as 


in case a, but is now equal to 


DE \(P|H¢|0)|? (13) 


and hence is logarithmically infinite. It cannot be renormalized away and can be 
interpreted only by limiting the momentum space available to the pairs in some 
way. Weshall return to (13) in the discussion of case e. 

The simplest Hamiltonian that corresponds to case ¢ 1s 


H, = (H4 + He) 1+ 6) + He + oH?" 
where 
Hp" = Hp — 2mc?>>— a,a,*. 
Again, a perturbation calculation shows that 


EB. = Ey (1+ ¢) — (+ ¢) DEp (P| He¢|0)? — 2me’o Ep | (P|He|0))?. 
(16) 
The quantity 176A of the preceding section, the amount by which the inertial‘mass 
exceeds the gravitational mass, is just the negative of the last term of (16), divided 
by c°@: 
M6A ~ 2m > Ee 2 (P\H-\0)\2: (17) 
this quantity is finite, and we now proceed to compute its magnitude. 
We substitute the interaction Hamiltonian He = ef y*Vyd*x, where V is the 
atomic electrostatic potential and y is the electron field operator, in (17), and per- 


form the spin sums. ‘The result is 
MbA = [2m/(2r)*| fidPki [dhe U(q)\? [Eikb, — h®et(k? + uw? — 
q?/4))/ {Eiki + F2)*)], (18) 


where AR, and hk» are the electron and positron momenta, w = me/h, Ey» 


he(ky,27 + p*)*!?, ¢ ki + kik '/,(k; — ky), and U is the Fourier transform of 
eV. Expressed in terms of the form factor / of the atomic charge distribution: 


U(q) = —(4nZe?/q?)F(q). (19) 
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Relation (19) between U and F is so chosen that F is equal to unity when q is 
somewhat larger inan the reciprocal atomic radius and somewhat smaller than the 
reciprocal nuclear radius; it becomes zero when qg goes appreciably outside this 
range. 

It is convenient to change the variables of integration in (18) from ky, Re to q, 2, 
y, Where 


(k,? + k.? ao 2”) y° - (2k? q°) + l oa, 
Cag = ty te") /q? 2(q: k)/q’, 
a= 1-+ (4p? q*). 


The sixth variable is the azimuthal angle of k with respect to q as polar axis and 
can be integrated over at once to give a factor 27. The x, y integrations are per- 
formed next, with the range of y from — (22 — a)? to +(24 — a)! 


2 and the range 
of cfrom!/,ato ©. The integration over q involves the atomic form factor and is 
done last. With the heip of (19) and (20), equation (18) can be written 


(4m /32*)(Z/137)? ( F(q) *flajdq/q, 
J WV 


(xv? — y*)'/F=+1—e 
t(a@) a dy Y t ’ 
(uv? — y?)"/? [ax + (a? — y?)"/?] 
where a slight simplification has been achieved by assuming that /'(q) is spherically 
symmetric. The y-integration in (22) can be performed analytically and the result 
cast in the following form: 


Ia) 
2x — a) 


1/9 
1/2 


X( (x7 — 2 


We wish to compute f(@), given by (23), for the range from a = 1(q = ©) toa= © 
(q = 0). It does not appear possible to evaluate the integral in closed form, but 
an excellent approximation over the entire range of @ can be obtained. For large 
a, the first four terms of the asymptotic expansion of f(a) are 


fla) = (9r/32a"/?) [1 + (1/12a@) + (3/128a?) + (5/512) | (24) 
The integral (23) can be evaluated exactly for a 1, and it is found that f(1) = 
1. Finally, the behavior of f(@) for @ slightly greater than unity can be found by a 
rather complicated expansion procedure; it changes linearly with a slope of ~*/,. 
Thus for @ slightly greater than, or equal to, unity, we have 


fla) = 1 — 3(a@ — 1)/4. (25) 


The leading term of the asymptotic expansion (24) is plotted on log-log co-ordi- 
nates as the dashed line in Figure 1, and all four terms as the solid curve. The dotted 
curve in the upper left corner of the figure is computed from (25) and drawn to 
join the solid curve so as to agree with a numerically computed point at a = 2. 
This combination gives an excellent representation of f(a). It can be fitted within 
the plotting accuracy by the following formula, which is chosen to have the con- 
venient form of the asymptotic series: 
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Fig. 1.—Plot of the leading term of eq. (24) (dashed line) and all four terms of eq. (24) 
(solid curve). The dotted curve in the upper left corner is computed from eq. (25) and 
is drawn so as to join the solid curve, in agreement with a numerically computed point 


ata = 2. 


f(a) a~"/?10.884 + (0.118/a) — (0.137/a?) + (0.1385/a%) | (26) 


The first term of (26) is the same as the first term of (24), and the other three co- 
efficients are chosen to fit the value and slope at a = 1 given by (25) and the numeri- 
cally computed value at a = 2. 

Substitution of the leading term of (24) or (26) in (21) gives a good approxima- 
tion for the mass difference: 

M6A = (3m/8r)(Z 137)" | F(q)\*dq/(q? + 4u?)!?, 
0 
which was quoted as equation (1) in note 11. It is apparent from Figure 1 that 
this underestimates the mass difference, but never by as much as 12 per cent. 

More accurate results than those quoted in note 11 can be obtained by sub- 
stituting (26) in (21) and choosing, for /(q), the Fourier transform of a uniformly 
charged sphere of radius R, where FP is given by equation (5). A simple expression 
for the integral in (21) is obtained if terms of order (uF)? and 1/ya are neglected 
in comparison with unity, where a is the atomic screening radius: 


} F(q)\*fl(addq/q ~ In (1/uR) + 0.338. (28) 
Jo 


Substitution of (28) in (21), with appropriate numerical values for the parameters, 
gives 6A/A for magnalium, copper, and platinum equal to 12, 25, and 53, respec- 


& 


tively, in units of 10 These numbers are somewhat larger than those quoted 


earlier! (10, 20, and 43); this is partly because the leading term of equation (24) 
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underestime: \f(@) and partly because the assumption made coneerning F'(q) in 
note 11 unde »stimates the large-qg Fourier components of the nuclear charge 
density. 

Thus, on the basis of case ¢, we would expect fractional differences of the ratios 
of gravitational to inertial masses of 13, 28, and 41 parts in 10°, respectively, for 
magnalium-copper, copper-platinum, and magnalium-platinum. Since the experi- 
ments show differences of less than !/2 part in 108, we can conclude that the gravita- 
tional rest mass of a positron is the same as that of an electron within about | per 
cent. Thus the principal result of note 11 is not altered by the more accurate cal- 
culation presented here. The remaining approximations, having to do with the 
neglect of (uf)? and 1/ua compared to unity, the simple choice for F(q), and the 
neglect of higher powers of Z/137 in the treatment of the atomic Coulomb field, will 
not significantly modify our conclusion. Further, we note that, since Zp 2 2me?, 
comparison of equations (13) and (17) shows that, no matter how interpreted, case 
b will always yield a larger value of 176A than will case e. 

The foregoing results depend on the specific assumptions made in equations (9), 
(10), (11), (14), and (15) concerning the forms of H, and H,. While these forms 
are simple, they are by no means unique. Tor example, the term ¢H¢ could be 
modified so as to alter its pair creation and destruction parts without changing its 
electron and positron scattering parts; this last requirement means that the atomic 
electrostatic binding energy included in equation (2) does not contribute to the dif- 
ference between gravitational and inertial masses. Another modification would 


. . . . . . n 
consist in adding pair creation and destruction terms to Hp’ and Hp ; these would 


not alter the effect of a gravitational field on a free positron. It is possible that 
such pair terms could mask the assumed anomalous gravitational properties of a 
positron and thus vitiate the conclusions of the last paragraph." 

If we multiply the pair creation part of ¢1¢ by some number 6 and the pair de- 
struction part by 8, in order to keep the whole Hermitian, we readily find that 


equation (12) is changed to 
Ey’ ~ Eo (1+ $) — [1 + 26Re(8))Q/ Ee ~'|(P|He|0)*, (29) 
and equation (16) is changed to 


BE! = Eyl + 6) — {1 + o[2Re(@) — 1]} DEPP \Al0 
- 2mc*o> Ep~*|(P|He(0)|?, (30) 


where Re(8) denotes the real part of the number 8. Comparison of equations (12) 
and (29) shows that the choice Re(8) '/, would completely mask the assumed 
positron gravitational anomaly in case 6. On the other hand, the form of equa- 
tion (30) shows that no such masking is possible in case ec. 

The second modification mentioned above—addition of pair creation and destrue- 
tion terms to Hp’ and H p”—leads toa more complicated calculation; contributions 
of order @Z, as well as ¢Z*, would be expected. However, it turns out that if these 
pair terms have the simple form of a sum of products of positive-energy creation 
operator and negative-energy destruction operator for the same momentum and 
spin, together with the Hermitian adjoint, there is no change in equation (12) or 
equation (16). This result is direetly attributable to the assumption of no momen- 
tum transfer in the newly added pair processes. If the spatial variation of @ were 
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taken into account, momentum transfer would occur; then the dependence on the 
radius of the earth would make such terms small and ineffective in masking. 

Other possibilities for masking of the assumed anomalous gravitational be- 
havior of a positron have not been explored. However, it seems likely, on the basis 
of the examples investigated, that it could occur only by a fortuitous combination of 
the assumed positron anomaly and a very special modification at some other point in 
the theory. Considerations of simplicity, taken together with the work of this 
section, lead us to conclude that positrons are very likely to have normal gravita- 
tional properties. 

Passive and Active Gravitational Mass.—In the last two sections we have dealt 
with the passive gravitational masses and the inertial masses of particles and anti- 
particles. We can conclude with a very high degree of certainty that the inertial 
masses of positrons and antinucleons are positive, simply because energy is required 
to create them and is released when they annihilate. Furthermore, the experi- 
ments on the creation and annihilation of real electron-positron pairs make it. vir- 
tually certain that electric charge is conserved in these processes; then the deflee- 
tion of positrons in magnetic fields also shows that they have positive inertial mass. 

A somewhat less definite conclusion can be drawn from the preceding sections: 
that the passive gravitational masses of antiparticles are also positive. One may 
then ask whether it is possible to account for a separation of matter and antimatter 
on a cosmological scale by means of gravitation.’ The only possibility lies in as- 
suming that the active gravitational masses of antiparticles are negative. As 
Bondi has shown,? such behavior may be consistent with the equivalence principle 
and the general theory of relativity, if the inertial and passive gravitational masses 
are precisely equal. All test bodies will be attracted by particles, and all will be 
repelled by antiparticles. An anomalous situation then arises when the motion of 
a particle-antiparticle pair under their mutual gravitational interaction is consid- 
ered; this case is discussed in some detail by Bondi. However, it is sufficient for 
our purpose to point out that an aggregation of antiparticles by themselves would 
behave in a normal but unstable manner, in complete analogy to an aggregation of 
electrically charged particles all of the same sign. 

Thus cosmological separation of matter and antimatter by means of gravitation 
requires that the passive gravitational masses of at least some kinds of antiparticles 
be negative. Then the equivalence principle is invalid, and general relativity can- 
not be used. While the interpretation of the E6tvés experiments discussed earlier 


in this paper probably precludes this possibility, it is of some interest to explore it 


further to see where it leads. When objects of both signs of passive gravitational 
mass are interacting, the anomalous behavior discussed by Bondi can be avoided 
by assuming that Newton’s third law of motion is valid; for positive inertial masses 
this is the same as assuming that the active and passive gravitational masses of each 
object are equal in magnitude and sign. Then if at least the most massive anti- 
particles have negative gravitational mass, matter and antimatter will tend to sep- 
arate from each other. 

An interesting consequence of this model is that gravitational experiments per- 
formed in a region of the universe in which all particles are replaced by their anti- 
particles (antilaboratory) will lead to results different from those obtained when the 
same experiments are performed in a normal region of the universe (laboratory). 
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This is because, as was shown in the first section, electromagnetic and nuclear bind- 
ing energy have positive gravitational mass. Now photons, which give rise to 
electromagnetic energy, and 7 mesons, which give rise to nuclear binding energy, 
are their own antiparticles. It follows that even if positrons, antiprotons, and anti- 
neutrons were all to have negative gravitational mass, there would be an asym- 
metry between the situations in the laboratory and in the antilaboratory. 

It may be possible to preserve complete symmetry between laboratory and anti- 
laboratory, while still having gravitational repulsion between the two, by drastically 
modifying the concept of gravitational field or by adopting an action-at-a-distance 
theory of gravitation. We could postulate that all photons, particles, + mesons, 
and their kinetic and intereaction energies attract each other; that all photons, 
antiparticles, « mesons, and their kinetic and interaction energies attract each 
other; and that all particles repel all antiparticles.'’ It seems unnecessary to at- 
tempt to develop such a model, however, since the preceding sections strongly indi- 
cate that particles and antiparticles actually attract each other. 

In conclusion, it seems very likely that all particles and antiparticles have posi- 
tive inertial and passive gravitational masses and that the equivalence principle is 
valid to at least very great accuracy. The active gravitational masses of anti- 
particles are then unknown. However, gravitational separation of stable aggrega- 
tions of matter and antimatter on a cosmological scale cannot then be achieved by 


assuming that they are negative, and the validity of Newton’s third law would re- 


quire that they be positive. 
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SPECIFIC INTRACELLULAR BINDING OF RUBIDIUM 
BY RAT DIAPHRAGM MUSCLE* 
By P. Menozzi,t D. Norman,t A. Pouuerti,§ G. Lester, AND O. HECHTER 
WORCESTER FOUNDATION FOR EXPERIMENTAL BIOLOGY, SHREWSBURY, MASSACHUSETTS 
Communicated by Theodore Shedlovsky, November 15, 1958 


Introduction.—How cells discriminate between K+ and Nat to maintain an 
unequal distribution of ions, potassium being high inside the cell and sodium 
outside, is a fundamental problem in biology. Currently, it is widely accepted that 
unequal ion distribution in muscle and nerve is the resultant of the activity of the 
cell membrane, postulated to possess energy-requiring pumps which transport 
K+ into and Na* out of the cell against electrochemical gradients.!. To account 
for ionic selection by these “active-transport”’ processes, the existence of specific 
carriers which discriminate between sodium (or lithium) and potassium (or rubid- 
ium) ions is postulated. It has alternatively been proposed that the selective ac- 
cumulation of potassium is not primarily due to the ‘“‘membrane”’ but is the resultant 
of intracellular mechanisms, such as “fixed charges’? of structural proteins,? a 
charged-ordered semirigid lattice structure*® or organized ‘water lattices.’’4 

In a preceding paper? it was shown that the selective accumulation of rubidium 
against an apparent concentration gradient in Neurospora crassa is not explicable 
in terms of a membrane-pump mechanism but appears to involve a specific eyto- 
plasmic binding system. This paper likewise presents evidence that the cytoplasm 
of rat diaphragm muscle exhibits properties of selective ion binding. Specifically, 
it will be shown that (a) the preparation of the usual hemidiaphragm involves cut- 
ting each muscle fiber of the diaphragm at two ends; (6) associated with this, a 
characteristic permeability feature of the cell is markedly altered, so that ions 
which are excluded from intact fibers readily penetrate into the cut fibers of the 
hemidiaphragm preparation; (¢) under these conditions Na+ and Rb* both enter 
the cell water of the cut fibers, but only the latter ion is accumulated against an 
apparent concentration gradient; and (d) the accumulation system for rubidium, 


which appears to operate with potassium as well, probably cannot be explained in 


terms of a membrane-pump mechanism. 
Materials and method.—-Two kinds of diaphragm muscle preparations, from male 
rats (Sprague-Dawley strain) weighing 100-180 gm. fasted 18-24 hours, were em- 
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ployed in these studies. One preparation was the isolated sectioned hemidia- 
phragm,® widely employed to study insulin action in vitro, where the diaphragm 
muscle is cut from the ribs on one side and the connective tissue attachments to the 
central tendon on the other, to yield, after trimming, two rectangular muscle 
strips. Since the muscle fibers run from the rib to the central tendon, each muscle 
fiber in the hemidiaphragm is cut at both ends and in our studies is about 5-8 mm. 
long. The other preparation, consisting of relatively uninjured fibers, was ob- 
tained by removing the diaphragm together with a section of the rib cage and the 
central tendon, in a manner similar to that deseribed by Kipnis and Cori,’ except 
that we included a larger portion of the rib cage, to minimize the possibility of in- 
juring fibers of the diaphragm muscle. 

The muscle preparations were incubated in Krebs-Hensleit Ringer’s phosphate 
buffer’ at pH 7.4 (henceforth referred to as “normal Ringer’s’’), potassium-free 
tinger’s, where Rb~*§ (as RbCl) replaced K*, or in modifications of these solutions, 
which will be deseribed. Unless otherwise stated, incubation was carried out 
at 37.5° (with shaking) in the absence of added glucose, with O» as gas phase. 
With sectioned hemidiaphragms the volume of incubation medium employed 
was generally 5 ml.; cage preparations were generally incubated in 50 ml. of 
medium. 

Experimental.—Cutting of muscle fibers results in marked changes in the per- 
meability characteristics to various ions and sugars. ‘T’o localize the basis for these 
alterations, we have examined the penetration of ferrocyanide into hemidiaphragm 
and cage preparations by histological means. When hemidiaphragm or cage prep- 
arations are incubated for 1 hour in a Ringer’s solution, modified by substituting 
sodium and potassium ferrocyanide for NaCl and KCI, subsequent histological ex- 
amination shows marked differences in ferrocyanide distribution between these 
preparations.? In the hemidiaphragm it was observed that ferroeyanide had pene- 
trated into every fiber from both cut ends, leaving a central zone which had not 
been penetrated; in cage preparations ferrocyanide did not penetrate into any of 
the fibers and was restricted exclusively to the intercellular ground substance.!° 
Study of ferrocyanide penetration in hemidiaphragms cut at only one end (the ten- 
don was left in place at the other end), incubated for varying intervals of time, re- 
vealed that the entry of this large anion occurred exclusively via the cut ends and 
increased progressively with time at a rate estimated to be 1.8 mm. per hour from 
each cut end, over a 2-hour period. 

Mannitol and sucrose, which are excluded from intact fibers, likewise penetrate 
into the intracellular water of the cut fibers of the hemidiaphragms, presumably via 
the cut ends. Thus the “space”’ of mannitol and sucrose was found to be 0.71 and 


0.50 ml/gm muscle (wet wt.), respectively, whereas the inulin “space” in hemi- 


diaphragms is 0.26 ml./gm muscle. The latter value is similar to those obtained for 


the intercellular volume of rat diaphragm in vivo, as determined by mannitol,” 
raffinose,’ and inulin,!? as well as to our own results in the cage preparation, where 
the ‘‘space”’ for sucrose and mannitol is equivalent to 0.24—0.28 ml/gm muscle. It 
would thus appear that the fibers ina sectioned hemidiaphragm possess two cut ends 
which permit non-selective diffusion of molecules like sucrose and ions like ferro- 
cyanide but still restrict the entry of a large molecule such as inulin. 

Despite these findings, which indicate that the cell membrane is no longer a con- 
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tinuous envelope in the fibers of the hemidiaphragm, the property of “active trans- 
port” of ions is retained. While hemidiaphragms incubated in Ringer’s solution 
lose intracellular potassium and gain sodium to a much greater extent than do cage 
preparations, an equilibration of these ions between cell water and the external 
medium does not occur in cut fibers; the intracellular concentration of sodium is 
maintained, to a variable degree, below that of the medium, while the intracellular 
concentration of potassium is higher than that of the medium.!* Following incu- 
bation in normal Ringer’s solution, intracellular sodium gain is roughly equivalent, 
on the average, to potassium loss. This apparent reciprocity between potassium 
accumulation and sodium exclusion could suggest a primary causal relationship 
between these processes, as has often been proposed. However, these processes 
can be dissociated. Table 1 shows that, when hemidiaphragms are incubated in 
Ringer’s solution where the external Na* concentration is varied from 10 to 148 
mM (substituting sucrose for NaCl), potassium loss is quite constant, independent 
of whether or not sodium is gained. When the external Na+ concentration is 


TABLE 1 


Na GAIN AND K Loss in HEMIDIAPHRAGMS* 


VOLUME OF 

DISTRIBUTION OF 

INTRACELLULAR - Na IN INTRACEL- 
Na Gaint LULAR WATER 


(CONCENTRA 
TION OF 
Nat In No. o1 
M epiumM H eMIDIA- K Losst 


(uM/ML1) PHRAGMS uM /Mu Ceti Water) (Per CENT) 
148 38 110 + 2.9§ +81.5 3 
50 6 108 + 4.0 +17.! 2 

] 
| 


0§ 66 + 2.0§ 
2 72+ 69 
{ 66 + 4.1 
0 


25 6 inl = 3.9 — 3.4 
100 + 8.0 


10 6 120 = 2.0 -10.0 + 
Incubated at 37.5° for 1 hour in Ringer’s solutions modified by substituting sucrose for the NaCl removed. 
Initial K was 170 + 2.0 uwM/ml cell water ‘ 


Initial Na was 20 + 1.1 «6M /ml cell water. 
Standard error of the mean. 


varied from 25 to 148 mM, the mean intracellular concentration of sodium achieved 
is equivalent to equilibration of Nat in about 70 per cent of the cell water. The 
sodium-exclusion process operative in cut fibers is exceedingly variable among in- 
dividual tissues exposed to similar conditions; the basis for this variability has not, 
as yet, been elucidated. 

Sodium exclusion in cut fibers appears to be an “active” process, in that it is in- 
hibited by dinitrophenol (0.2 mM) or low temperature and to this extent resembles 
a “sodium pump.” It should be noted from Table 1, however, that if a “sodium 
pump” actually were operative, the activity of this pump, but not its capacity, 
must be assumed to vary directly with the Na* concentration of the medium. 
Alternatively, it is possible to explain sodium exclusion on the basis that a certain 
fraction of the cell water is not available for Na* and that sodium gain represents 
an equilibration process in an available fraction of the total cell water. On this 
latter view, the maintenance of cell water in a state unavailable to Na* would re- 
quire energy; a mechanism for sodium exclusion along these lines has been suggested 
by Baird et al.4 Our present data do not permit a definitive decision between these 
dissimilar alternatives. 

The “active-transport” properties of cut fibers are clearly revealed when the up- 
take of rubidium by hemidiaphragms is examined. Table 2 shows that hemidia- 
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TABLE 2 
Tue ACCUMULATION OF RB IN HEMIDIAPHRAGMS* 
Incubation Time Intracellular Rb 

A11n.) ue ) 

15 23.5 2& 1.17(21) 

30 +1 + 0.2 (75) 

60 62 + (0.4(118) 

120 65.0 + 4.0 (17) 


* Incubated at 37.5° in K t-free Ringers containing 5.1 mM Rb*. 


+ Standard error of the mean; the figures in parentheses refer to the number of determinations. 


phragms incubated in potassium-free Ringer’s containing Rb* (5 mM) accumulate 
rubidium against an apparent concentration gradient, so that, after 1 hour, the 
ratio of Rbj/Rb, is about 12:1. Rubidium accumulation can be inhibited by addi- 
tion of K* to the Ringer’s solution containing Rb*, as shown in Table 3, and it ap- 
pears that the inhibitory effect of external K+ is competitive in nature. As with 


TABLE 3 
tT oF EXTERNAL K* CONCENTRATION ON RB ACCUMULATION 
IN HEMIDIAPHRAGMS* 


INTRACELLULAR Rp 
uM/ML. Cell Water) 
MeEpDIUM Excess 
‘ONCENTRATIO over 
th HEMIDIA- External PeR CENT 
uM/M)) PHRAGMS Concentration INHIBITION 
5 38 


0 
19 


2 
o 


3 
9 


0 


m DS OS ee OF 
aE ES 


) 
) 
) 


l 
I 
l 
l 
I 
l 


5 
20 5 
122 5 


* Incubated at 37.5° for 1 hour. 


ft Standard error of the mean. 


other active-transport systems, rubidium accumulation in excess of the external 
coneentration is no longer observed when tissues are incubated in the cold (at or 
near 0°) or with dinitrophenol (0.2 mM) at 37.5 

The rubidium-accumulation system in cut fibers has been sufficiently character- 
ized to permit some conclusions concerning its nature. Thus rubidium accumula- 
tion is not dependent upon a ‘‘sodium pump,” which, by expelling Na+ against a 
gradient, results in the entry of Rb* against a concentration gradient. Figure 1 
shows the results obtained in individual tissues where rubidium accumulation is 
plotted against sodium entry, measured in terms of equilibration of sodium in the 
cell water. It can be seen that the accumulation of rubidium is not dependent 
upon Na* exclusion. Indeed, the data demonstrate a tendency for rubidium aec- 
cumulation to increase as the effectiveness of the sodium-exclusion process de- 
creases. 

It might be possible to account for rubidium accumulation by assuming a mem- 
brane mechanism which specifically pumps Rb*+ (or K*) into the cut fiber at a rate 
sufficient to maintain an intracellular concentration of free Rb*, 8 to 15 times in 
excess of the external concentration, despite efflux from the cut ends as well as the 
membranes. On this basis, it might be expected that, as the number of cut ends 
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RELATION BETWEEN Rb ACCUMULATION 
AND Na EXCLUSION IN HEMIDIAPHRAGMS 
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Fic. 1.—Each dot represents an individual experiment in a 
hemidiaphragm incubated for | hour at 37.5° in K *-free Ringer’s 
containing 5.1 mM Rb 





per fiber length are increased progressively (decreasing the ratio of the surface area 
of membrane to cut end), rubidium accumulation should decrease proportionately, 
since the cut ends may be regarded as offering no barrier to the efflux of free Rb*. 
As shown in Table 4, when hemidiaphragms (2 cut ends) are further cut perpendic- 
ular to the fibers by a razor blade to give halves (4 cut ends) or thirds (6 cut ends), 
the expected proportional decrease in rubidum accumulation does not occur. 


TABLE 4 
Krrecr or AppITIONAL CurtTinG OF HEMIDIAPHRAGMS Upon Rs ACCUMULATION * 


Approx 
Fiber Length Intracellular Rb 
Preparation Mim.) uM MI. Cell Water) 
8 ‘7.0 2 1.717 C36) 
1 o1.0 & 2.708) 
Ae 30.4 + 1.97(8) 


“4 cut”’ 2.6 


“6 cut” 7 


“2 cut” (normal hemidiaphragm 5- 
5 


* Incubated for 1 hour in K-free nger’s containing 5.1 mM Rb 
t Standard error of the mean e figures in parentheses refer to the number of observations. 


These observations indieate that a Rb* pump mechanism is unlikely, unless one 
assumes that progressive cutting of the fiber results in an increase in the activity or 
efficiency of such a membrane pump. 

If the accumulation and maintenance of high rubidium levels in the cell are 
actually the resultant of a rubidium pump, then factors which completely inhibit the 
uptake of rubidium in excess of the external concentration should have an equivalent 
effect on maintenance of previously accumulated rubidium. Nevertheless, as 
shown in Table 5, cold can prevent the initial aeeumulation of rubidium by hemi- 
diaphragms, but, once rubidium has been accumulated intracellularly by prior ineu- 
bation at 37.5°, it is retained intracellularly for several hours at 0°, so that the 
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TABLE 5 


KeFECT OF COLD AND DINITROPHENOL ON ACCUMULATION AND MAINTENANCE OF 
INTRACELLULAR RB IN HiMIDIAPHRAGMS 
Intracellular Rb 
Incubation Conditions M/M1) Rbi/Rbo 

Rb Ringer’s* at* 37.5° 
60 min. 56 E 2.2T(11 
120 min. j 6.1 (6 

Rb Ringer’s at O° 
60 min. - 0.: 
120 min. Oo! 

Rb Ringer’s + DNP (0.2mM) at 37.5 

60 min. 

120 min. 

60 min. incubation at 37.5 
Plus 5 min. at O0°f 
Plus 15 min. at 0 
Plus 60 min. at 0 
Plus 120 min. at 0 

60 min. incubation at 37.5 
Plus 5 min. DNP (0.2 mM 
Plus 15 min. DNP (0.2 mM): 
Plus 30 min. DNP (0.2 mM 
Plus 60 min. DNP (0.2 mM 


HHH HH HH 
Dm DH K 


HH HH 
a 


*5.1 mM Rb 


+ Standard error of the me t figur parenthe refer to the imber of rminations 
+ 
t 


Incubated in ice batl 


ratios of Rbj/ Rb, are about 7:1 and 5:1 at | and 2 hours, respectively. Dinitro- 
phenol (0.2 mM), on the other hand, causes a more rapid exit of previously ac- 
cumulated rubidium than does exposure to cold. The dissociation of rubidium ae- 
cumulation at low temperatures into separate energy-requiring processes for “for- 
mation” and ‘‘maintenance”’ excludes the possibility that a pump mechanism is 
‘esponsible for rubidium accumulation in cut fibers. 

The active-transport mechanisms which have been considered in relation to 
rubidium accumulation are based on the assumption that ions inside the cell are in 
the same “state”? as in the external medium. Since these mechanisms do not ade- 
quately explain the observed findings with rubidium, some other mechanism must 
be involved. Asan alternative, it may be proposed that rubidium accumulated in 
the cell in excess of the external concentration exists in an “associated” or “bound” 
state rather than as the free ion. Such a hypothesis, which emphasizes the im- 


TABLE 6 


ErrFECT OF EXTERNAL RB* CONCENTRATION UPON ACCUMULATION 
or Rp BY HEMIDIAPHRAGMS* 

EXCESS 

OVER 

MepIuM CONCENTRATION oO. OF NTRACELLULAR EXTERNAL 
uM/Mi MII Ri CONCENTRA 

Rb -HRAGM Porat TION 

uM /Mu. Cetin Warer 

0O.5t 15.6 

0.4 57.3 
3.0 2.U 
3.9 3.3 
o.] S85 

* Incubated at 37 


+t Standard error « 
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portance of cytoplasmic factors rather than the ‘‘membrane,” would be consistent 
with the observed accumulation of rubidium by fibers which lack an intact cell 
membrane. The characteristics of the selective accumulation system for rubidium 
may be interpreted in terms of a selective cytoplasmic binding system. Table 6 
shows that hemidiaphragms accumulate a relatively constant amount of rubidium 
in excess of the external concentration, when the Rb* concentration of the medium 
is varied from 10 to51mM. The binding system may thus be assumed to involve 
a limited number of sites (B), which can combine with Rb+ to form RbB. The 
number of these sites may be quantitatively designated as the number of uM Rb 
per ml. cell water, in excess of the external concentration. In cut fibers there thus 
are about 70-90 sites per ml. cell water; whether or not all these sites are occupied 
depends on the external concentration of Rb*. It may be assumed that B also 
interacts with K+ to form KB, since the inhibition of rubidium accumulation by 
K+ appears to be competitive. If K+ and Rb*+ compete for sites, the question 
arises whether the rubidium accumulation observed in K+-free Ringer’s solution 
containing Rb* represents an exchange of external Rb* for intracellular potassium 
associated with B, or occupation of “vacant” sites, made available as a consequence 
of potassium loss from the injured fiber. Preliminary experiments indicate that 
rubidium accumulation occurs as a resultant of both processes. 

The specific binding system is characterized by a marked sensitivity to low tem- 
perature and dinitrophenol, both of which completely inhibit the accumulation of 
rubidium; this suggests that the formation of RbB has an energy requirement. 
Maintenance of RbB likewise appears to have an energy requirement, since pre- 
viously accumulated rubidium dissociates, to different degrees, when dinitrophenol 
is added or when the system is chilled to 0°. The fact that, at 0°, intracellular 
rubidium is maintained for 2 hours against an apparent gradient, whereas initial 
accumulation of rubidium is completely blocked at this temperature, indicates 
that the energy requirement for the maintenance of RbB is less than that required 
foritsformation. The energy requirements of the binding system may be assumed 
to be relatively modest. The association of Rb* with B proceeds at 37.5° in the 
absence of an exogenous substrate, and addition of glucose to the medium has no 
significant influence; moreover, fluoride (1.0 mM) and cyanide (1.0 mM), both of 
which interfere with cellular metabolism, have little effect. Azide (1.0 mM), while 


possessing definite inhibitory action, does not abolish the association of rubidium 


with B as does dinitrophenol. 

The accumulation of rubidium in rat hemidiaphragms has certain characteristies 
which are similar to those observed in NV. crassa.6 Thus, in both eases, it can be 
postulated that either Rb+ or K+ interacts with a specific binding system and 
there is competition between these ions for binding sites. The major difference be- 
tween these two systems appears to relate to the stability of the RbB or KB, once 
it is formed. In Neurospora, RbB or KB does not appear to be readily dissociable 
or available for exchange; in diaphragm muscle, however, dissociation of RbB is 
readily effected by inhibitors such as dinitrophenol, azide, or, to a lesser extent, by 
lowering temperature. Moreover, as previously indicated, Rb*+ exchanges with a 
fraction of intracellular potassium which may be assumed to exist as KB. It is 
possible to explain these differences between muscle and N. crassa by assuming dif- 
ferent types of B, as has been previously discussed.® 
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Our present information concerning the nature of the accumulation system in 
muscle is limited to the finding that it appears to be labile. Freezing and thawing 
or heating to 50° for a few minutes results in irreversible destruction of the system. 
Also, injury associated with cutting fibers reduces activity. There is insufficient 
information to permit profitable speculation upon the fundamental mechanism in- 
volved in the postulated intracellular binding system for Rb* or K+. The limited 
data available could fit any hypothesis involving a metastable system permitting 
selective association of K+ with binding sites, without necessitating the exclusion of 
Nat. To this extent, the present data are consistent with either the theories of 
Ling? or those of Simon and Shaw.? 

The results presented in this paper have an obvious bearing on the broad prob- 
lem of the active transport of K*+ and Nat. If the present findings on rat dia- 


phragm muscle and the previous findings on Neurospora,> which indicate that most 


of the intracellular potassium may be selectively bound by a cytoplasmic system, 
should prove to be applicable to cells generally, certain widely accepted views con- 
cerning the role of the membrane in permeability processes may require basic re- 
vision.!4 
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After incubation of cage preparations for 2 hours, ferrocvanide penetrates a small number of 
fibers, indicating that, with prolonged incubation under our conditions, a deterioration of a cellular 
barrier to ferrocyanide occurs in some fibers. 

1 ©, R. Park, L. H. Johnson, J. H. Wright, Jr., and H. Batsel, Am. J. Physiol., 191, 13, 1957. 
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tracts and in the incubation medium were measured as follows: Rb** was measured with a flow 
gas counter with suitable corrections for self-absorption; Na* and K* were determined with a flame 


photometer. The results for Rb, Na, and K are expressed in terms of «M per ml. cell water cal- 


culated on the basis of the following considerations. The total water content of both types of 
muscle preparations was found to be 0.80 ml/gm muscle. The inulin space in hemidiaphragms 
(0.26 ml/gm) and the sucrose space in the intact cage preparation (0.24 ml/gm) were taken to rep- 
resent the intercellular volume in the two preparations, respectively, and were assumed to be in 
equilibrium with the external medium. 

After this paper has been submitted, E. Bozler, M. E. Calvin, and D. W. Watson (Am. J 
Physiol , 195, 38, 1958) reported studies which indicate that a large proportion of the intracellu- 


lar KC of the smooth muscle of frog stomach exists in a “‘bound”’ state 
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ABOUT THEORIES OF COLOR VISION 
By Yves LEGRAND 
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In this paper I shall try to give some kind of introduction to the following papers, 
which are dealing with recent experimental developments in the study of color 
vision. Color vision is a very old chapter of science, and it may be interesting to 
summarize briefly its main historical developments. 

It seems that the French philosopher Descartes was the first, about the year 1630, 
to seek a physical explanation for colors. Before him, even those who thought that 
light was a modification of matter did not suppose that light might be endowed with 
any coloring property; they imagined color as a metaphysical quality of objects, a 
quality which could be perceived only when light fell on the object. but was other- 
wise independent of light. Descartes understood that the color of an object is a 
change of the light after its reflection on the object, and he thought that this change 
might be related to the “spin” velocity of light particles. 

On the basis of experiments on the mixing of pigments, several writers in the 
seventeenth century thought that all colors were a mixture of only three of them 
say, red, yellow, and blue. This “trichromacy” was also derived by Mariotte from 
un examination of the solar spectrum, but he did not try to explain that curious 
fact. One century later, in 1756, a Russian physicist, Lomonosow, put forward the 
first theory of the three principal colors. He supposed there were three kinds of 
light particles, but this idea does not seem to have been widely accepted, and 
D)’Alembert wrote in 1764 that there is a nearly infinite number of colors and not 
three, “‘as has been supposed by several physicists.” 

So trichromacy was in the air long before the beginning of the nineteenth century ; 
but the great merit of Thomas Young was to establish in his celebrated 1802 paper 
that trichromacy was not a physical property of light but a physiological one of the 
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retina. The reason given by Young is very simple: he was convinced of the truth 
of Huygens’ undulatory theory of. light. In Newton’s emission theory, it was 
quite simple to apply some form of “quantification” to the number of kinds of light 
particles, but it seemed impossible to do so at that time with the undulation theory, 
in Which wave length is evidently a continuous variable. ‘As it is almost impos- 
sible,” Young writes, “to conceive each sensitive point of the retina to contain an 
infinite number of particles, each capable of vibrating in perfect unison with every 
possible undulation, it becomes necessary to suppose the number limited, for in- 
stance, to the three principal colours, red, yellow and blue.” As regards the choice 
of “principal colours,” Young was, five years later, as a result of Wollaston’s de- 
scription of the spectrum, to modify this in favor of red, green, and violet. 

Young’s hypothesis was not followed up until it was revived by Helmholtz in 
1852, and Helmholtz remarked that the choice of the “principal colours”? remained 
arbitrary to a large extent. It was also under the influence of Young’s ideas that 
Maxwell undertook in 1857 the measurements which gave to the trichromacy of 
vision an experimental foundation, though it was already implied in a rule given by 
Newton in 1704 for predicting the results of color mixtures. But Newton thought 
that his rule was not exactly valid, and, in fact, trichromacy was firmly established 
by Helmholtz and Maxwell a century ago. 

The experimental fact of trichromacy may be expressed as follows: given any 
four lights, it is always possible to place three of them in one half of a photometric 
field and one in the other half, or else two in one half and two in the other, and, by 
adjusting the intensities of three of the four lights, to make the two halves of the 
field appear indistinguishable to the eye. This ability to match perfectly mixtures 
which are spectroscopically very different is a very striking fact in color vision. 

These “matching experiments” are the basis of colorimetry. It is important to 
understand that in such measurements the observer is not asked what he is seeing 
but only which stimuli are equivalent to him: the eye is used only as a null instru- 
ment for establishing equalities. The mathematical form of any color theory re- 
stricted to the explanation of these matching experiments is very simple, owing to 
the algebraic character of the quantities used to define the stimuli and to the linear- 
ity of equations between them. These properties are known as “Grassmann’s 
laws.’ The simple basis of all color theories obeying Grassmann’s laws may be 
stated as follows: In the retina there are at least three different kind of receptors; 
any number greater than three, say n, is possible if we suppose that there may be 
n —3 additional linear relations between the m responses. The simplest possibility 
is evidently n 2 

The first question is: What may be the physical basis for this spectral difference 
between retinal receptors? The most probable answer is that these differences may 
be due to various photosensitive pigments. Another possibility would be the exist- 


ence of one single pigment but of different physical filters (colored oil droplets 
like those of birds and tortoises, or laminated structure of the outer limb producing 


interference filters); but the multiplicity of pigments is in accordance with the 
broad variety of visual substances for which there is evidence from studies of retinal 
extracts in different animals. 

The second question is how to explain the n —3 linear relations between the re- 
ceptors. It is easy to understand why « linear law governs the initial action of 
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light on any photosensitive substance, but linearity is lost as soon as a nervous 
message is sent from the cone to the bipolar cell; then it seems necessary to suppose 
that our n —3 relations are of a physical, not a physiological, kind. The simplest 
way of obtaining such relations by physics is the mixture of three different photo- 
sensitive pigments only; if the maximum absorption of each pigment remains low, 
the relations are very nearly linear. 

It. is easy to understand why all the experimental laws of matching experiments 
result immediately from the existence of three pigments, whatever the mixture made 
eventually with them in the retinal cones may be: equivalent stimuli are such that 
the three initial photochemical actions are identical, and anything which oecurs 
afterward cannot alter this identity. Colorimetry may be affected neither by any 
mixture of pigments nor by nervous phenomena. 

If three photopigments exist in the retina, it would be very important to deter- 
mine their absorption curves. But, unhappily, matching experiments alone are in- 
adequate to determine a unique set of 3 spectral sensitivities which would corre- 
spond to the spectral absorptions of the pigments, after correction for selective 
losses through the ocular media or the macular pigment. Owing always to linear 
Grassmann’s laws, any trilinear transformation of such a set gives another set which 
is exactly equivalent from the colorimetric point of view. In practice, a very arti- 
ficial but convenient set has been chosen by the CIE under the name “XYZ sys- 
tem.” In this system of spectral sensitivity curves (Fig. 1) the ¥ curve is identical 
with the photopic efficiency curve of 
cones; the * curve has two maxima 
and a minimum, almost zero, near 500 
mu; it does not indicate any particular 
characteristic of vision and comes from 
the arbitrary choice of this wave length 
for the nearest approach of the side YZ 
of the triangle of reference stimuli to 
the spectrum locus. 

In the past, a great number of at 
tempts have been made to add supple 
mentary conditions that might limit the 
infinite variety of possible spectral sen 





sitivities of the retinal photopigments. 
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These conditions were sought in a ied 


study of dichromatism, — assuming, 

with Helmholtz and his pupil Konig, that the reduction from trichromatism 
to dichromatism results from the suppression of one kind of pigment. For the 
protanopie sort of dichromatism, this hypothesis is probably true, but for the 
deuteranopic one it is rather uncertain, and Fick’s interpretation seems to be 
better: the reduction would take place by unification of the responses of two kinds 
of retinal receptors which are independent for a normal observer. If so, the study of 
dichromatism is not able to give sufficient information about the actual absorption 
curves of the photopigments, and the only ways to obtain these curves should 
be the extraction and the chemical study of the pigments, as Dr. Wald did in his 
celebrated work on iodopsin, or the physical study of light absorption by pig- 
ments in situ, as in Dr. Rushton’s beautiful experiments. 





PHYSICS: Y. LEGRAND Proc, N. A.S8. 


We have said that the simplest way of explaining matching experiments was 
to suppose the existence of three photopigments, colorimetry being unaffected 
by any mixture of these pigments or by nervous interactions. But in theory it 
would also be possible to suppose a number of pigments larger than three, 
if supplementary conditions of linkage between the responses were assumed. 
Brindley has shown that this last hypothesis is in conflict with experience, by 
studying the breakdown of color matches in the case of intense adaptation. There 
is direct’ evidence, given by Rushton, that adapting lights, if bright enough, 
can bleach away a large fraction not only of the rhodopsin in the reds but also of 
the cone pigments of the retina. In fact, conspicuous changes in the appearance 
of colors can be produced by adaptation to colored lights without upsetting color 

matches between lights of different 
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pensation mechanisms, which are highly 


} 
Ae 
fe) ~- 

| 

1 


_- 


improbable from a physiological point 
of view. 
reap 
o 


Al If luminance has only a small effect on 
| “$10G 


1.5 9 hasa much larger one in the extrafoveal| 


+ 
| 
| 








matching experiments in the fovea, it 
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log 9 adapting luminance retina at low levels, because the scotopic 
Fic. 2 efficiency curve of rods differs from the 
photopic one of cones (Purkinje shift). 
As Brindley has shown, this is a direct proof that rhodopsin is not one of the three 
pigments of color vision. 

In conclusion, all work done in matching experiments and their breakdown at very 
high or very low levels of illumination is in good agreement with the very simple 
assumption of three, and only three, independent photopigments in the human 
retina. Ifn = 3, does it mean that there are also three types of cones in the human 
retina, each of them containing one of the three pigments? Not necessarily. In 
fact, apart from the regular change of form from the center of the fovea to the ex- 
treme periphery, all the cones seem to be alike, whereas one might expect to observe 
anatomical differences either causing or caused by the difference in the pigment con- 
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tained in the cones. There is another argument against the three types of cones, 
given by Abney’s law of luminance additivity: it is not possible to suppose that cer- 
tain bipolars or ganglion cells may add the impulses of the three types of cones 
and so respond to every wave length, because in these combined impulses any linear 
characteristic is lost. On the contrary, a fourth type of cone containing a mixture 
of the three pigments would give a very simple way of interpreting Abney’s law; 
if this fourth type is the most frequent, it also explains the fact, known from Aubert 
onward, that a very small colored object on a neutral background of the same lumi- 
nance appears colorless and, more generally, that color vision is fully developed with 
large field tests only. In the three-cones theory, on the contrary, the opposite effect 
should occur, and a white point should look colored in the neighborhood of the foveal 
threshold because it would excite only the cone with the lowest threshold. 

If it seems necessary to suppose at least four types of cones, the simplest assump- 
tion about the mixing of the three pigments is that it is made at random. Then 
the majority of the cones would contain the three pigments mixed approximately 
in equal concentrations, and their spectral response would be roughly the photopic 
efficiency curve; a smaller number of cones would contain a relatively high con- 
centration of one or two pigments, and it would be possible to classify them in six 
groups, although, in fact, their spectral sensitivity varies in a continuous manner. 
This hypothesis would seal the reconciliation between Young and Newton, between 
the sacred numbers 3 and 7, three pigments and, roughly speaking, seven kinds of 
cones in the retina! 

As matching experiments are of no value in teaching us anything on this subject 
of retinal receptors, it is necessary to find another line of approach. The first may 
be given by electrophysiology, and the thresholds for spike discharges obtained by 
the microelectrode methods have led to Granit’s familiar dominator and modulator 
curves. Such results are unhappily difficult to interpret, because they reflect 


neural interactions which modify the spectral responses. 


A second kind of information is given by the size of color discrimination steps. 
Let us reeall Dr. MacAdam’s well-known ellipses at a constant level of luminance, 
which are of great practical value but rather difficult to interpret theoretically, be- 
cause we do not know at what part of the visual pathway this form of discrimination 
is obtained. Dr. Stiles’s results obtained by the two-color method are simpler, 
but here also the assumption made that the threshold of each mechanism is inde- 
pendent of the activity of the others may not be valid, perhaps, when different 
mechanisms are very differently active (such interactions between sound and light 
have been observed); nevertheless, Dr. Stiles’s experiments give us the hint that 
nis at least 5. 

A third line of approach is obtained through artificial dichromatism (and even 
monochromatism) produced by adaptation to very bright lights. The most interest- 
ing result is obtained through adaptation to intense yellow light because, besides 
the loss of color discrimination, there appears a severe defect of luminance diserim- 
ination. This seems to indicate that the ‘blue’ pigment plays only a minor role 
in the luminance receptor where the three pigments are mixed. Another difference 
between the “blue” cones and the others is the low visual acuity associated with the 
blue mechanism, perhaps because synaptic convergence is greater for these recep- 


tors. 
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In conclusion, if colorimetry is very simply and completely explained by the 
existence of three photopigments, the nature and even the number of different ret- 
inal receptors is still unknown; but there are some reasons to think that this num- 
ber is greater than three. 

Finally, it remains to explain how colors look to us, and this psychological part 
of color theories is perhaps the most difficult to place on a physical and physiologi- 
cal ground. Young’s trichromatic idea meets here with its two main difficulties; 
the first is the nature of yellow sensation, which does not at all appear as an addition 
of green and red; the second is the total neutralization of the color attribute of 
sensation, in complementary pairs, for example. Both necessitate some kind of 
inhibition, which may be found perhaps in the retina itself. We may imagine a com- 
pound receptor containing a mixture of green and red pigments, the excitation of 
which by yellow light would deliver a special message to the optic nerve; but it is 
necessary to suppose also that in “green and “red” receptors, the messages cancel 
one another when they are equal. An inhibition of the same kind would suppress 
the responses of red, green, and blue cones when they are in equilibrium and leave 
only the message of the white compound receptor, as in Pieron’s theory of “recep- 
tor tetrad.” All this is physiologically possible, but far from being proved. 

More than eighty years ago, Hering proposed a general principle of retinal in- 
hibitions, produced by opposite actions of light on three substances; it is difficult to 
find a physiological basis for such opposed responses given by one retinal mecha- 


‘ 


nism; perhaps it might be looked for in the duality of the “fon” and “off” responses 
that certain recent researches in electrophysiology ascribe to a single retinal recep- 
tor. 

Another possibility is postulated in ‘‘zone” theories, such as those of Miiller and 
Judd, which combine Young’s peripheral trichromacy with more complicated struc- 
tures at another stage of the visual pathway. The difficulty is that a large number 
of such theories may be constructed with any of these psychological data but that 
none of them can be proved as being the right one. It would be necessary to know 
more of the mechanism of vision in the cortex because experimental evidence ob- 
tained at this level, say by electrical records, has until now been rather poor. 

Some fifteen years ago, LeGros Clark thought that the geniculate body might 


play an important part in color vision, its six layers being related with trichromacy, 
but now his ideas seem to be abandoned, or nearly so. In the visual cortex of the 


cat, Chang had registered electrically a complex response to the excitation of the eye 
by a flash of light; here again he had thought that this fact was related to trichro- 
macy, but unhappily we are now almost certain by psychological experiments that 
the cat does not perceive colors. 

In fact, the only thing which is quite sure is the extraordinary complication of the 
visual cortex; each optic nerve contains about one million fibers, while the visual 
cortex has near 1,400 million neurons. This multiplication is specially high for 
the foveal center, each fiber which comes from the fovea being connected with 5,000 
cortical neurons, according to Sholl. Therefore, it is easy to imagine the complex- 
ity of interactions and inhibitions that may arise in so intricate a network of neural 
synapses. It is easy to imagine it, but it has been impossible until now to apply 
any cortical theory even to the simplest psychological laws of color vision. 

For example, color contrast is perhaps explainable through retinal adaptive effects 
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at the line of separation between two areas the colors of which differ, but there exist 
also opposite phenomena of interaction between contiguous colors, which lead to 
a diminution of the chromatic difference instead of its exaggeration. Fuchs, who 
discovered this effect, named it “equalization,” and it has recently been studied by 
the psychologist Musatti. Let us give an example of such an equalization effect: 
red circles which are actually identical in color everywhere seem more purple on a 
blue ground and more orange on a yellow one. Professor Kanizsa has described 
another similar effect, the action of the marginal gradient of the figure upon its 
chromatic appearance: the lack of abrupt separation between a red ground 
and yellow bars changes the red tonality of the ground which orange equalization 
makes less different from the yellow of the figure. Whether such phenomena are 
retinal or cortical or both, we do not know. Though they seem rather simple when 
we try to explain them by words as psychologists do, they disclose probably the 
astonishing complexity of synapses in the brain (or in retina—but the retina is a 
peripheral part of the brain). 


A 
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And how ean we explain color constancy, one of the major facts in the psychology 
of color vision, which allows us to identify objects in spite of the physical variations 
in the intensity and composition of the light which falls on them? It is an aspect 
of the great problem of generalization and abstraction, which was analyzed by 
Professor Hebb in terms of cortical connections in his book, The Organization of 
Behavior; by Pitts and McCulloch by formal mathematical methods; and by 
Beurle by means of a probability theory of the aggregates of interacting neurons. 
But I agree entirely with Sholl’s conclusion that “at this time it would be premature 
to attempt to put forward a complete theory of cortical organization” which would 
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explain in particular the psychological laws of color vision. I expect that in the 
future anatomists and physiologists will add greatly to our knowledge in these very 
interesting problems, but as a conclusion I think best to remind you of an old 
picture in Deseartes’s book on man (Fig. 3). Our modern representations would 
not differ much from Descartes’s, except that the “pineal gland,” which he supposed 
to be the“seat of imagination and common sense,” is now out of date and replaced 
by the cortex. As to the transmission of visual messages to the brain and storage of 
visual engrams in the cerebral substance, we have not added a great deal to this 
picture of the seventeenth century. It is sometimes good, for the modern scientist 
who is so proud of the progress of science, to learn modesty in an old French 
book which was written three hundred years before the beginning of this so- 


called “atomic age.” 


STUDIES OF COLOR BLINDNESS: A UNILATERALLY 
DICHROM ATIC SUBJECT 


By C. H. GRAHAM AND Yun Hsia 
COLUMBIA UNIVERSITY, NEW YORK 


Since the time of Thomas Young! considerations of color blindness have played 
in important and possibly primary role in determining the forms of various theories 
of color. The literature of color blindness reserves an important place for studies 
of unilaterally color-blind subjects, that is, persons who are color-blind in one eye 
but normal in the other. Judd? states that thirty-seven cases of unilateral color 
blindness have been reported in the literature, but, because of inadequate procedures, 
equipment, or background information, only about eight have proved useful for 
theory, and then not in all desirable details. No ordinary color-blind subject can 
tell us how the colors he sees compare with those seen by a normal subject, but this 


is precisely what a unilaterally color-blind subject can do. He can make a direct 
comparison of colors seen by his color-blind and normal eyes and thus provide us 
with a basis for relating the color perceptions of the color-blind eye to those of the 


normal, 

We have recently had an opportunity to examine some visual discriminations of a 
unilaterally dichromatic subject. The results of that investigation form the con- 
tent of this discussion. 

Color Mixture.—The fact that our subject is unilaterally color-blind was estab- 
lished well enough in preliminary tests, but its precise specification was unequi- 
vocally established by her data on color mixture. 

The color-mixture experiments show that our subject’s normal eye can match 
any wave length of the spectrum by a mixture of three primaries—460, 530, and 
650 mu. In addition, they give data that are in accord with Wright’s standard 
results in 10 normal subjects. The luminance units involved in the settings fall 
within the distribution given by Wright’s subjects. In a word, our subject’s 
right eye has normal color vision.* 

The data on the color-blind eye are different from the data on the normal eye. 
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In her color-blind eye our subject can match any wave length of the spectrum 
with a combination of two primaries, 460 and 650 my. Judging by this per- 
formance, we can say that her left eye is unequivocally dichromatic. Her results 
are in accord with Pitt’s standard data* on dichromats. 

The color-mixture curves for our subject’s left eye show that that eye is di- 
chromatic, but they do not tell what type of dichromasy is represented. Two 
rare types of dichromasy are excluded, but the data do not differentiate between 
two forms—deuteranopia, the most frequently encountered form (sometimes 
inaccurately referred to as “green blindness’), and protanopia, the next most fre- 
quently encountered form (sometimes inaccurately referred to as “red blind- 
ness’’). 

Spectral Sensitivity. Differentiation between the two forms may be accom- 
plished in a number of ways, all depending upon a testing of the sensitivities of the 
eye to red light. If a dichromatie eye is relatively insensitive to red light, it may 
be classified as protanopic. If the eye is sensitive to red and at the same time 
gives the same kind of color-mixture curve as our subject’s dichromatic eye, then 
it is deuteranopic. Some data by 
Hsia and Graham?’ and by Hecht and 





Hsia’ on groups of protanopic, deu- 
teranopic, and normal subjects indi- 


cate that deuteranopes usually lose Mermet bye 


sed ie lor- 
sensitivity in the blue and green re- saa ves td 


gions of the spectrum. 
The spectral sensitivity curves of 





Log Relative fovea! Sensitivity 


our unilaterally color-blind subject,’ 





shown in Figure 1, demonstrate that sal 


her dichromatic eye loses sensitivity Wavelength ~ Md 

in the blue-to-green region of the Fic. 1.—-The spectral sensitivity curves for a 
spectrum and has normal sensitivity unilaterally color-blind subject. The upper 
: / < : ie curve for the normal eye; the lower one, the 
in the red. She is therefore classifi- curve for the color-blind eye. Data of Graham 
able as deuteranopic. In addition, and Hsia (these PROCEEDINGS, 44, 46-49, 1958). 
her data on spectral sensitivity 

strengthen the interpretation that deuteranopes usually show a loss of sensitivity 
in the blue-to-green region of the spectrum. Her loss in this region is well marked. 


Binocular Color Matches.-As part of our analysis of our subject’s vision we 


performed an experiment on binocular color matching. This experiment had to 
do with the question: What colors does our subject see in her color-blind eye? 
The subject regarded the color of a vertical slit presented to her left eye in the 
visual field and simultaneously observed the color given by a horizontal slit in her 
normaleye. The wave length of light stimulating the normal eye could be changed 
until a narrow wave-length band was obtained that, in the opinion of the subject, 


matched the color seen in the dichromatie eye. 

The results of the experiment are summarized in Figure 2. In general, it seems 
that in her dichromatic eye the subject matches all wave lengths greater than 
the neutral point (.e., the narrow wave-length band that a dichromat sees as 
gray) with a wave length seen in the normal eye lying near 570 mu. Wave lengths 
shorter than the neutral point in the dichromatie eye are matched in the normal 
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eye by a blue of about 470 my. This result is in line with Judd’s analysis? of earlier 
data. 

Discussion.—Our subject shows a loss of spectral sensitivity in her color-blind 
eye and at the same time sees all wave lengths above 502 my as yellow. How ean 
she (or, in fact, any deuteranope) see yellow (which we take to be a mixture of red 
and green) in the total or even partial absence of a green brightness mechanism? 

To begin with, it seems that, in order to account for the seeing of yellow (we 
shall not deal with blue!) by deuteranopes, we may appeal to the Fick* and Leber® 
theory. This theory says that yellow is seen by deuteranopes because, in such 
persons, the red and green receptors become indistinguishable so far as sensory 
input is concerned. However, they still maintain their usual central connections, 
and thus, on stimulation by wave lengths longer than about 500 my, they provide 
a mixture of yellow due to excitation of the usual red and green centers. This 
latter consideration sets a requirement. (In discussing it, we shall refer to the 
normal red and green processes as /? and G and their possible transformations, sen- 
sory or central, as R’ and G’.) (1) It is required that spectral luminances for a 
mixture of yellow determined by R’ and G’ must be in nearly the same ratio 
approximately unity—at all stimulating wave lengths. (2) In addition, a second 
requirement must be met. For the case of the deuteranope (or, for that matter, 
the protanope) the sum of the R’ and G’ ordinates at any wave length must add 
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to give deuteranopic (or protanopic) luminosity, as specified, for example, by the 
data of our groups of dichromats. 

These relations are shown in Figure 3. The full curve shows color mixture 
(minus the contribution of a blue fundamental) in the case of a deuteranope, com- 
parable to one of our subjects, who shows no luminosity loss. (The curve is also 
nearly the luminosity curve for normal subjects.) The right dashed curve repre- 
sents average deuteranopic luminosity with appropriate loss. In the latter curve, 
the R’ and G’ ordinates are taken to be in a ratio of unity, as dictated by color- 
mixture data. The area under the curve is assumed to be proportional to the 
number of R’ and G’ cones (both being identical in all details except their central 
connections). Each type contributes half the total number. Similar considera- 
tions could apply to the logical extension of Fick-Leber ideas to protanopia, as 
represented in the left dashed-line curve. 

The results of the present experiments concerning luminosity losses in dichro- 
mats may be at least qualitatively in line with results described by Rushton in 
another paper of this symposium. For the moment we may regard the main 


merit of the present accounts as stating how trichromatic theory may account 


for luminosity losses and the colors seen by dichromats."” 
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discussion of this research see C. H. Graham and Yun Hsia, Science, 127, 675-682, 1958. 
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COLOR VISION: THE APPROACH THROUGH INCREMENT- 
THRESHOLD SENSITIVITY 
By W.S. STILES 
NATIONAL PHYSICAL LABORATORY, TEDDINGTON, MIDDLESEX, ENGLAND 


I am most sensible of the honor of addressing this meeting and of the generosity 
of the National Academy of Sciences in bringing me here as their guest. The in- 
vitation to contribute to the present symposium came just as I had finished a rather 
elaborate study of the color-matching properties of the average observer. How- 
ever, that work I believe to be of more colorimetric interest, and I have elected, 
therefore, to deal with a quite different approach to color vision by a method which 
does not rest on actual color-matching—the method of increment thresholds. 

Consider, to take first the simplest case, a subject whose eyes are fully adapted 
to a completely dark field. The increment threshold is then just the smallest 
perceptible quantity of a certain test light which is applied at a particular part of 
the visual field. For convenience, I use the expression “smallest perceptible 
quantity,”’ although, of course, thresholds are definable only statistically, for ex- 
ample as critical values of the stimulus variable at which the chance of a perception 
equals 50 per cent. In the general case, the increment threshold is still the smallest 
perceptible quantity of a certain test light, but now this is added to another stimu- 
lus—the conditioning stimulus—which consists of any distribution of brightness 
and color to which we may choose to expose the eye. The conditioning stimulus 
need not be an unvarying one: the increment threshold is, in fact, often measured 
while the eye is adapting itself to a change from one constant stimulus to another 
constant stimulus, the whole spatial and temporal pattern of stimulation being 
considered to form the conditioning stimulus. The increment threshold (or, 
better, its reciprocal, the threshold sensitivity) is a measure of visual sensitivity 
which can be varied in different ways to bring out particular properties of the visual 
response system. We may, for example, vary the angular size and exposure time 
of the test stimulus, to determine how well the visual system integrates over area 
and over time. Most important, we may change the spectral composition of the 
test stimulus, in order to determine the spectral sensitivities that are active in 


the response system. So far, only changes in the parameters of the test stimulus 


have been mentioned. In addition, any number of different conditioning stimuli 
may be used. However, there is evidence that the modifications of the increment 
threshold by more complex patterns of brightness and color are similar to those 
obtained with a very simple pattern, namely, a uniform stimulus which extends 
over and around the area of the field where the test stimulus is applied. We should 
regard the action of a more stimulus, including stimuli applied at some distance 
from the test area, as equivalent to that of a suitably chosen uniform stimulus of the 
kind mentioned if, for all possible test stimuli (of various sizes, colors, exposure 
times, ete.), the observed increment threshold was the same for the two condi- 
tions. The equivalence is probably not generally complete, but it has been found 
to hold approximately in some cases, and we are justified in concentrating, in the 
first instance, on the simple stimulus pattern. This has the advantage that much 


100 
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bigger effects on the increment threshold can be obtained with ‘covering fields’’ 
than with conditioning stimuli confined to the surrounding retina. 

For probing the color mechanisms of the eye, the method that has been ex- 
tensively used is to apply, in flashes, a small monochromatic test stimulus of wave 
length \ to a uniform monochromatic conditioning field of wave length uw and in- 
tensity M, extending over the test area (the two-color threshold technique). Most 
work has been done for the case when the eye is fully adapted to the conditioning 
field, although several recent investigations (particularly by Boynton!) as well 
as some earlier ones have considered non-equilibrium situations. If all the other 
possible experimental variables (size and exposure time of test stimulus, angles of 
incidence on the retina of test, and conditioning stimuli, ete.) are kept constant, 
the increment threshold Ny, i.e., the smallest perceptible intensity of the test 
stimulus, will be a function f(A, w, 7) of the three variables A, u, W,. The fune- 
tion f is certainly not simple, and its determination for the whole range of possible 
values of the three variables is a considerable undertaking. Some progress has 
been made, however, which suggests that the function may be put in a form in- 
volving separate terms identifiable, on certain reasonable assumptions, with the 
contributions to the increment threshold of component visual mechanisms having 
different spectral sensitivities and some other differences in their properties. The 


” 


possibility of ‘‘separating”’ the function f in this way does not, of course, prove 
that there exist in the eye separate structures and processes corresponding to the 
different terms. The psychophysical analysis has to be related to objective ob- 
servations on the visual pigments, on the types of end-organ and the ramifications 
of the nerve pathways through the higher neurons, on the electrophysiological 
response, etc., before the “component mechanisms” can begin to acquire a real 
existence. Despite the wealth of information on the visual system generally that 
is provided by these objective techniques, there has until recently been little which 
could be applied directly to the psychophysical data of color vision. The impor- 
tant new development is the objective detection and measurement of bleachable 
visual pigments, other than rhodopsin, in the living human eye (Rushton). Some 
repercussions of this work on the increment-threshold analysis will be mentioned 
later. 

The guiding principle in interpreting increment thresholds is that in the retinal 
area on which the test stimulus is imaged there are several component mechanisms, 


‘seen”’ when the 


any one of which may be the means of triggering off the response 


test flash is applied. Each component mechanism is conceived of as an association 
of a selected proportion of the end-organs (rods and cones) in the area, the 
association being effected somewhere in the neural system through which the 
nervous activity, initiated by light absorption in the end-organs, is transmitted 
to the brain and finally to the effector nerves of the external response. It is not 
assumed that the end-organs of a component mechanism are necessarily all of the 
same kind or contain just one photopigment, nor that end-organs belong exclu- 
sively to one component mechanism, but the mechanisms are assumed to have 
sufficient independence for the idea of a particular component mechanism acting 
alone to have a meaning. It is assumed, in fact, that the response of all but one 


component mechanism could be blocked (ideally; it may not be practicable to do 
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it), leaving the remaining mechanism to respond in its normal way. For a given 
test stimulus and given conditioning stimulus, the observed increment threshold 
in the retinal area is assumed to depend only on the respective increment thresholds 
of the component mechanisms that would be observed under the same conditions 
if they were acting alone. In particular, it is assumed that the observed increment 
threshold will not exceed any of the component increment thresholds and _ will 
generally approximate to the smallest of these. Finally, each component mecha- 
nism is assumed to have a characteristic relative spectral sensitivity function, so that 
the effect on the mechanism of an intensity /; (in energy units) of one wave length 
\,, forming the whole or part of the test or conditioning stimulus, is just the same 
as that of an intensity /. of any other wave length d» that may be substituted for it, 


provided that /; and /, are in the inverse ratio of the relative spectral sensitivity 


values at these wave lengths. 


5° PARAFOVEA 
A= 58Omp 
kh = SOOmip 


LOG Mu 
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! | 
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Fic. 1.-Inerement threshold, Ny, versus field intensity, Mp 
for 5°’ extrafoveal vision. Test stimulus: 1°, 0.06 sec. Subject: 
W.S.S.) (Va and My expressed in erg sec~! [deg. of are] ~?.) 


The analysis of increment-threshold measurements on the above lines may be 
illustrated by examining briefly the case of extrafoveal vision, for which there is 
the ample evidence of the duplicity theory of vision for the operation, in some 
sense, of distinct rod and cone mechanisms. Figure 1 gives an example of a so- 
called threshold-versus-intensity (t.v.i.) curve, in which the logarithm of the inere- 
ment threshold is plotted against the logarithm of the uniform conditioning stimu- 
lus for a particular pair of test and conditioning wave lengths, A = 580 my, u 
500 mu. In this case the test stimulus (1°, 0.06 sec.) was imaged 5° from the 
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Fic. 2.— Diagram illustrating how t.v.i. curves like that of Fig. 1 result from 
the t.v.i. curves of two component mechanisms, A (rod) and B (cone), whose 
positions with respect to the axes of log Ny and log M, for any test and field 
wave lengths A and w are determined by the respective spectral sensitivity 
curves of the mechanisms shown in the auxiliary diagrams 


fovea. Figure 2 shows diagrammatically the tentative interpretation of the curve 
of Figure | in terms of two component mechanisms of the kind contemplated above. * 
Curves A and B (Fig. 2) represent, respectively, the t.v.i. curves assumed for the 
A (rod) and B (cone) component mechanisms when acting alone. At each field 
intensity the observed increment threshold is put equal to the lower of the values 


appropriate to the A and B curves. The crucial question now is What happens 


if \ and w are changed? According to the assumption about the spectral sensi- 
tivity curve of a component mechanism, a change in A, keeping » constant, should 
merely displace the curves A and B—without change of shape—parallel to the 
ordinate axis, by amounts determined by their respective relative spectral sensi- 
tivities. A change in yp, keeping A constant, should, on the other hand, displace 
the curves parallel to the axis of abscissae by amounts also fixed by the respective 
spectral sensitivity curves. The spectral sensitivity curves are also shown in the 
diagram, and it is easy to visualize the effects of changing \ and uw. The experi- 
mental t.v.i. curves obtained for various combinations of test and field wave lengths 
show that the branch of the curve attributed to the A mechanism conforms satis- 
factorily to these displacement rules, the displacements parallel to the two axes 
being controlled by the same spectral sensitivity curve. For the B branch, how- 
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ever, the displacement rules fail to hold, and it is impossible to represent the incre- 
ment threshold for extrafoveal vision as arising from just two component mecha- 
nisms in the way proposed. The next step would be to see whether the B branch 
could be represented as the resultant of two or more branches, but this is difficult 
because the form of the B branch for many of the interesting (A, ~) combinations 
is not accessible to observation: it is masked by the A branch, which has the lower 
threshold. 

We now examine the various component mechanisms which we are led to postu- 
late for the fovea, the area of the retina where color discrimination is most highly 
developed and where the absence of rod end-organs eliminates (probably) any 


participation of a “rod”? component mechanism. (This last simplification is not 


quite a certain a priori, as the adapting field extends into the retinal area containing 
rods, but it seems to be so in fact.) Consider, first, the t.v.i. curves obtained with 
test stimuli of short wave length (A < 510 my) and with conditioning stimuli of 
medium or long wave length (u > 520 my). These show a division into two branches 
(u < about 560 my) or of three branches (u > 560 my), an example of the latter being 
given in Figure3. The three 
branches are provisionally 


ascribed to three component 
FOVEA 
A= 476 mb 
=578 mh neutral symbols a4, m4, 3 


marked in the figure. The 


mechanisms, denoted by the 


changes in the t.v.i. curve 
when \ and yw are altered 
establish that the branch 7, 
obeys the displacement rules. 





But if A becomes too large 


ia) — 
ZERO i or uw too small, the a, branch 

LOG Mu 
Fic. 3.—Foveal t.v.i. curve. Test stimulus: 1°, 0.2 : 
sec. Subject: EK. K. (Ny and My, in erg see~! [deg. of thus the displacements ob- 


are|~? servable by changing \ and yu 


is masked by the as branch; 


correspond to non-overlap- 
ping spectral ranges of the spectral sensitivity curve of the mechanism, and it is 
not possible immediately to check that a single spectral sensitivity curve 
determines displacement of both kinds. The difficulty can be largely over- 
come by the “auxiliary field method.’”’ In this we use as conditioning 
stimulus a mixture of a fixed auxiliary field of long wave length sufficiently 
intense to bring the increment threshold to the a, branch of the curve, and a main 
field of wave length uw, which may be as short as desired. The comparative effects 
of main fields of different wave length in moving the increment threshold higher up 
the m branch can then be determined.‘ It is found that in the overlapping region 
(A < 500 my) the \ and yu displacements correspond to a common relative spectral 
sensitivity curve and that this curve has a maximum at approximately 440 mu 
(energy expressed in quantum units). The evidence on the a; branch is more 
limited because of the difficulty in obtaining the very high field intensities needed 
to reach it. The results again indicate a branch obeying approximately the dis- 
placement rules and with consistent \ and yw displacements in the overlap region 
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(A < 500 mu). As can be seen in Figure 4, the spectral sensitivities of the a, and 
mr; component mechanisms (derived from u-variation displacements) agree closely 
in shape in the region of their maxima (A < 500 my) and differ only at long wave 
lengths. (The curves of Figure 4 and later figures refer to the mean results of the 
same four subjects.) But for test stimuli in the short-wave region the two mecha- 
nisms still differ in their absolute thresholds, i.e., the values to which their inere- 
ment thresholds tend as the conditioning stimulus is reduced to zero: the absolute 
threshold for 73 is about four times that for a;. It was unexpected to find two com- 
ponent mechanisms with maximal 


aa fs WAVE LENGTH 
and closely similar spectral sensitiv 100 my 150 500-550 600 650 700 


ities in the blue. But there is 
also some evidence of a third com 
ponent mechanism (a:) with maxi 
mal sensitivity in the same spectral 
region. This is indicated by a 
breakdown in the displacement 
rules for the branch a, when, for 
conditioning fields of long wave 
length (e.g., 600 my), the wave 
length of the test stimulus is re 
duced below about 460 mu. The 


FIELD SENSITIVITY 


discrepancy would be consistent 
with a m2 mechanism with the spec- 


LOG 


tral sensitivity curve (by the u-var 
iation method) inthe range 510-700 
mu, as shown in igure 4. For u - 

510 mu, the curve must lie above 
the points of the arrows shown in 
the figure. Unfortunately, the 


25, 000 20,000 5,000 em 
auxiliary field device is not appli- FIELD WAVE NUMBER 1/u 
cable in this ease, and the form of Fic. 4 Field spectral sensitivity of foveal cone 
mechanisms 7m, zm, and 3. Field sensitivity is 
. : : : defined as the reciprocal of the field intensity 
in the neighborhood of the maxi- in quanta see~! (deg. of are)~? required to raise 
the increment threshold of the mechanism to ten 
2 ‘ times its zero-field value. For \ < 510 mg, the 
lie outside the range 440-480 mu. arrow points indicate values below which the 
As the wave length A» of the undetermined part of the am. curve cannot fall 
A A ~ Means for our subjects, r, and 73; 1Ol three sub- 
test stimulus is increased above — jects, x 


the curve has not been determined 


mum, although the latter can hardly 


about 450 my, the branches of 

the t.v.i. curve associated with mechanisms having maximal sensitivity in the 
blue move upward with respect to the ay branch, and for A greater than 520 
mu they move completely above it and are inaccessible to measurement, whatever 
the wave length of the conditioning stimulus. For these test stimuli of longer 
wave length, the a, branch does not obey the displacement rules, although the well- 
marked transitions obtained at shorter wave lengths are not in evidence. It was 
originally thought that the results would be explained by admitting that there were, 
in fact, two branches—-a ay and a 75 branch —with associated spectral sensitivities 


not very widely separated along the spectrum, so that sharp transitions from one 
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branch to the other would not be expected. Further measurements indicated, 
however, that, while this explanation is fairly satisfactory if attention is confined 
to conditioning stimuli below about 150 trolands, it leads to discrepancies with 
observations made at higher intensities. The spectral sensitivities of 2, and 7; de- 
rived by the w-variation method from measurements in the lower-intensity range are 
shown in Figure 5. These sensitivities, like those of Figure 4, represent at each 
wave length the reciprocal of the field intensity (in quantum units) required to 
raise the increment threshold of the particular mechanism to ten times its zero-field 
value. If we also know the shape of the branch of the t.v.i. curve appropriate to 
the mechanism and the absolute threshold of the mechanism for some test wave 
length, we can derive the curve showing the variation of the increment threshold 
or its reciprocal, the threshold 
sensitivity, with test wave 
length for any intensity and 
wave length of the field. This 
calculation has been made in 
Figure 6 for the five compo- 
nent mechanisms 2,~—7; for the 


case of a field of wave length 
667 muand intensity 29 tro- 
lands. The observed thresh- 


old sensitivity at each test 
wave length should equal ap- 


mt 

- 
t 

vA 
L 


proximately the highest of the 
component threshold sensitivi- 


LOG 


ties, and it is apparent that 
only wm; and m5 are concerned 
in this instance. The circle 
points represent the directly 
| observed threshold sensitivities 


a ee: ee . nO 
22,500 20,000 17,500 15,000 em. ~! and are in satisfactory agree- 
FIELD WAVE NUMBER 1/uz . + Le ry 

ae decay ieee ment with the prediction. The 
Fig. 5 Field spectral sensitivities as in Fig. 4 for cone 
mechanisms mw, and 7 The inset diagram shows the 
log (ratio of the sensitivities: high-intensity level/low- and observed curves when the 
intensity level) for the mw, and 2; mechanisms. Means for 


5,000 


comparison of the predicted 


four subjects field intensity Is raised bya fac- 
tor of nearly 100 is shown in 
Migure 7. The discrepancy is not only in the absolute values (which might arise 
from uncertainty in the shape of the zy branch) but also in the shape of the observed 
curve in the wave-length range where it should be at least parallel to a. The 
indication is that, with a high-intensity red field, the a, mechanism is relatively 
more sensitive on the short-wave side of about 580 mu than the calculation treating 
T, as a single component mechanism with its low-intensity properties would pre- 
dict. This deviation is independently confirmed by t.v.i. curves taken to high 
intensities. An analogous, although smaller, effect for a; is shown by observations 
with low- and high-intensity green conditioning fields (igs. 8 and 9) and is again 
confirmed independently by t.v.i. curves. The a; mechanism under high green 
stimulation is relatively more sensitive on the long-wave side of about 580 my than 
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the calculation treating a5 as a single component mechanism with its low-intensity 
properties would predict. Insufficient is yet known of these deviations from single 
component properties of a, and m5 to say whether they could be explained by as- 
suming additional component mechanisms. For the present, it is advisable merely 
to recognize modified high-intensity conditions of these mechanisms, denoted by 
my’ and m5’, respectively. The differences in spectral sensitivity in the low- and high- 
intensity conditions of mor m5 are not large compared with the differences between 
mand m5. A first determination of the ratio of the high- to the low-intensity spee- 
tral sensitivity is shown in the inset diagram of Figure 5. A tentative deduction 
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(29 


15,000 em. 
ROLANDS) 
Fig. 6.—Points: Inerement-threshold sen- 1 Nok 
sitivity by 1/.V\ versus test stimulus wave 20.000 17.500 15,000 em 
number 15,000 em~! (wave length 667 mu 500 571 667 my 
and intensity 29 trolands. Test stimulus: 1/d 
0.2 see ans ur subjects. Curves 1/p 15,000 em 
.. . ec. Means of four ubjec irve 3500 TROLANDS) 
predicted threshold sensitivity versus wave- 
number curves at this conditioning stimulus Fic. 7 As for Fig. 6, but for conditioning 


for the component mechanisms m—7 field of intensity 2,500 trolands. 


from other observations is that the transitions 2 a 74 aud a3s;—> mr, occur ina range 
of conditioning intensities of the order of 1 log unit centered on an intensity of about 
300 trolands. 

The scheme of component visual mechanisms indicated by increment-threshold 
analysis is summarized in Table 1. 9 The Fechner fraction quoted in the final 
column is defined as the ratio of the increment threshold to the field intensity when 
the conditioning field has the same color as the test stimulus (A uw) and the inere- 
ment threshold is raised by it to ten times the zero-field value. There are inter- 


esting differences between the Fechner fractions of different component mecha- 
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| 
| 
20.000 17,500 15,000 em . J — 
1/d 20,000 17,500 15,000 em 
1/u 20,000 em 1/r 
3 TROLANDS 1 20.000 em 
Fig. 8 As for Fig. 6 but for conditioning 1680 TROLANDS) 


field of wave number 20,000 cm =! (wave hia. 9.—As for Fig. 8 but for conditioning 
length 500 my) and intensity 43 trolands field of intensity 1,680 trolands. 

nisms: the “blue”? cone mechanisms have values four to five times those for the 
“green”? as “red”? cone mechanisms. In applying the present ideas of component 
mechanisms to the explanation of other psychophysical measurements on color 
discrimination and brightness matching, this difference plays an important part: 


it corresponds to the small contribution to luminance of short-wave stimuli, despite 


their strong coloring value. 


TABLE 1 
Wave Length 
of Maximal Fechner Fraction 
Sensitivity 1°, 0.2 Sec 
Mechanisn \ Remarks (my) Test Stimulus) 
tod T Absent at the fovea 503 (30) 
Ratio of blue to approx. 10 8.7 
vellow ' 2.6 log 
sensitivit\ ' units 
(440/590 
T atio of blue order of f f 
Blue’ come vellow | log (Between 440 (Probably as for 
sensitivity unit and 480 my 7, and x 
(440/590 
atio of blue to approx, 440 8.7 
vellow : 1.0 log ; 
sensitivity units 
(440/590) 


» 9 


2 


T 4 540 
“Green” cone ‘ar; Modified high-intensity 540 
stute 
1 575 
‘Red”’ cone ‘ (very flat 
f Modified high-intensity 587 
state 
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The discussion so far has been concerned with the breakdown of the three-vari- 
able function f(A, uw, W,), all the parameters of the test stimulus other than color 
being assumed constant. If one of these is allowed to vary, the effects on the 
different branches of the t.v.i. curves can be compared, and more information on 
the properties of the component mechanisms can be obtained. For example, it is 
known that the response of the cone end-organs of the retina is highly directional, 
and, by determining the change in increment threshold on altering the angle of 
incidence on the retina of the rays forming the test stimulus (keeping, of course, 
the retinal area stimulated the same), differences in the directional properties of 
different mechanisms can be studied. The differential effects are small, and more 
measurements are needed, but the early work indicates that z, is more directional 
in its reponse than zy, that a4 is more directional than a; for stimuli in the green and 
yellow, but that this position is reversed for red stimuli. Changes in the relative 
response of different mechanisms as the angle of incidence is varied may be ex- 
pected to produce changes of color even of monochromatic stimuli. The small 
changes, observed are, in fact, in general agreement with the relative directional 
properties of the mechanisms obtained by the threshold observations. 

Two outstanding questions are: (a) What is the relation of the increment- 
threshold mechanisms to the three processes of trichromatic color-matching? and 
(b) Do the spectral sensitivities of the component mechanisms represent photo- 
sensitivities of visual pigments? We limit ourselves to foveal vision, so that the 
complication of the rod mechanism does not arise. The two questions are closely 
related, as the three small-field color-matching functions which sum up the color- 
matching properties of the foveal retina (conditions of very high intensity excluded) 
are almost certainly linear combinations of the spectral photosensitivities of three 
visual pigments, although the coefficients in these combinations are not determin- 
able by straightforward color-matching measurements. Thus, if the spectral 
sensitivity of a component mechanism could be represented as a linear combination 
of the color-matching functions, it would also be a linear combination of the pig- 
ment photosensitivities. This would indicate that the end-organs of the mecha- 
nism contained the pigments represented in the linear combination (i.e., those hav- 
ing non-zero coefficients) and that the effects of light absorption in these different 
pigments added up linearly. Electrophysiological evidence makes it unlikely that 
this last state of affairs could come about unless all end-organs of the mechanism 
contained the pigments in a mixture of the same composition, including the spe- 


cially simple case of a “mixture”? comprising just one pigment. The rigorous 


application of the argument is hampered by the difficulty of obtaining “hard” 


data from increment-threshold measurements. However, tests made on the mean 
spectral sensitivity curves (ligs. 4, 5) indicate that a, or a,’ can be expressed fairly 
well as linear combinations of the small-field color-matching functions. For a5 the 
linear representation is less satisfactory, particularly for the high-intensity condi- 
tion a5’. For a; and for m in the main (short-wave) region of the curve, a linear 
representation is acceptable, and even with the inclusion of the long-wave lobe on 
the m, sensitivity curve it could not certainly be excluded. The actual shapes of 
the sensitivity curves also have a bearing here. The spectral absorption and spec- 
tral photosensitivity curves that have now been obtained for many visual pigments 
from animal retinas, by measurements on retinal extracts or complete excised 





110 PHYSICS: W.S. STILES Proc, N. A. S. 


retinas, differ little in shape from the absorption curve for the pigment of rod 
vision, rhodopsin; to a first approximation, they correspond to the rhodopsin curve 
(log absorption factor against wave number) displaced along the wave-number axis 
so that its maximum occurs at the appropriate position. Before physical absorp- 
tion and photosensitivity curves can be compared with psychophysical sensitivity 
curves, the latter must be corrected for the selective light losses in the eye that 
precede the actual absorption by the photosensitive pigment. The correction is 
uncertain in the violet but is prob- 

WAVE LENGTH u 


100 mu 150 500 550 600 650 700 


ably not seriously in error in the 
blue and green: it is negligible at 
longer wave lengths. It is found 
that the corrected a, 24’, and ms 
spectral sensitivity curves and the 
7, Curve ignoring the long-wave lobe 
are fairly well represented by dis- 
placed rhodopsin absorption curves 
(Fig. 10: the broken lines are, in 
fact, the corrected spectral sensiti- 
vity curves of the rod mechanism, 
but the latter has substantially the 
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mechanisms T1, %2, 4, 7 corrected for pre-receptor 

light losses and compared in each case with the the spectral sensitivity. of the 

similarly corrected sensitivity curve of rod vision 7 

displaced parallel to the two axes to give the best fit. 


along the spectrum) will represent 


“red”? cone mechanism in low- 
or high-intensity conditions. — It 
seems certain that the end-organs of a; cannot be equipped either with a single 
rhodopsin-type pigment or with a simple mixture (the same in all end-organs) 
of two such pigments. It has already been noted that the 25; and a5’ spectral 
sensitivities are not adequately represented by linear combinations of the small- 
field color-matching functions, and the failure to fit rhodopsin-type curves might 
result from the fact that the increment-threshold spectral sensitivity is not 
built up linearly from the spectral absorption curves of the pigments contained in 
the end-organs of the a; mechanism. It is most improbable, however, that this 
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is the sole cause. Other proposed fundamental (pigment) spectral sensitivities for 
color vision obtained as strict linear combinations of the color-making functions 
for which the coefficients have been derived from the properties of the three kinds 
of dichromatie color defectives and other psychophysical measurements demand a 
similar non-rhodopsin-type “red”? curve, usually including still more striking devi- 
ations from the rhodopsin shape at short wave lengths (lig. 11). 

Rushton’s? objective measurements of bleachable pigments in the human fovea 
have established the presence in protanopes, deuteranopes, and color normals of 
a pigment he has named chlorolabe with a difference spectrum and action spectrum 
having a maximum at 540 my. As he has pointed out, the spectral sensitivity of 
the mw; mechanism is a fair approximation to the spectrum of chlorolabe. The 
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Fic. 11.--Thomson-Wright colour fundamentals corrected for 
light losses and compared with the displaced rod curve as in Fig. 
10. 


r,’ spectral sensitivity curve lies still closer to the objective curve. In fact, we 
have here a convergence, onto a common curve of rhodopsin type, of four spectral 
sensitivities: the objective spectrum of chlorolabe, the ‘“‘green’’ fundamental 
derivable from the small-field color-matching functions of normal and dichromatic 
subjects,> the increment-threshold sensitivity m4’, and, finally, the sensitivity of a 
normal eye in a state of “artificial monochromasy” produced by exposure to very 
intense violet and red stimuli (Brindley®). The implication for the “green” cone 
mechanism is that, in the high-intensity condition, its end-organs contain, or at 
least are activated by the light absorption of, just one pigment—-chlorolabe. Why, 
then, should the spectral sensitivity a, attributed to the mechanism in the low- 


intensity condition, be different? If exactly the same end-organs were concerned 
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at all intensity levels, either these would have to contain—probably in relatively 
small amount-—a second pigment more red-sensitive than chlorolabe, or, alter- 
natively, the effective spectral absorption curve of the latter pigment in its situation 
within the end-organs would have to vary with the intensity of stimulation. A 
variation of the kind required would occur if at low levels of stimulation the optical 
density of the pigment in the end-organ was high enough for self-screening to occur 
and if the reduction of density by bleaching at high levels of stimulation substan- 
tially reduced the self-screening. “‘Self-screening’’ denotes here the selective 
filtering of the light reaching the deeper pigment layers in the end-organ, in its 
passage through the upper pigment layers. Brindley has considered the possibility 
that the breakdown of metameric color-matches after the eye has been exposed to 
high intensities may be explained by the reduction of self-screening present at low 
and moderate intensities.6 The difference between the a, and 2,’ spectral sensi- 
tivities would correspond—very approximately—to a change in self-screening if 
in the low-intensity condition the chlorolabe in the end-organs had an optical density 
at the maximum wave length 540 my of about 1.0 and if this was reduced to about 
0.2 in the high-intensity condition. However, the existence of self-screening in the 
end-organs of the green cone mechanism is contrary to all the evidence: Brindley’s 
results show that the breakdown in metameric color-matches in the green-yellow- 
red range is associated with an anomaly not of the ‘“green’”’ process but of the “red”’ 
process only; the transition from a, to a4’ occurs in an intensity range lower than 
that required to produce appreciable breakdown in color-matches; the requisite 
density of the chlorolabe under zero-field conditions and the necessary degree of 
bleaching in passing from the a, to the a,’ intensity level are both improbably 
high, aecording to Rushton’s objective data on chlorolabe. As a more red-sensitive 
pigment, erythrolabe, has been shown objectively to be present in the retina (Rush- 


ton),? the suggestion that the end-organs of the “green’’ cone mechanism might con- 


tain a small proportion of it mixed with chlorolabe is, at first sight, promising. ‘The 


difficulty here is to understand why the red-sensitive pigment should be bleached 
or rendered ineffective by an increase in intensity corresponding to the transition 
my tom’. The objective data on erythrolabe would poind to a comparatively small 
reduction—not more than 10 per cent—in its density by bleaching in this intensity 
range. We are led, therefore, to the other main possibility that, at low intensities, 
additional end-organs, not containing only chlorolabe, have contributed to the 
spectral sensitivity derived for the ‘‘green’’ cone mechanism. Many ways can be 
imagined in which this might come about, but it would be pointless to enlarge on 
them, as there is at present little evidence to distinguish different suggestions. One 
general observation may be made: whatever the way in which the non-chlorolabe 
end-organs are associated with the chlorolabe end-organs in determining 7,4, the 
association is broken as a result of increasing the level of stimulation. 

It is possible that in three blue-sensitive mechanisms—ze, 7, 73—we have some- 
thing like the same principle in operation. . If we accept a; as the “blue” cone 
mechanism with end-organs containing just one blue-sensitive pigment of rhodopsin 
type (cyanolabe? not yet detected objectively), the spectral sensitivities of the 
other mechanisms would correspond to associations of cyanolabe with nou-cyano- 
lube end-organs, these associations being broken at higher levels of stimulation by 
long-wave light. On this view, only one of the three mechanisms would be in 
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being at any moment, and this represents formally a radical change in our working 
hypothesis. But, while it can be shown from the steady-state increment-threshold 
measurements that the mechanisms m, and a4 are present together, the relation- 
ship of the spectral sensitivity curves of 7, m2, and 3 prevents a similar proof that 
any two of these are simultaneously present. 

Comments generally similar to those made on a4 and 2,’ apply to the change 
from 75 to m5’, but there are some differences. There is an additional objection to 
the use of self-screening as an explanation of the change; bleaching of the self- 
screening pigment would lead to a change in the spectral sensitivity curve even 
qualitatively different from that observed. This is not necessarily in conflict with 
Brindley’s conclusion that, at higher-intensity levels, a change in self-screening of 
the end-sensitive pigment accounts for the breakdown of metameric color-matches, 
although the objective data on the erythrolabe in the retina do not support this 
theory. There is no convergence of the different spectral sensitivities with maxima 
at long wave length on a common rhodopsin-type curve. The 75 curve does not 
agree with the objective spectrum of erythrolabe: Rushton has shown that 75 can 
be fitted by the absorption curve of a mixture of erythrolabe and chlorolabe in the 
proportion 1:0.6.2 For a5’, the corresponding proportion would be about 1:0.3, 
so that a,’ is nearer the objective curve of erythrolabe. The fact that both 
erythrolabe and chlorolabe have been shown to be present in the deuteranopic eye 
favors the Fick hypothesis of red-green fusion for the explanation of deuteranopia. 
The character of dichromatic color vision as a reduction form of normal color vision 
means that the spectral sensitivity of the red-green fusion in the deuteranope must 
be a linear combination of the normal’s color-matching functions and hence almost 
certainly of the spectral absorption (strictly photosensitivity) functions of the 
normal’s pigments. Deuteranopia would then be explained by the presence in all 
the end-organs of the fused red-green system of a mixture in an invariable ratio of 
the chlorolabe and erythrolabe pigments; any other explanation would encounter 
linearity difficulties. If, as is sometimes assumed, the normal also has the fusion 


process of the deuteranope but, in addition, the simple ‘green’? process (and, of 


course, the ‘blue’? process also possessed by the deuteranope), then the appearance 
in the normal eye of a a5-75’ mechanism supplied by end-organs containing a mix- 
ture of chlorolabe and erythrolabe would fit in very well. But there are difficulties: 
on the one hand, the pigment mixture appears unexpectedly stable when very high- 
intensity colored stimuli are used in increment-threshold and artificial mono- 
chromasy measurements, and, on the other hand, the breakdown of color-matches 
at high-intensity levels could not correspond merely to an alteration of the com- 
position of a pigment mixture in the end-organs of the ‘‘red”’ cone mechanism. 
This review has been devoted mainly to the spectral sensitivities of the color- 
vision mechanisms of the fovea arrived at by increment-threshold measurements 
made under steady-state conditions. The spectral sensitivities are perhaps of 
special interest at this time in the light of the new objective data on cone pigments. 
The summation properties of different color mechanisms, the way the mechanisms 
are modified in moving out from the fovea, and their behavior when the condition- 
ing stimulus is suddenly changed have all been studied by the increment-threshold 
method to a limited extent. There is an interesting field here for further work. 
1R.M. Boynton, J. Opt. Soc. Amer., 46, 172, 1956. 
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VISUAL PIGMENTS IN THE INTACT HUMAN EVE 
By W. A. H. Rusuron 
CAMBRIDGE UNIVERSITY, CAMBRIDGE, ENGLAND 


It. is generally agreed that the first step in the visual process is the absorption of 
light by one or more photosensitive visual pigments contained in the rods and cones 
of the retina. The pigment of the rods has long been recognized and studied under 
the name rhodopsin, but rod vision is the colorless vision of twilight, and thus, to 
understand the photochemical basis of color vision, we need to measure the visual 
pigments in the cones. No extract of mammalian eyes has ever yielded a visual 
pigment other than rhodopsin; this is not surprising in view of the observations of 
Wald, Brown, and Smith! that the chick’s eye, where cones relative to rods are a 
thousand times as abundant as in man, still contains twice as much extractable 
rhodopsin as iodopsin (the cone pigment). If in man we are to measure the pig- 
ments in cones, we must take advantage of two exquisite features of retinal organ- 
ization. First, there is a tiny area of the retina—the fovea centralis—0.6 mm. in 
diameter, where there are no rods and where the cones are close-packed. If meas- 
urements can be confined to this region, contamination by the enormous preponder- 
ance of rhodopsin will be absolutely excluded. 

The difficulties in dealing with so small a sample are somewhat offset by the sec- 
ond feature of cone structure, for cones are beautifully organized to catch quanta 
and secure photolysis. By measuring the pigments zn sztu within the cone, we take 
advantage of this organization and not only detect physically the change in pig- 
ment density during bleaching and regeneration but measure this change to an 
accuracy of about 5 per cent by a procedure which takes about 7 seconds to perform. 

In principle, the method is to analyze the light reflected from the eye as seen in a 
cat’s eye on the road by night or in man through an ophthalmoscope. This re- 
flected light has passed twice through the retina and so must have suffered absorp- 
tion by the retinal pigments traversed. If the incident light falls only upon the 


tiny fovea centralis, the pigments will be cone pigments, and if these are bleached to 
transparency by a strong dazzling flash, the reflected light will return brighter. 


Thus measurement of the increase in intensity of lights of various wave lengths will 
give the difference spectrum of the cone pigments involved. 

Preliminary reports of equipment and results have already appeared,? and some 
C.antitative results upon the kinetics of cone pigments have been published.* It 
v iJ] therefore be unnecessary to repeat here the results presented at the National 
‘cademy, though they show some improvement upon the earlier experiments. 
] it the main conclusions may be summarized. 
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1. Measurements upon the protanope.—This is the type of red-green color blind 
who cannot see the far red spectrum at all. 

a) He exhibits a change in foveal reflectivity upon “bleaching” with a bright 
light and subsequent “regeneration” in the dark, such as would be expected if what 
we were measuring was the density of a visual pigment. 

b) The amount of bleaching produced is related to the energy of the bleaching 
light by the curve expected from photochemical theory. 

c) Bleaching by red or green lights caused identical difference spectra; thus the 
protanope has only one visual pigment (in the red-green range). This would ae- 
count for his monochromacy in this range if the pigment measured were that upon 
which vision depends. 

d) The pigment measured is a visual pigment because red and green lights, 
which appear to the protanope identical, are found to bleach equally. And the 
difference spectrum coincides well with the protanope’s luminosity curve. 

Measurements upon the normal eye. 

a) Adeep red bleaching light, which has no effect at all upon protanope pigments, 
gives, with the normal, a difference spectrum, indicating the presence of a red- 
sensitive pigment which the protanope lacks. 

b) The total bleach in the normal shows the presence of a second photosensitive 
pigment, and this appears identical with that found in the protanope. Names for 
these two human cone pigments are ‘“‘chlorolabe” (= green-taking), with maximum 
absorption at 540 my, and “erythrolabe,” wi h maximum at 590 mu. 

c) Inthe normal there is somewhat less chlorolabe than erythrolabe, but in the 
protanope the amount of chlorolabe is about equal to the sum of all the pigments in 
the normaleye. It therefore seems that the protanope who has the normal number 
of cones is peculiar just in this—that those cones which in normal men contain 
erythrolabe have this pigment replaced by chlorolabe in the protanope. 


1G. Wald, P. K. Brown, and P. H. Smith, /. Gen. Physiol., 38, 623-681, 1955. 
2W. A. H. Rushton, V.P.L. Symposium, Paper I (London: H.M. Stationary Office, 1958). 
W. A. H. Rushton, Proc. N.Y. Acad. Sci., 74, 291-304, 1958. 


COLOR VISION AND THE NATURAL IMAGE. PART I* 
By Epwin H. LAnp 
POLAROLD CORPORATION, CAMBRIDGE, MASSACHUSETTS 


We have come to the conclusion that the classical laws of color mixing conceal 
great basic laws of color vision. There is a discrepancy between the conclusions 
that one would reach on the basis of the standard theory of color mixing and the 
results we obtain in studying total images. Whereas in color-mixing theory the 
wave-lengths of the stimuli and the energy content at each wave-length are signifi- 
cant in determining the sense of color, our experiments show that in images neither 
the wave-length of the stimulus nor the energy at each wave-length determines the 
color. This departure from what we expect on the basis of colorimetry is not a 


small effect but is complete, and we conclude that the factors in color vision hitherto 
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Fig. 1.—Diagram of camera and projector 


regarded as determinative are significant only in a certain special case. What has 
misled us all is the accidental universality of this special case. 

It was scientifically proper to base our earlier hypotheses about the functioning 
of the eye on our knowledge from this special case. Now, however, these recent 
experiments show that the special case is extraordinarily restrictive: The eye will 
see color in situations entirely unpredictable on the basis of older hypotheses. ! 
Thus it may be necessary to build a new hypothetical structure to explain the 
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Fig. 2.—Kodak Wratten filter No. 24. Diagram of filter 
transmittance reprinted, with permission, from the copyrighted 
Kastman Kodak Company publication, Kodak Wratten Filters for 
Scientific and Technical Use. This applies also to Figs. 8, 9, 10, 
and 11. 
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mechanism for seeing color. The purpose of this paper, however, is not to discuss 
internal visual mechanisms, but rather to begin the description of the laws for a 
large general case-—laws describing the relationship of the sense of color in the 
human eye to the structure of total images, laws in which wave-length of stimulus 
and intensity of stimulus will play entirely new roles. 


APPARATUS 


For these experiments two black-and-white photographic positive transparencies 
are made in a split-beam cam- 
era, so that they are taken at 
the same time from the same 
viewpoint through the same 
lens (Fig. 1). One is taken 
through a Wratten No. 24 
filter, passing wave-lengths 
longer than about 585 my (Fig. 
2), and one through a Wratten 
No. 58 filter, passing wave- 
lengths shorter than about 585 
mu. These two photographic 
records will be referred to 
hereafter as the ‘“‘long record” 
and the “short record,” re- 
spectively. The picture is 
exposed so that the gray scale 
appears to have the same den- 
sities in both records. 

The two records are pro 
jected in a double-lens pro 
jector and superimposed ex- 
actly on a screen by means of 
aregistration system. In front 
of each lens is a pair of polar 
izers mounted so that the in- 
tensity of the light in each 


beam can be varied separately, 


In projecting, the long record _Fic. 3.—Blonde girl (Experiments 2, 3, 4, and 5 

is projected ordinarily through Upper image projected with longer wave-length stimulus 
; E Pe ge lower image with shorter wave-length stimulus. This 

a red filter (Wratten No. 24), applies also to Figures 4, 5, 6, 7, 12, and 13. 

and the short record through a 

neutral density filter of about 0.3 (when the polarizers are not being used). The 

half of the projector containing the long record will be referred to hereafter as the 

“long projector,” and the other as the “short projector.””? 


CHROMATIC PHENOMENA 


The experiments demonstrated were the following: 
Kaperiment 1. The double projector is turned on, a filter transmitting red 














“ores 


Frc. 5.—Numbered board of standard objects (Experiments , 10 
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(Wratten No. 24) is placed over the lens of one projector and the polarizers are 
rotated to change the observed color of the screen from full white to full red, travers- 
ing the various shades of pink in between. Thus the screen is seen in the colors 
expected on the basis of classical color theory. 

Experiment 2._-The red filter is removed, and the two components in Figure 3 
are shown separately by occluding first one of the projectors and then the other, 
so that the audience may see the two black-and-white images produced by the 


beam-splitting camera. The differences to be expected in any color-separation 
photographs between the picture taken through the red filter and the picture taken 


through the green filter are pointed out to the audience. 

Experiment 3. With the 
short projector occluded, the 
long record is shown projected 
through the red filter. Then 
the lens of the short projector 
is uncovered so that the black- 
and-white picture is combined 
with the red picture. Instan- 
taneously the portrait: appears 
in full color: blonde hair, pale 
blue eyes, red coat, blue-green 
collar, and strikingly natural 
flesh tones. (While it is not 
the primary purpose of these 
experiments to achieve verisi- 
militude, one is impressed by 
the frequent acceptability of 
this and many of the other 
pictures as lovely color photo- 
graphs. ) 

Experiment 4..-The room 
lights are turned on full, flood- 
ing the sereen with light, 
and it is shown that there 

Groceries (experiment 12 is no visible time lapse be 
tween snapping off the lights 
wnd the appearance of full color in the picture. 

Experiment 5, The room lights are brought up from darkness to a medium level 
at which objects in the room can be clearly seen in color, but at which there is not so 
much light that the picture is washed out entirely. It is shown that the eye sees 
color in the room and color on the screen at the same time and that, furthermore, 
the colors on the screen are not changed by having the surrounding room illuminated 
in white light. This experiment gives the first premonition that multiple color 
universes can coexist side by side, or one within another. 

Leperiment 6. -In this picture (Fig. 4) the audience is to compare the pink states 
on the map with the yellow pencil; and the vellow pencil with the white pencil 
and the white marmalade jar which holds the pencils. It is pointed out that since 
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classically all the colors on the screen would be pink, it is particularly important to 
note the difference between the pinks and the yellows and equally important to 
notice the strong sense of whiteness in the white objects. This is the first demon- 
stration that it is not necessary to use traditional complementary colors to produce 


white, a fundamentally important fact in the image situation. The green-blue color 
of the large book entitled Bartleit’s Familiar Quotations is contrasted with the deep- 
blue stripe on the air-mail letter, as well as to the blue tax stamp and the pale blue 
ocean on the map. The gamut of colors in this picture includes also the brown 


wooden box, the green blotter, the black inkstand, and the red lettering on the 
magazine. It is pointed out that this slide demonstrates dramatically that the 
transition from one stimulus 

(here, either red or white) to 

two (here, red and white) is the 

tremendous step: the full 

gamut of color appears. (The 

addition of a third stimulus 

does enrich the colors some- 

what, but we do not know 

whether what is missing is due 

to inadequacy of technique or 

to the need of a little help from 

a third stimulus. ) 

Experiment 7.—A_ colorful 
picture, notable for the auburn 
red hair and the green 
sweater, is used to show that 
the field that counts in the 
image situation is not de- 
termined by angular subtend. 

The observer finds the sweater 

to be green even when he 

stands so close to the screen 

that the sweater fills a large 

portion of his field of view; 

and he also finds the sweater Fic. 7.— Bowl of fruit (Experiments 13, 14, 15, and 16) 
to be green when the image is 

made small and the observer is at a distance from the screen. 

Experiment 8.—The same picture is shown with the red filter placed on the short 
projector, and with no filter on the long projector, except that the polarizers are 
turned to give the best color. The sweater is now reddish and the hair greenish, 
and the lips an intense blue-green, illustrating the phenomenon of hierarchical 
reversal, usual, to be sure, in any color photography system but particularly im- 
portant here because it demonstrates that the colors are independent of what the 
observer expects, and that the system, in spite of use of the new “arbitrary” pri- 
maries, behaves like any system of color reproduction. 

Experiment 9.—A slide is shown in which there are twelve colored objects photo- 
graphed against a gray background divided into twelve numbered squares (Fig. 5). 
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It is pointed out that, except for the orange on square No. 2, each object could 
appropriately be any color, so that there can be no prior association between the 
name of the object and the name of the color of the object; thus one would not know 
that the telephone is red, the paper stapler gray, or the design on the teacup, on 
square No, 12, dark blue. It is pointed out that this slide had been used to in- 
terrogate naive observers, and that, except for the color-blind, marked consistency 
was found in their reporting of the color names. 

Experiment 10.--Using the same standard-objects slide as for Experiment 9, 
a pair of polarizers is placed in position in front of each projector, the relative 
brightness is varied from the position in which only the red image shows, to the 
position where just the white image shows. The audience is asked to indicate the 
two extremes at which color first appears; that is, with only the red image on the 
screen, the white is slowly added, and the colors at the point where they first appear 
are noted. Then, with only the white image on the screen, the red is slowly 
added, and as the colors appear they are again noted. It is pointed out that when 
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Fic. 8.— Kodak Wratten filter No. 70 


color does appear at both extremes it appears at the same time for all the objects 
on the slide. It is demonstrated that the ratio of stimuli can be varied enormously 

with some pairs of stimuli by ratios as high as 100 to 1—without altering the 
color pame given to each one of the objects. 

Experiment 11.—-In this experiment a duplicate of the short record is super- 
imposed on the short record already in the projector, thus doubling its contrast. 
A single long record is used, so that it has the original contrast of the previous experi- 
ments. Thus for the short record in places where, in the previous experiments, a 
fifth of the light was transmitted, the slide now transmits a twenty-fifth; in places 
where a tenth was transmitted, the slide now transmits a hundredth; in places 
where all the light was transmitted, the slide now transmits all the light. Since no 
change was made in the long record, the same amount of light is transmitted by it 
as Was transmitted in the earlier experiments. It follows, therefore, that for each 
object on the slide there is now a new ratio of long to short stimulus, and it would be 
expected classically that each object would now be a new color. Actually, the 
objects retain their original color names, demonstrating the stability of the colors 
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of an image with respect to the contrast of its component images, and illustrating 
further that color is independent of the quantity of the stimulus defined locally. 
It is pointed out that we must not let the casualness of the statement that we have 
doubled contrast in one image without doubling it in the other obscure in our minds 
the amazing fact that utterly new pairs of ratios give the same colors in the same 
places. 

Experiments 12, 13, 14, and 15 are designed to demonstrate further that color at a 
point in an image is independent of the wave-length composition of the radiation 
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coming from that point. While with changes of filters there is a slight over-all 
wash of color in extended background areas, objects retain their original color 
identity. (Many of the audience sitting off the axis of the projector are unable to 
witness the following four experiments because of the extremely low light level re- 
sulting from the dark filters. ) 

Experiment 12.—A slide showing multicolored groceries on a table top (Fig. 6) 
is shown with Wratten No. 24 and no filter on the one hand, and on the other hand 
with Wratten Nos. 24 and 58. 
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Experiment 13.-A bowl of fruit on a table top (Fig. 7) is shown with Wratten 
No. 70 (Fig. 8) and Wratten No. 73 (Fig. 9). 

Experiment 14.--The slide in Experiment 13 is shown with Wratten No. 73 and 
Wratten No. 49B (Fig. 10). 

Experiment 15. The same slide is shown with Wratten No. 47B (Fig. 11) and no 
filter. 

Experiment 16.—It is demonstrated to the audience that the filter Wratten No. 73 
which passes a narrow band of yellow may be used as a shorter member of a wave 
length pair, and as a longer member of a wave-length pair, depending on what is 
picked to go with it. (Note: Experiments 13 and 14.) 


ACHROMATIC PHENOMENA 


Experiment 17.—Three identical black-and-white transparencies are used, two in 
superposition on each other in one projector, one in the other projector. These 
are projected with red and green filters, and the question arises: What colors should 
one expect? One expects color for two reasons, either classically because of the 
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varieties of ratios of red to green or in terms of our new approach because the 
numerical value of many of the ratios are equal to ratios which existed in the 
image which showed color in Experiment 9. Actually, one finds that the whole 
slide is uniform in color and nearly colorless, even though many of the ratios of 
stimuli for the various objects are the very ratios that existed when the slide ap- 
peared in full color. 

Experiment 18.—The short record is placed in the short projector and no record 
is placed in the long projector. One sees the wash of red light over the image from 
the short projector. The polarizer on the long projector is varied over its range. 
In spite of the fact that all possible ratios of stimuli have been on the screen, there 
appears only a uniform wash of the colors in their classical mixtures. 

Experiment 19.—Another class of achromatic phenomena is introduced by this 
experiment in which a photographic step-wedge, running horizontally, is placed in 
one projector, and a step-wedge, running vertically, is placed in the other projector. 
Thus, in superposition on the screen, the wedges create 256 possible ratios of the 
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two stimuli, arranged, however, in a geographically systematic way. The colors 
mix only in a classical sense, red and white giving a vague pink, and red and green 
showing some red, some green, and some yellow. 

Still another class of nearly colorless situations is exemplified by Experiments 20 
and 21. 

Experiment 20.—The map slide (Fig. 4) is projected with the long and short 
record far out of registration. The gamut of colors is strongly restricted. While 
going beyond the classically expected colors, there is only the smallest suggestion 
of the new colors which here are few and unsaturated. The images are slowly 
brought into registration. At the instant of perfect registration, the full gamut of 
color snaps into appearance. 

Experiment 21.—A picture of one scene is placed in the long projector, and a 
picture of a different scene is placed in the short projector. The limitations on 
color are similar to those in Experiment 20: there is almost none. 
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A NEW CO-ORDINATE SYSTEM 


A new co-ordinate system is shown which describes color in images and which 
predicts the achromatie situations (Fig. 12). A remarkable thing about the new 
co-ordinate system is that it is physically dimensionless, involving a ratio of ratios 
and holding constant an arbitrary pair of wave-lengths. In this co-ordinate system 
the percentage of available short-wave stimulus is plotted against the percentage of 
available long-wave stimulus.* The first message of the co-ordinate system is that 
color in an image is indeed independent of the over-all flux in the individual com- 
ponent images and, second, that the color at a point in an image is independent of 
the wave-lengths of the radiation at that point. This is true of what is seen with 
a large variety of pairs of wave-lengths and relative brightnesses. Part II of this 
paper will show that the gamut of colors is restricted when some pairs of wave- 
lengths are used, and that there are levels of relative energies at which there are 
significant color shifts. However, even in these situations, the color is not at all 
the classical function of wave-length and relative energy. 
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We now study our twenty-one experiments in terms of the co-ordinate system. 
Slides such as those in Figures 3, 4, 5, 6, 7, and 17 are measured as indicated above, 
and the results are shown in Figure 13. Here we see that in this new dimensionless 
co-ordinate system the grays fall on a straight line at an angle of 45° to the axes, 
warm colors being to one side of the gray and cools to the other. 

Doubling the contrast of the image, as in Experiment 11, rotates the gray line, 
which carries with it the family of warm and cool colors, disposed as they were 
originally on either side of the gray line; in Figure 14 imaginary lines have been 
drawn through some of the colored areas, the unbroken lines representing points in 
an image of normal contrast, the dotted lines indicating the new position of those 
points when the contrast of one of the images has been doubled. The significant 
fact is that the colors themselves remain constant for the observer, in spite of this 
translation in the co-ordinate system, as if the eye were indifferent to the change 
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in position of the array as a whole as long as the relative positions of the color 
points in the array remain unchanged. 

Achromatic Situations..-When the points in the images in Experiments 17 and 18 
are plotted on the diagram they are found to fall on straight lines (Fig. 15), suggest- 
ing that this is the criterion for whether or not an image is achromatic, and possibly 
suggesting further that if the points measured in an image fall on any smooth 
continuum in the diagram, such an image may also be achromatic. For three 
identical images in Experiment 17 all the points fall on a line (which is in the same 
position as the gray line in the double-contrast image of Experiment 11). When a 
single image in one projector is combined with a wash of light in another, as in 
Experiment 18, the measured points will fall on the straight line parallel to one of the 
axes. As the observer’s eye moves over the field of the crossed wedges in Experi- 
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ment 19, in any direction, the points traversed, when measured, will fall on simple 
smooth curves on the diagram of the co-ordinate system, 

If the wedges were continuous instead of being step-wedges, then one might 
think of the crossed wedges as being the continuum which is the co-ordinate dia- 
gram as a whole. 

The last class of situations, which are nearly colorless, includes the out-of-register 
images of Experiment 20 and the two different images of Experiment 21. If one 
thinks of the fusion of the components of the colored image on the retina as being 


analogous to the central fusion of the two members of the stereoscopic pair on the 


two different retinas, then it seems plausible that the sense of color might be greatly 


Fic. 16.—Spectrum of natural objects, negative and positive 


diminished when the component images on one retina are out of register, just as 
the three-dimensional sense will be greatly diminished, indeed destroyed, if the 
members of the stereoscopic pair do not lie in “registry” on the corresponding points 
of the two retinas. A second way of thinking, which will not be discussed in this 
paper, has to do with the relationship between the way the eye moves over the 
geography of the image and the trace of the point on the geometry of our diagram. 
Experiment 22. -To test the hypothesis that continua in general, and straight lines 
in particular, have a special significance in our co-ordinate system, one may ask the 
following question: If the measured points in the image fall on a straight line at 
right angles to the 45° line in the co-ordinate system, what colors will appear? 
Now if a negative image is placed in one projector and the positive from that 
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negative is placed in the other projector (Fig. 16), then the measured points will 
fall on a line at right angles to the 45° line. The red filter is placed on the long 
projector and no filter is placed on the short projector; the negative and positive 
images are registered carefully. Is the co-ordinate system a good prophet? All 
ratios of stimuli are present in these images which are geometrically identical and 
in register. Will this be an image in good color, or will it be essentially colorless? 
The answer is that the co-ordinate system 7s a good prophet. There are only reds 
and whites, and a little pinkness. This is further support for the hypothesis that 
stimuli whose measurements are disposed on a straight line in our co-ordinate system 
do not “mix” to give color. (The control picture, the two positives taken and 


Fic. 17—Spectrum of natural objects, control. 


shown in our new “standard” way [Fig. 17], is vividly colored.) This negative- 


positive pair promises to be one of the great clues to the nature of color vision. — It 


will be discussed in the next part of this paper. 

Because of the limitations of flexibility imposed by using filters, a dual mono- 
chromator has been built since the time of the symposium. The work with this 
instrument has supported the experiments shown at the symposium and has enabled 
us to arrive at a meticulous description of the separation of wave-lengths, as a fune- 
tion of wave-length, that is required to see the gamut of colors for experiments like 
those described above. This instrument and the detailed results obtained with it 
will be described in Part If. There will also be described new apparatus demon- 
strated at the meeting of the Optical Society of America on October 9, 1958, which 





Vou. 45, 1959 PHYSIOLOGY: E. H. LAND 129 


enables an audience to see that the same monochromatic stimulus can be either the 
long or short member of a pair of stimuli. The demonstrations given at the 
Physics Colloquium, Harvard University, on November 17 and 18, 1958, showing 
how the eye functions to see full color in a universe involving only a narrow wave- 
length band (‘“Johnny’s Magic Bookcase’) will be described, and the generalization 


will be drawn that our everyday world is just a special case among these universes. 


Part IIT will continue the discussion of the nature of the field in the image situa- 
tion and the relationship of the spectrum to color vision. 


* This paper will be in several parts, the first of which is this symposium. Part I then contains 
the experiments up to the time of the symposium, April 28, 1958. 


' The work in this paper should not be confused with experiments such as those of Hess, which 
describe how a stimulus of a particular wave-length looks differently colored when viewed after 
the eye has been exposed to another wave-length for a significant period of time (Arch. Ophthalmol., 
36, I [1890] ). Nor should it be confused with those theories which explain the phenomenon which 
Hess measured, and other allied phenomena, as the result of fatigue of the eye to various states of 
chromatic adaptation. 

The colors described in our experiments appear immediately, and do not alter appreciably with 
time, and we do not seek to explain them as effects of adaptation. To the best of our knowledge 
the gamut of color is much larger than that noticed by Hess or predicted by any theory concerning 
chromatic adaptation; our newly discovered achromatic situations are not predicted by any pre- 
vious theories; nor is there any basis in those theories for our dimensionless co-ordinate system 
(see below). Furthermore it seems that none of the previous work in any kind of color matching 
or system of complementary colored shadows (Goethe) can be the basis for the observed fact that 
the colors in our images, arising from unorthodox stimuli, correspond in hierarchical order to the 
colors that the observer would have seen had he been at the criginal site. 

? Karlier studies in two-color projection were done by Fox and Hickey who found that moving 
pictures which showed some color could be taken by alternately photographing through no filter 
and through a filter which transmitted red light, and by projecting correspondingly (British Patent 
No. 636, issued July, 1914, to William Francis Fox and Wm. H. Hickey, /mprovements in Kine- 
matographic Apparatus). Some years later Bernardi used red and green filters for taking movies 
and red and white for projecting (Adrian Cornwell-Clyne, Colour Cinematography {London, 1951}, 
p. 261). 

3This means that one measures the amount of light on the screen for the long wave-length 
image in the brightest place and takes that as 100; then one measures the amount of light for 
that long record alone at another place in the image and makes the percentage. Then that per- 
centage is plotted against the percentage of available light for the short record alone at the same 
point on the screen. The ratio of these two percentages describes the color at that point. 
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